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PREFACE 


III tile present work our mean efforts have been to 
study the static properties of light nuclei (A ^.40) which are 
obtained from Hartree-Pock calculations with realistic inter- 
actions. 'Tiiis work mainly consists of five chapters. Chapter I 
is mainly an introduction giving a brief account of the indep- 
endent particle model of the nucleus, and the Hartree-Pock and 
Brueclcner-Haa^tree-Pock approximations of solving many-body 
problem. Several Hartree-Pock and Brueckner-Hartree-Pock 
calculations that have been done in the past are also summarized 
and the problems taken up for the present work have been 
mentioned. Chapter II describes the nucleon-nucleon inter- 
actions which can be used in nucleaa structure studies. We 
have chosen the following two interactions for our purpose; 
effective interaction of Shakin, Waglimare and Hull derived 
from Yale potential by using a unitaa-y -model-operator approach 
and the Sussex interaction of Elliott et. al. obtained directly 
from experimental phase shift data by using an auxiliary potent- 
ial method. Chapter III contains the Hartree-Pock theory in 
detail and the matrix method of solving the Hartree-Pock 
equations. Ihe formulae used in the calculations of various 
ground state properties are also given. Two important cor- 
rections due to the coulomb energy and centre-of-mass motion 

of the nucleus are estimated. Both these corrections are 
included in a self-consistent manner in our calculations. 



unlike ilie ea.nliex Haxtree— Pock calculations of various other 

authors. These corrections are estimated, by those authors 

externally after solving the Hartree-Pock equations. Single 

oscillator approximation limit to the Hartree-Pock solutions 

has also ‘been discussed. In Chapter IV we present the results 

of spherical Hartree-Pock calculations for the nuclei ^He, ^Be, 
12 ^ , 40 ^ 

0, 0 and Ca. The properties calculated are the binding 

energies, the single particle energies, r.m.s. radius and 
charge radius, matter and charge densities, spin-orbit split- 
tings , etc. It is 'found that the binding energy results, 
obtained with effective Yale interaction for all the nuclei 
considered are in satisfactory agreement with the correspond- 
ing experimental values if all the second-order terms of the 
Yale potential are accounted for. Sussex matrix elements givg" 

better results for the binding energies in the first order 

40 

compared to the first order in Yale potential. Por Ga, 

however, Sussex matrix elements give binding energy in excess 

of the experimental value. Rest of the properties calculated 

are in fair agreement with the experimental data. The results 

obtained in the single oscillator configuration approximation 

to Hartree-Pock wave functions are also presented. Single 

particle v/ave functions and self-consistent single particle 

potentials obtained in the Hartree-Pock calculations with 

effective Yale interaction have been studied in detail. The 

16 

Hartree-Pock single particle wave functions for 0 have been 
compared with the conventional harmonic oscillator wave 
functions. The general nature of the wave functions is the 



same for all tiie nuclei, liie Hartrce-Pock single particle 
potential is non-local and state-dependent, file static limit 
to this non-local potential has been discussed. 

finally, in Chapter V wc have studied the energy 
IS 18 10 

levels of y F and ^Be in the Hartree-Pock representation. 
It is assuiaod that the nucleons move independently of each 
other in a Hartree-Pock self-consistent field and the low- 
lying states of these nuclei may arise due to the excitations 
of only the last two neutrons to the higher states. The single 
particle energies and wave functions describing the behaviour 
of these two nucleons are obtained from Hartree-Pock calcul- 
ations. The nucleon-nucleon potentials used are the effective 
Yale interaction of Siiakin et, al. and Sussex interaction. It 
is seen that the HP spectra are compressed compared to experim- 
ental levels in all the cases. Probable reasons for this 
discrepancy have been discussed. 


(xii) 



CHAPTER I 


IITRODUGTIOI 

I . 1 Independent Particle Model 

The study of the nucleus is a many-body problem which 
cannot be solved ene.ctly . Instead, studies of complex nuclei 
have been carried out by various phenomenological models. There 
are several nuclear models which can be broadly divided into 
tvro classes. First class consists of so-called ‘strong inter- 
• action models’ which treats the nucleus as an assemblage of 
closely coupled particles, e.g., liquid drop model, Bohr and 
Mottelson model, etc. These models have been successful in 
explaining many collective properties of nuclei but we shall 
not discuss them here. The second class is of ‘independent 
particle models'* which assumes that the nucleons move independ- 
ently of each other in an average potential. Recent evidence 
suggests that the independent particle model provides a better 
description for low-energy phenomena though the coupling of 
nucleons by virtue of their mutual interaction can by no means 
be neglected. The most successful of these independent particle 
models are the shell model, the optical model and the Rilsson 
model. The uniform potential is spherical for shell model, 
complex for optical model and spheroidal for Nilsson model. 

Shell model, proposed by Mayer in 1948 and developed since 
then, has been very successful in explaining some regularities 
in the binding energies of nuclei, spin and parities of some 
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of tile low-lying nuclear leizels, certain multipole transitions, 
P-decay, etc* Ihe optical model was primarily designed to 
explain the nuclear reaction phenomena and it successfully 
accounts for the scattering cross-section data of low and medium 
energy reactions. Nilsson model can predict well the ground 
state spins, moments of inertia, energies of excited state ro- 
tational hands, equilihrium deformations, and Tih rational pro- 
perties of deformed nuclei , which are not described well on the 
basis of simple shell model. In all these independent particle 
models the nucleons fill up the successive lowest quantum levels 
of the mean potential well and move in definite orbits. Olhis 
picture provides the qualitative features of many nuclear pheno- j 
mena, lo get the quantitative agrsement with experiments the { 
residual interaction between the particles is introduced. It ! 
mixes up the orbits to some extent but the basic character of 
the model is essentially the independent particle motion. 

With increasing addition in our knowledge of nucleon- 
nucleon force from the study of two, three and four nucleons 
systems, it was thought that there ^ould be a way of correlat- 
ing the average potential field of the nucleus with- the nucleon- 
nucleon interaction. The strong objection to this was the fol- 
lowing. Most of the realistic nucleon-nucleon potentials 
which fit the two-body data contain a repulsive core. Such a 
repulsive core will give rise to strong short-range correlatio.ns, 
which are not described well by an independent particle model. 
However, there are many experimental evidenoes in favour of 



irLdependent particle motion, e.g* the nuclear reaction experi- 
ments, in which an incident particle loiooks out a nucleon from 
target nuclei, show that single particle structure does prevail 
for the nucleons inside the nucleus. Besides, elastic scatter- 
ing sho-we re sonance s which are explained in terms of optical 
potential model. Thus, the fact that independent particle model 
is successful in explaining a large number of experimental data 
implies that short-range correlations do not play an important 
role in many nuclear properties and the nucleons inside the 
nucleus behave- to a certain extent as independent particles in 
a common potential. Recent advances in the theory of nuclear 
matter by Moszkowskl and Scott (Moi 60) and Oomes et, al. 

(Gos 58) have explained the mechanisms which reduce the effect 
of strong interaction on the nuclear motion. In particular, the 
essential features of IJoszkowski- Scott method could be described 
as follov;s: 

(i) The scattering of two nucleons inside a nuclear medium 
is characteristically different from the collision of free 
nucleons. The reason for this is that inside the nuclear 
medium Pauli principle forbids scattering to states which are 
already occupied and thus, the number of available final states 
is limited. The nucleons far apart can interact through low / 
Pourier components of the force. Such interactions are gene- 
rally forbidden inside the nucleus since 1±Le momentum transfer- 
is insufficient to raise the particles above the Permi sea. 

The short-range part of the potential eorcesponds to. high 
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momentum transfer and scatters the particles to unfilled states. 
Thus, "beyond a certain healing distance the relative wave 
function of two nucleons approaches that of noninter acting 
particles- This is in sharp contrast to scattering of free 
nucleons where the wave function never heals up and is phase 
shifted relative to free particle wave function even at infi- 
nity. The healing distance is of the order of 15* and the ave- 
rage separation of the nucleons in nuclear matter is approzi- ; 
mately l.SSF. Consequently, the effect of two-hody correlations 
would he small. Besides, the fact that healing distance is 
smaller than the average separation implies that the proh ability 
that more than two particles will simultaneously he at distances 
smaller than the healing distance would he small and three-hody | 
correlations would he unimportant. | 

(ii) The effect of the attractive part of the potential on | 
the nuclear phase shifts is opposite to that of the repulsive 
part- Thu^ a part of the attractive potential can he comhined^ 
with the repulsive core so that the two together give zero phase 

If 

shifts. Consequently, 'iiie phase shifts due to the strong nucleon- 
nucleon interaction are the sane as that due to weak attractive I 
long-range part of the force heyond the healing or cut-off 
distance, . 

As a consequence of these two effects particles move 
almost independently with very weak interactions. This implies 



that the nuclear Hamiltonian H of a system of A nucleons can 
he written as 


H 





( 1 - 1 ) 


where is the independent particle Hamiltonian and v is the 
residual interaction, llie average potential fields which each 
nucleon feelsj must he some sort of average of the interactions 
it has with all the other nucleons. The residual interaction 
between two nucleons must he the interaction left over after 
this average effect is removed. A systematic way of determining 
Hqs starting from fundamental nucleon-nucleon interaction, is 
the self-consistent Hartree-Fook (HF) method. The ground state 
of the nucleus is represented, in the first approximation, hy 
a nondegenerate eigenstate (D^ of II i.e. 



( 1 . 2 ) 


where is the lov/est eigenvalue of H^. We shall briefly 
describe the HE method of determining^^ and E^. 


1'2 Ha rtree-Eo c k Met hod 

The Hartree-Eock (EE) method consists in assuming tha.t 
the wave funGtion^^ of the nucleus is an anti symmetrized pro- 
duct of single particle wave functions. To determine the best 
single particle wave functions there are many eqtrl valent criteria-' 
Eor example, one can minimize the total energy or one can require 
self-consistency between the single particle potential field of 
H and the average potential seen by the particle as a result 
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of its interaction mtii aXl the other nnoleons. The single 
particle wave functions v/hich minimize the total energy are the 
solutions of following Hf equations 

<a|T|p)> + <aiUlp> = , (1.3) 

where a and jB specify the single particle states, T is the 
Icinetic energy and h bhe HF single particle potential which is 
related to the nucleon-nucleon potential V hy the following 
e quati on, 

A 

<afujp^ = ^ <aYlv|P'j^ . (1.4) 

Y=1 

Thus, the HF single particle potential U is calculated in teims 
of Y hy a sum over occupied single particle states, hut the lattei 
cannot he obtained without a prior knowledge of U. This is the 
HF self-consistency prohlem and is solved hy an iteration 
interpolation procedure, starting from a trial guess and improv- 
ing the wave functions until self-consistency is achieved. The 
HF field may he spherical or deformed depending upon the nuclei 
studied. Thus, for example, the HF field is spherical for the 
douhly closed shell nuclei Sle, ^^0, ^*^Ga and ^^^Ph. The ground 
statei^Q is a Slater determinant wherein all the single particle 
states are occupied up to Fermi level e|,. For . deformed nuclei 
single particle states are described in a deformed HF field. 

The ground stateQ^ is again a Slater determinant in which all 
the deformed orbitals are occupied up to Fermi level. 
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Tlie HF method not only descrihes the independent parti- 
cle picture _of the nucleus hut can also he used in the discus- 
sion of collective motion of the nucleons, e.g., the vibration 
and rotation of the nucleons as a whole can he explained hy let- 
ting the self-consistent HIP field become time -depen dent. Besides, 
the pairing correlations of the two-body force can he taien into 
account in an extended HB' theory, known as Hartree-Fock- 
Bogoliuhov-Valatin theory, Ihis method extends the variational 
principle to a class of trial wave functions of independent 
excitations of quasi pairtioles which include pairing correla- 
tions. In this case the total Hamiltonian H can he written as 

H = + v*^ , (1.5) 

OP 

where is the independent quasi particle Hamiltonian and 
OP 

V “ is the residual interaction between these quasi particles. 

HP -theory has been successfully applied in atomic and 
solid state physios. In fact, the idea of self-consistent field 
was first intutively introducod by Hartree and was successfully 
applied to electronic structure of atoms. Later this was 
extended by Pock who included the exchange effects and estab- 
lished the. HP equations using a variational method. Ihus,the 
HP method in nuclear physics is essentially borrowed from 
atomic physics, but the application of this method to nuclear 
structure calculations is much more tedious than in the case 
of atoms because of the followiiig reasons: 
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(i) Tile micleon-xiu Cleon interaction y is not knovm 
uniquely. 

(ii) Most of the realistic nude on- nucleon potentials which 
fit the two-hody data contain a hard core. The matrix element 
of such a potential, in a basis of independent particle wave 
functions is infinite, 

(iii) Complete HP calculations are very time-consuming and 
cannot be cariied out without high-speed computors. 

■ I 

, Inspite of all these difficulties several spherical 
and deformed HP calculations have been done during the last 
few years. A detailed survey of all these oaloulations is not 
possible here, Por an excellent review of defomed HP calcu- 
lations we refer to an article by Bipka (Bia 68),' Since our 

j 

main interest lies in the spherical HP calculations, we. shall ; 

briefly review only such calculations tlncu have been carried 
out in the past. Spherical HP calculations by Davies et. al. 

(Das 66) have been carried out using simple velocity- dependent 
force acting only in relative s-state. They have studied the 
nuclei ^^0, ^^0, ^®Si, ^®Ca and ^'^Zr and obtained qualitative 
agreement for various calculated ground state properties with 
experiment. In the calculations of Muthukrishnan et. al. 

(Mun 65), Yamaguchi potential has been used and the results 
obtained for various nuclei are ve3^ similar to those of Davies 
et, al. (Das 66). Abgrall and Monsonego (ibl 66) have used a 
central force with a nonsaturating exchange mixture, Alas, 



there are oalculations done Toy Kileger et. al. (Krr 66), whose 
effective interaction includes three-body contact term propor- 
tional to cube of density. In their later calculations, Davies 
et. al. (Tan 68) have eraployed a more complete two-body poten- 
tial including both even and odd noncentral potentials and . 
giving a fair fit to two-body data. They have studied in detail 
a number of closed shell nuclei in the range 16 4 208. 

Although^in some of these calculations the calculated HF proper- 
ties are in reasonable agreement with experiments, -the mayor 
drawback of these calculations is that they were carried out 
with nucleon-nucleon potentials especially designed for use in 
HF calculations. The parameters of these potentials were . 
adjusted to fit the nuclear matter data so that they gave a 
fair fit for the binding energies of various nuclei (particu- 
larly heavy nuclei) with the experimental data, but they did 
not fit the two-body data. In fact these effective potentials 
do not contain any feature of the realistic nucleon-nucleon 
potentials, such as hard core or nonlocality, one-pi on-exohango 
behaviour at large distances, etc. First time a realistic poten- 
tial (Tabahin potential) in the HF oalculations has been used 
by Kerman and ec -workers (Ken 66), These authors have studied 
the nuclei ^^0 and in considerable detail and find that 

Tabakin potential fails to give adequate binding in first-order 
and gives too large a value for the spin-orbit splittings. 

Later on,Pal et, al. (Pal 66) have improved the agreement for 
the binding energy by calculating the contribution of second 
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order term. In all the H]? calculations descrihed above the 
nucleon-nuGleon interaction was smooth and did not contain, any 
hard core singularity. Ihe first theoretical foundation for 
the use of singular potential in Hf calculations has been laid 
by Shakin et. al. (Shn 66, 67a). Iheir theory is based on the 
unit ary-model- operator approach of Villars (?is 65) and is 
applicable to any general singular interaction. Ihe fact that 
short-range correlations do not play an important role in many 
nuclear properties implies that independent particle EP wave 
functions are significant aJid can be obtained from an effective 
interaction from which the effect of hard core has been Removed. 
This was the philosophy adopted by Shakin et. al. in dedijicing 
a smooth effective interaction for use in HP calculations. 

Thus, all the problems from (i ) to (iii) have now largely been 
solved. There is another approximation method proposed oy 
Brueckner and oo—workers for dealing with the problem of hard 
cores. We describe it very briefly below. 


1.3 Brueckner Theory 

In the Brueckner theory motion of a pair of nucleons 
is described by means' of a correlated wave function^^pC^, rg) , 
which vanishes inside the core region [^-r^ | r^ and asympt- 
otically approaches the unperturbed wave funotion^^p: 

( 1 . 6 ) 
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is the solution of a modified form of two-nuoleon Schrodinger 
eq[Uation Icnown as Bruecloier-Bethe-Goldstone equa'^ion (Bee 57), 
aXLd is related to (} througii a reaction matrix (K-matrix) defined 
as 


K(b 




Y' 




(1-7) 


Because of the strong short-range correlations, generated hy 
the two-nucleon potential, it is assumed that the motion of the 
pain will not he affected hy the other nucleons through any 
multipartiole clustering, the other nucleons only modify the 
energy of the interacting pair through HB potential, and Pauli 
principle forbids scattering to states which are occupied hy 
other nucleons. This is known as independent pair approximation. 
The interaction between the nucleons inside the nuclear medium 
is determined hy a reaction-matrix wiiich differs from the two- 
body potential over the short-range of correlation distances in 
the nuclear wave functions. The defining equation for K-matrix 
is 

K = V + 7 -S K , (1.8) 

Q is the Pauli operator which takes care of exclusion principle 
and e is the energy denominator which takes into account the 
binding effects of average nuclear field. In a basis of indep- 
endent particle wave functions 9 Js together with associated 
eigenvalues the matrix elements of E oan.be written as 
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^a|3,Y6 ^a~\'^ ^^0t,y6 * (1*9) 

< 3 % 

The sum over a and i is to he carried out over empty states only. 
The single particle energies ocnuxcihg in equation (1.9) are- 

~ ^ ^^^~a.y ,a.y~^oiy ,ya^ * ( 1 - 10 ) 

y(occ) 

The second term in equation (1,10) defines the single particle 
potential U, i.e. 

= Z1 ^^ay,ay-\y,ya> ■ 
y(ooc) 

In the limit of a weak interaction IC can be replaced by the 
potential V and the single particle potential defined by 
equation (1.11) is the HF potential- Starting from a given 
set of eigenstates 9 ^ (n means n"^^ iterate) the new set 
obtained from the single particle potential will in general 
not agree with the input. The problem of obtaining agreement 
is- the self-consistency problem of HF, In the Brueckner theory , 
a new self-consistency problem arises since a change in the 
representation 92 _>(p^'''^ not only changes , the single particle 
potential, but also through the change in the energy spectrum 
changes the reaction-matrix. This in turn reacts back on the 
single particle potential and eigenfunctions. This is the 
Brueckner self-consistency problem of obtaining self-consistent 
energy spectrum and reaction-matrix. In Infinite nuclear 



matter the eigenstates are plane waves and only Brueokner self- 
consistency problem exists. In finite nuclei Brueokner and HP 
self-consistency problems are coupled together and exact 
Bnueckner-Hartree— Pock (BHP) calculations "become quite compl- 
icated. Besides, in finite nuclei exact treatment of Pauli 
operator causes some difficulties and various approximation 
methods are used to calculate the reaction-matrix. Many authors 
have Bypassed the HP self-consistency pro"blem hy using pure 
oscillator wave functions ("yog 67, 68s Kor 66, 6?: Mar 66). Such 
calculations encounter only Brueokner self-consistenoy pro"blem 
and are called Brueokner- type calculations. In his calculations 
Wong (Wog 67) has studied qualitatively the ground state prop- 
erties of ^*^Ca, ^^0 and "^He using a num"ber of hard core potent- 
ials. It is found that realistic potentials under"bind these 
nuclei "by 2-4 Me"7 per particle, a result similar to that 
o"btained in pure HP caloula,tions with realistic potentials 
(Ken 66: Shn 671). Some BHP calculations for finite nuclei 
have also "been reported in the literature. Pirst few BHP 
calculations (Brr 58, 61s Man 65: Kor 65) have "been done in 
coordinate representation, "but none of them are exact. Calcul- 
ations "by Brueokner and coworkors have "been done in the local 
density approximation. In tlsis approximation, the reaction- 
matrix elements for finite nuclei are taken to he simple 
averages over matrix elements for nuclear matter, calculated for 
the density at the centre- of-mass position of the two inter- 
acting nucleons. In the calculations carried out with Gammel- 
Thaler potential, light nuclei are found to he underhound hy 
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2-4 MeT per particle but discrepancy d.ecreases with increasing 
A. However, the equilibrium nuclear radius becomes too small 
for heavy nuclei. - The results for ^^0, '^^Ca, and ^®^Pb 

are summarized in references (3rr 61) and (Man 63). BHH calcul- 
ations in an oscillator ba,sis have been recently reported by 
Davies et. al. (Das 69). As in the pure HD calculations, the 
single particle wave functions are expanded in terms of harmonic 
oscillator wave functions and the calculations have been carried 

out with Hamada-Johnston potential for the nuclei ^^0, "^^Ga, 

48 2 08 

Ga and Pb. Their main conclusions are that nuclei are too 
small and Underbound. 

We thus see that irP is very much alive even when the 
nucleon-nucleon interaction is sn3.gular, if not the simple HP 
something more sophisticated like BHP is valid. In fact, it has 
been shown (Grg 68) that the effective interaction of Shakin 
et. al., obtained by the method of canonical transformations, 
can also be expressed in terms of K-matrix. However, in the HP 
calculations carried out vd-th this effective interaction 
( Shn 67b) there is no Bruec'cner-self-consistency involved since 
in the calculation of second-order terms in the effective inter- 
action these authors have used an aver, age value for the energy 
denominators and angle averaged Pauli oper.ator Q. Por the 
purpose of present work we come back to simple HP calculations 
which do not involve any Bruckner- self-consistency problem. 





I *4- Problems for Present Work 

The following problems have been undertaken in the 
present works 

(i) Spherical HP caLoulations with effective Yale inter- 
action have been carried out for ^e, ^Be, ^^0 and ^^0. This 
interaction has been widely used in several spherical and 
deformed HP calculations of vo,rious nuclei. Spherical HP 
calculations for ^^0 and ^^Ga. with this interaction have been 
done by Shakin et» al, ( Shn 67b), Present work is on extension 
of their caloulations to other lighter nuclei, ¥e have also 
included coulomb and centre- of-mass corrections in our calcul- 
ations in a self-consistent .manner whereas in the calcul- 

ations of Shcokin et. al. the effect of these corrections on 
the total ground state energy was estimated externally. The 
effect of these corrections on various ground state properties 
ha's also been studied in detail. The results are presented 
in Section. IV, 2 of Chapter IV, 

(ii) Elliott et, ol. at Sussex (Elt 68a) have recently 
piihLisihed relative matrix elements of the nucleon-nucleon inter- 
action in an oscillator basis obtained directly from experim- 
ental phase shift data. As a test of their matrix elements 
Elliott et. al. (Elt 67) have calculated the spin-orbit split- 
tings of light nuclei (A$41) and obtained reasonable values. 

A n-umber of deuteron (Elt 68b) and triton and alpha particle 
properties (Jon 69) have also been calculated and satisfactory 



agreement vd-th experiments lias lieen olitained. The success of 
those calculations suggests that Sussex matrix elements are 
reliahle and can he used in the calculations of more complicated 
nucl-ei. We have used these matrix elements in the spherical 
HI calculations for the nuclei ^He, ^^0, and The 

results are given in Section IY,3 of Chapter lY. A detailed 
comp axis© a of the Sussex matrix elements results with those 
of several other oaloula-tions, particularly with those of poten- 
ticols obtained from a fit to experimental phase shifts, has been 
made in Section IY,4 of Ohaptcr lY. Since the HY formalism is 
well-known to give a good account of the ground state properties, 
our calculations also serve as a test of the saturation prop-, 
erties of Sussex matrix elements. 

(iii) We have calculated the low-lying excited stales of 
1 0 1 R 1 R 

Be, 0 and 1 using HF wave functions. As in the pure shell 
model calculations, it is assumed that low-lying states of these 
nuclei arise due to the interaction of last two nucleons and 
remaining nucleons ane treded as inert, Ihe single particle 
parameters needed in the calculolion a.re obtained from self- 
consistent HP calculations. 

Several shell model spectroscopy calculations for. the 
energy spectra of these nuclei have been reported in the liteiv 
ature. In most, of the earlier calculations the residual inter- 
action between the last two nucleons was taken to be purely : 
phenomenological, e.g. some simple exchange force with Gaussian 
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or Yukav/a radial dependence. The adjustable parameters of the 
potential were determined from best fit to experimental data on 
one or more nuclei. Such calculations have been done by Waghmare 
( ¥ae 64 ) > , Elliott and Flowers (Elt 57)? and many others. 

The residual interaction so obtained depends upon the config-* 
urations taken into account and the energy levels which are 
included in the calculation and cannot be compared with any 
standard realistic potentials. A realistic interaation in the ; 
energy level calciilations for A = 18 nuclei was first used by 
Kuo and Brown (Kuo 66). These authors? starting from Hamada- 
Johnston potential* have obtained an effective interaction which 
when used in the standard spectroscopy calculation, gives a good 
description of the energy levels of these nuclei. Clement and 
Baranger (Git 68) using TaboZ^ln potential have obtained results 
for A = 18 nuclei very close to those of Kuo and Brown. In all 
the calcxilations described above pure harmonic oscillator wave 
functions and experimental single particle energies have invari~ 
ably been used. The use of harmonic oscillator wave functions . 
is preferred because in this representation separation between 
relative and centre-of-mass coordinates can be made exactly and 
two-body matrix elements of the nucleon-nucleon -potential can be 
easily calculated. Very recently, Kahana (Kaa-Preprint) has rep-* 
orted energy level calculations for A = 18 nuclei using Woods-Saxon 
Y'iSiVe f-unctions. Eo attempt till now has been made to obtain the 
single particle energies and single particle wave functions 
from some self-consistent calculations. HE calculation gives 
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us a picture very close to sliell model, i .e., it provides us with 

* 

a self-consistent single particle level spectrum and associated 
with each level is the single particle wave function. Ihese 
single particle parameters serve as a "basis for any shell model 
spectroscopy calculation and the residual interaction would "be 
some as that used in the self-consistent HF calculation to 
obtain the single particle parameters. We have carried out the 
calculations with effective Yale intera,ction and Sussex inter- 
action- The details of the calculations and the results are 
presented in Chapter V. 


In. a somewhat different approach followed by some 
authors one tries to dispense with any potential picture and the 
two-body matrix elements which enter any shell model spectroscopy 
calculations are treated as arbitrary parameters and adjusted to 
give best possible fit to spectra, of several nuclei. Such 
calculations ha,ve been done by Arima, (ira 68) for oxygen isotopes 
^^0, ^*^0. In his calculations ^^0 is regarded a,s inert 
core- and low-lying states of these nuclei are assumed to arise 




1 matrix 


from neutrons in (Td^y-g configuration, 

elements have been calculated from least square fit to these 
levels. * 1 ^ = 0 ma.trix elements have been caloulated by extend- 
ing the calculations to cover the nuclei F, ^F and Ue as 
well. Similar calculations have also been done by Federman. 
and lolmi (Fen 65). It would be interesting to compare the 
empirically determined two-body matrix elements with those 
calculated from a realistio potential for this configuration 
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us a pioture very close to siiell model, i.e., it provides us with 

* 

a self-consistent single particle level spectrum and associated 
with each level is the single particle w<ave function. These 
single ponticle parameters serve as a basis for any shell model 
spectroscopy calculation and the residual interaction would be 
same as that used in the self-consistent HF calculation to 
obtain the single particle parameters. ¥e have carried out the 
calculations with effective Yale interantion and Sussex inter- 
action. The details of the calculations and the results are 
presented in Chapter V. 


In a somewhat different approach followed by some 
authors one tries to dispense \vith any potential pioture and the 
two-body matrix elements which enter any shell model spectroscopy 
calculations are treated as arbitrary parameters and adjusted to 
give best possible fit to spectra, of several nuclei. Such 
calculations have been done by Arima (Ara 68) for oxygen isotopes 
ISq^ 19©^ ^*^0. In his calculations ^^0 is regarded a.s inert 
core and low-lying states of these nuclei are assumed to arise 


n-' _ 


1 matrix 


from neutrons in ^^1/2^ configuration 

elements have been calculated from least square fit to these 
levels. = 0 ma.trix elements have been calculated by extend- 
ing the calculations to cover the nuclei F, ^F and Ne as 
well. Similar calculations have also been done by Federman 


and Talmi (Fen 65). It would be interesting to compare the 
empirically determined two-body matrix elements with those 
calculated from a realistic potential for this configuration 
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space. Such a comparison, has heen made in Chapter 7. 

In Chapter II vre describe various types of nucleon- 
nucleon interactions available for use in any nuclear structure 
calculations. The prescription used to obtain the effective 
Yale interaction from Yale potential by Shakin et, al. and 
Sussex matrix elements directly from phase shift data by Elliott 
et. al. has also been described briefly. Tho details of the HE 
theory and various formulae used in the calculations are given in 
Chapter III. The results of HE calculation are presented in 
Chapter IV. The results obtained for the energy levels of 
various nuclei obtained with HE wave functions and different 
types of residual interactions are given in Chapter V. 

All calculations have been carried out on IBM 7044 
computer at Indian Institute of Technology, Kanpur. 



CHAPTEE II 


NUGLEOE-ITUOLEOE BTIERAGTIONS 

The exact nature of the nuoleon-nuoleon interaction is 
still not known and it is not possible to present it in a^ close 
deductive form. However, the nucleon-nucleon interaction in 
its nonrel ati Vi sti o version can be represented in terms of a 
potential. The nucleon-nucleon potentials, available for use 
in nuclear structure calculations, can be broadly classified 
into following categories: 

w 

(i) Realistic potentials 

(ii) Effective potentials 

(iii) Potential matrix elements obtained directly from 
phase shifts. 

We shall describe these potentials in the following sections. 

11*1 Realistic Potentials 

Realistic potentials are obtained directly from an 
analysis of the deuteron and two-body scattering data. The 
available data (Marr 68: Eos 67) consist of in pp and np scatter- 
ing: total cross-sections, angular distributions, and different 
types of polarization measuremaits and also s-wave scattering 
lengths. The two-body bound state data are deuteron binding 
energy, quadrup ole moment, D-state probability, etc. In the 
energy range 0 to 500 Me? there are about one thousand pieces 
of data left to be fitted, after discarding ijuite a few of them 
considered unrealistic. There are two ways of handling the 


pro'blem. One can either try to fit the data as they come from 
the experiments or one can use a phase shift parajaeterization of 
the data and compare these phase shifts with those extracted 
from potentials. The latter procedure is generally preferred. 
Phase shift analysis of the experimental data has been carried 
out mainly by tv/o groups: one at Yale led by Breit (Brt 62) and 
the other at Livermore (Marr 68). Signell and co-worhers 
( Sil 65) have also made a number of useful analyses. Only one 
set of phase shifts emerges and the phases are fairly well 
determined. It is to these phase shifts that the potential 
parameters are fitted in most of the models. Ho attempt is made 
to fit the data at energies above about 500 MeV. At hi^er 
energies pion production and other relativistic effects become 
important and two-nucleon Schrbdinger ec^uation is inadequate. 
However, the energies of nucleons in collision in nuclei are 
expected to be of the order of 160 MeY or so. Thus, we have 
reasonably accurate information about most of the partial waves 
which would play a dominant role in nuclear phenomena. 

The long-range part of the interaction is rather 
unambiguously Imown from meson- theory . It is determined by the 
one-pion exchange potential (Hun 57) corresponding to the light- 
est meson which can be exchanged between two nucleons. The pi on- 
nucleon coupling constant is measured from pion experiments or 
by nucleon-nucleon scattering data. However, at the pion range 
other terms corresponding to the exchange of heavier mesons, 
such as p 5 -me sons (Bm 67; Lon 68), dominate one-pion 



exchange potential (OPEP) hy a factor of two or three. Again, 
the one hoson exchange strength (OBEP) is determined hy the 
nucleon-nucleon scattering or alternatively "by calculating the 
meson-nucleon coupling constants, for example from SU^ symmetry. 

The nature of the short-range part of the interaction 
is not exactly known. The most general two-body interaction is 
a non-local one, which, to be translational invariant, must have 
the form 

. ( 2 . 1 ) 

The relative and centre-of-mass coordinates are defined by 

r = ^-r|, ^ = |(r^+r>) etc. (2.2) 

Meson-theory indicates that the correct two-body interaction 
may be non-local at short-range, but most of the potential 
models are local potentials, for which 

Y(r[r^) = Y^r) . (2.3) 

II. la' Local Potentials 

The first group of local potentials contair^ a repulsive 
core of infinite magnitude — the so-called hard core. The 
assumption that nucleon-nucleon potential becomes repulsive at 
short distances was first made by J as trow (Jaw 51 ). It was 
strongly supported by the t?/o-nucleon scattering data through 
the change in sign of the ^SQ-state phase shift from its low- 
energy positive sign (attractive potential) to a negative sign 
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oharaoteristic of a repulsive potential alove 200 MeV. Iwo suoh 
most widely known potentials are of Yale group of Breit and ool- 
laborators (Laa 62) and of Hamada and Jolinston (Haa 62). Yale 
potential is of the following form; 

(2.4) • 

where x is the intemuoleon distao.ee measured in units of pion 
Oompton wave length, p,i.e. 

X = p“^r , (2.5) 

( 2 ) 

yv >^(x) = one-pion ezoiiaJige potential^ 


•^p,( T- 




6(x)l 


( 2 . 6 ) 


where d' and are the spin end isospin vectors for the two 
nucleons, ^2 usual tensor-operator defined as. 




otherwise. 


( 2 . 8 ) 


Here (i=G, 1, t S, . q.) are the central, tensor, spih-orhit 
and quadratic spin-orhit potentials defined hy 



+ gO f or X / s - 

^ G * 

(2.9) 

is 0.35. She radius of the hard core is 

,3:^ = 0.5116F for singlet even and triplet odd, 

= 0,5002F for singlet odd and triplet even 

states. (2.10) . 

Ihe parameters of the potential can he found in referenoe 
(Mar 66). In some oases the idealized hard oore is replaced hy 
softer repulsive oore which is quite strong hut not infinitely 
large. The examples of soft core potentials are those of 
Bressel, Kerman and Lomon (Brl 65) and of Reid (Red 68). Outside 
the repulsive oore BKL potential has nearly the same form asihatof 
Hamada-Johnston (HJ ) with a few changes in parameters. 

The next group of local potentials is typified hy the 
boundary- condition model of leshhach, lomon aad Tuhis (feh 64). 

In this the one-pion and two-pion exchange potentials are used 
outside a radius r^ at which an energy dependent boundary condi- 
tion is speoified to fit the hi^-energy scattering data. 

inother group consists of me son- theoretic potentials 
which tahe into account exchange of several bosons 'in':, . , 

addition to -jt-meson,. such as p jio,Y^and o-mesons. Different 
combinations of masses and coupling constants give a variety 


Yj_(x) = 



for X >x^, 





of potentials. Recent attempts in developing such potentials 
Rave been done by Biyan aid Scott (Bm 67) and Green and col- 
laborators (Gm 67) i 

II, lb Ron-local Potentials 

A non-local potential is a non-diagonal operator, whioh 
implies that the' potential energy is a matrix in coordinate 
space. !Ehe Sohrb dinger equation then has the form, 

" 'li + ^¥(r^r’^) dr^ = E , (2.11) 

where m is the nucleon mass. Because of hermitioity, ¥('?■{'?■') 
must be a symmetrio matrix, real in the energy range O^^OO 
MeV, Such potentials are developed by Yamaguchi (Yai 54-), 
labakin (Ian 64), Mitra and co-workers and others. Yamaguohi 
potential is a simple separable potential of the form. 


= -?iU(r) UiT'). 


( 2 . 12 ) 


labakin potential is a series of separable terms of Yamaguohi 
type and is of the forms 


V(rfFO = 


f ^ ?.-f[-g(u)g(r0+h(r)h(r/)]^^(r)'l^^T3j(^0 


..■Ifej'TM 


m 


thl;42 . 


(2.15) 


Here, 


K 

m 


® (r) are the vector spherical harfflonios. 


-s jm 


defined in terms of tiie spiieriGal harmonios Y^(r) and spin 
states by 

r UlA 

(laii,^mg!JH) Y^*^(r)7(g®, (2«14) 

^6^'S ' ' 

( JM) is tiie GlebsiL-Gordbn ooefficient defined on page 37 
of reference (Eds 57). In expression (2,13) 'because the sum- 
mation is done OTer t and it allows t-mixing, i.e., especially 
S-D mixing. P.^ is the isospin projection operator. (Ehe g-term 
is attractive and h-term is repulsive. la’bakin has given the 
form of g and h in momentum space* iBiis potential gives only 
a qualitative fit to the two-'body data. Moreover, it does not 
have the one-pion exchange feature for the long-range part of 
the potential', Ihus, the labaMn potential is not very real- 
istic 'but it has been very useful in nuclear structure calcul- 
■••ations because of the ease with which its matrix elements can 
be calculated, 

Il^lc Velocity-Dependent Potentials 

Such potentials are of the form 

V = V(r, p) + constraints for hermiticity, (2.15) 

where p is’ the momentum operator. Ihe most commonly, used 
velooity-dependent potential is that due to Green (Gm 62) and 
is of the form 




W(r) and U(r) iiave tlie Gaussian radial dependence. 

II. 2. Effective Potentials 

It. is well lmo\m that the approxination in which one 
considers nuclear parti clos moving in some average potential 
field and .interacting with a weah residual interaction. . . is • 
■highly successful in the study of nuclear properties.. ..In an 
independent particle model oaloulation one u sually... considers ■ 
the effects of only few shells, Ihe neglected shells renorm- 
alize the interaction in the space of. the shells considered in 
one ’ s cal cul ati on , Reno rm ali z ati on o f the nuol e on-nucl e on 
interaction also takes. place due to admixture of veiy hi^ 

...energy, orhitals into the nuclear wave functions, Ihis results 
from strongly- .repulsive nature of the nuclear force at short 
distances . and also from specific characteristic of tenser 
force*- In fact^ the singular nature of the nuol eon-nucleon 
.interaction-, at short distances gives rise to strong short-range 
correlations which are not properly described by an independent 
particle . wave function, Thus, in the nuclear structure calcul- 
ations one must . use the unmodified or the free nucleon-nucleon 
interaction in, a '.configuration space described by wave functions 
which properly take into account short-range correlation effects, 
iltematively, one can use a modified or effective interaction 
■ in ...a- truncated, configuration s.pace oonhL sting of ..independent 
particle wave functions... Ihus^ the problsa is to find ah 
effective interaction which would give the same results 

when used with uncorrelated wave functions in a limited 
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configuration space, as the true interaction Y would give with 
correlated wave functions in an unlimited configuration space, 
For many years, in some simple minded calculations was 

chosen in an empirical manner and only a vague idea existed 
about its possible dependence on the choice of sin^e particle 
orbitals and configuration space and its relationship to the 
realistic two-body interaction. It was only after the develop- 
ment of Brueckner theory (Brr 5 5) that several different appro- 
aches to the relationship of Y^^^ and Y were given. In the 
Brueckner theory the nucleon-nucleon interaction Y is replaced 
by an effective interaction, commonly known as K-matrix and is 
defined in equation (1.8). Several different approaches to 
relationship of Y^^^ and Y, ba.sed on Brueckner theory, have 
been given by various authors. lo name a few are the separation 
method of Moszkowski and Scott (Moi 60), reference- spectrum 
method of Bethe, Brandow and Petsohek (Bee 65) s Eden-Emery 
method (Edn 59), etc. 


The reaction-mo-trix differs from the two-body potent- 
ial over the short-range of correlation distance in the nuclear 
wave functions. Ihere is an alternative way of introducing 
the effect of short-range correlations in the effective inter- 
action, namely the method of canonical transformation. Ihis 
approach (which in essence is close to Moszkowski-Scott 
technique) has been followed by Shakin, Waghmare and GO-workers 
(Shn 67a). We shall briefly describe below their method of 
obtaixiing the effective interaction. In their method* the 



correlations in the wave function are introduced via a unitary- 
mo del- ope rat or such that if the long-range properties of the wave 
function, are specified by a function "r^ . , . the 

corresponding correlated state is , 


(p(^, ^ ^n) "" "^n^ * (2.17) 

An effective Hamiltonian is defined in the space of the uncor- 
related functions as 


^eff " 


e 


iS ggiS 


( 2 . 18 ) 


where H is the true nuclear Hamilto2aian.. Since a transformation 
of the form described in equation ( 2 . 18 ) leaves the eigenvalue 
spectrum uncharged, H .and arc completely equivalent with 

respect to energy calculations. But,,besides liLe desired two- 
body terms, contains nany-body terms induced by the unitary 
transformation. Carrying out the expansion of and neglect- 

ing three and many-body terns, one finds that the effective 
Hamiltonian is. 


H 


eff 






n^n^n^n^ 


l“3V V 


4 3 

(2.19) 


where H^^^ has been expressed in the second-quantized notation. 
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The operators a'*’ and a are the nsual creation and annihil- 
ation operators for the single particle states. The sub- 
scripts n^ and 112 refer to the quantities (n, h , jj m) 
necessary to specify the orbitals for particle motion in 
hamonio oscillator potential. T is the kinetic energy 
operator and is the nucleon-nucleon interaction which 
may or may not ha.ve a hard core. The potentials and Ug 
are used to generate the single particle wave functions. 

The form of U is goneral-ly taken to be of harmonic oscillator 
type to simplify the calculations, i.e., U = ■^ kr^. The 
operator S is assumed to introduce only central short*- range 
two-body correlations in miich case equation (2.17) becomes 


^^ 1^2 


H^n n 


( 2 . 20 ) 


whereas are the unco rr elated wave functions and(|/s are the 
corrected wave functions occurring in equation (2.19). Y ^2 

can be written as a sum of two po.rts, v^g which is diagonal 
in relative orbital angular momentum I and v^^ which is off- 
diagonal in l>, i.e. 

h2 = ^2 + (2-21) 

Clearly v^^ receives contributions only from tensor part 
of the nucleon-nucleon interaction, v ^2 is separated into 

8 B 

a long-range part v^'g soil n. siiort-rmge part v^g such that 
the short-range part produces no energy shift in the pair 
state. This is a natural extension to finite nuclei of the 
zero phase shift condition of Moszkowski and Scott (Moi 60) 



for nuclear matter calculations. The main effect of 
is to admix high-'Ciomentum orbitals into nuclear wave function. 
This admixture is not readily treated by the separation method 
used for v^ 2 ‘ Thus, the tensor correlations are treated on 
a different footing than the correlations due to short-range 
part of v ^2 fhe separation of v^,^ into v^^^ v^-'g 

achieved by solving the following two equationss 


(T^+T2+U^+U2+v®2)q}n^n 

_L- 1 

and (li+T2+"Ji+h2-)5 


= ( S ,- n 2 ^$ n , n 2 

= -^_^ ni + S 2 ^$ nin 2 


( 2 . 22 ) 

(2.23) 


and s_ are the energies of single particle states Ut 

n^ n2 vj - j_ 

and n 2 respectively, for hard core potentials such as the 
Yale or the HJ potential it is necessary to introduce a short 
range pseudopotential (VP) to carry out this separation 
procedure. In that case 


Vi2 = (vj2 + ^^^12 " » 


(2.24) 


and one requires that 


(Tj^+l2+U3_+ll2+v=2+VP)^ 


(e„ +e„ )tl'„ ^ • (2-25) 


“ l “2 '"n 


In general, therefore, one has 


^eff t^2)> ^2 


^4* _4“ 


n^n2 


a„ a. 


nnLn2^3^4 


1 ^2 ^4 ^3 




( 2 . 26 ) 



If the shell-model or Hartree-Fock calculations are limited 
to orbitals having only fairly low-momentum components, almost 
all of the contribution of v^^ to the binding energy and to 
the effective interaction would be missed. Most of the 
contribution of v^^ may be talcen into account by renormalizing 
the effective Hamiltonian in the independent pair approxii^ 
mation (Kuo 65): 


H 


eff 


^ I ^ 1^2/ ®n. ^ ^ 


V2 






<? 




- TO) + vf+vf I _ > (2.27) 

5 4 


If the second-order term in the p sen. dope tent i al is also 
included, a term VP VP should be added to in equation 


(2.27). 


The prescription described above is quite general 
and is applicable to any singular interaction. Using 
expression (2.27), Shakin et. al. (Shn 6 7a, 67b) have evaluated 
the matrix elements of the Yale potential for various states 
of relative angular momentum. An attractive square-well 
type of pseudopotential (VP) has been used for those states 
of relative motion where the interaction is repulsive and 
also in Sq and states (to keep the separation (hLstance 
constant at 1. OF) . For the treatment of intermediate states 
in the calculation of second-order terms the approximation 
of Kuo and Brown (Kuo 65) is used, which consists of using 


plaJie-waTe intermediate states and an angle averaged operator 
to taJse account of the Pauli principle in the intermediate 
states. A simple parameterization of the energy denominator 
e is made, fhe details of these calculations can he found in 
references (Shn 67a, 67h). Results for the relative matrix 
elements of effective Yale interaction are tabulated (Shn 67h) 
and will he used in the present work. 


It can he, seen that the unitary-model-operator 
method of obtaining the effective interantion is in some 
sense analogous to Bruecknor theory. In fact it is possible 
to find a unitary transformation (Grg 68) for -vdaich the 
transformed two-body potential is 


Vgff = f(K + K+) 


(2,2S) 


Expansion of defined in equation (2.18) is 


H 


eff 


= + .... , ( 2 . 


29) 


where is an n-body operator. In first-order 


^eff = 


C?.30) 


when many-body terms are neglected. Applying the Hartree- 
Eock formalism to in equation (2.30), one finds the 

same expression for the ground state energies as in the 
first-order of Brueckner theory, except for the replaGement' 
of K by the hermitian combination ^(E+K’*’) . The higher- order 
corrections in both oases arise due to many particle effects, 



more accurately, the many-hody terms in the result of a 
unitary transformation have their counterparts in the cluster 
corrections of Brueckner theory. 

There are quite a few other effective interactions 
oh tained directly from realistic interactions, for example, 
the effective interaction of Kuo and Brown (Kuo 66) derived 
from HJ potential and later improved by Kuo (Kuo 67)* Their 
effective interaction includes the contribution of renorm-j^.-,. 
alization effects due to core-polarization and has been 
widely used in nuclear-structure studies. To name a few 
other effective interactions used in the literature, are 
those of Wong (Y/og 68) obtained from HJ potential and of 
Clement and Baranger (Olt 68) derived from Tabakin potential* 
These effective interactions have been obtained following 
the lines of Brueckner theory. 

Of late, there has been a tendency to develop 
potentials primarily for the use in nuclear structure calcul" 
ations. The principal aim of these works is to determine- 
an effective potential having no hard core singularity and 
giving the correct energy, density and perhaps the symmetry 
energy in nuclear matter. Only a moderate fit to the two- 
body data is required, in additional requirement is that 
second and higher-order contributions to the energy of 
nuclear matter be small in order to improve convergence of 
perturbation expansion. These potentials are effective 
potentials in some sense and it 'is presumed that they are 



not very diiferent from the realistic potentials so that two- 
body data still have some significance even for the effective 
ones. We mention in this category the v/orh done by Davies , 
Krieger and Ba^raiiger; hestor and co-worlcers and Donelly, etc. 

Ihe potential of Davies et. al. (Das 66) is a simple 
velocity-dependent potential of the form given by Green (see 
equation 2.16) . Iheir potential acts only in relative S- 
state. The parameters are determined from a fit to nuclear 
matter data but such a potential is in strong disagreement 
with the two-body data. In an attempt to improve the agree- 
ment with the two-body data Krieger and co-workers (Krr 66) 
have introduced a. density-dependent term in Davies potential. 
The strength and range parameters of their potential are 
determined from si fit to the two-body data? while the density- 
dependent term tal-ies care of nuclear matter saturation prop- 
erties. The potential of Destor et. al. (Her 68) is also of 
the form given b 3 r Green but contains additional spin-orbit 
and tensor terms. These potentials have been used in the 
calculations of finite nuclei (Das 66: Krr 66; Her 68) in 
the Hartree-Hock approximation and found to give reasonable 
results. 

Brink and Boeker (Brk 6?) have investigated various 
effective interactions of different forms, e.g. a sum of two 
Gaussians, a sur:i of a Gaussian and a delta function, and a 
velocity-dependent interaction. The parameters of each potent 
ial are adjusted so as to give the e3rperimental binding ^ ^ ^ 



energies of ^He pjid nnclenr matter at about the experimental 
densities. With the resulting forces the binding energies 
and densities of intermediate nuclei are calculated. But 
no force could reproduce all the experimental data at the 
same time. 

II. 3 Potential Matrix Elements Obtained Directly 
from Phase Shifts 

In most of the nuclear structure calculations the 
primary input needs matrix elements of nuclear force between 
states of two nucleons in nuclei, Ihe usual procedure 
described above is to construct a potential model from 
nucleon-nucleon scattering data and then use this potential 
to calculate matrix elements needed for nuclear structure 
studies. In an alternative approach relative matrix elem€nts 
of the two-body potential are obtained directly from experim- 
ental phase shifts data without attempting to establish a 
realistic potential form as an intermediate step.: This method 
has been used by Kallio (Kao 65) , Koltun (Kon 67)'» Elliott 
and co-workers (Elt 67, 68a). Elliott et, al. have obtained 
the relative matrix elements of the two-body potential in an 
oscillator basis by using an auxiliary potential method. 

Their matrix elements are commonly knovm. as Sussex matrix 
elements . We shall briefly describe below their prescription 
used to obtain these matrix elements. 

The potential IT is regarded as a most general non- 
singular potential consistent with the general conservation 
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laws of spin, parity and isospin. It may le written as 
ir = yStL'j (r, ^)|sl'j>< sEjI , (2.31) 

S t^J 

where S6j denotes the particular spin and angular momentum 
state of the pair and n the oscillator radial quantum number 
describing the relative radial wave function E^p(r) of the 
pair. The .Slj labels are referred to as a particular 
channel. Thus, the set of the matrix elements 
<n^ Sl!jfVlnS?.J> axe 

<n* S rJlVlnS.{iJ> = J (r, . 

( 2 . 52 ) 

In a particular channel Y (the label StJ can be dropped for 

brevity if one is confined to a chosen channel) ^ written 

as the sum of tv/o terms V = where is the auxiliary 

potential and V-, is small in the sense that it can be treated 

■ >, 2 , 2 

in Born approximation at some energy E = — , where k is the 

relative momentum and m is the nucleon mass. If 6 is the 
observed phase shift and 6^ that calculated with the auxil- 
iary potential at the same energy E, then it is shown by 
Morse and Eeshbaoh (Moe 55) that 

^(^Uj^(r) Vi(r, u^(r) r dr = - tg(6-6^) , 


(2.35) 



where Uj^(r) is trie radial v/ave function for scattering by the 
auxiliary potential with the usual asymptotic normalization 

u^(r)— s_L/i (hr - + 6^) • (2-. 34) 

Ihe auxiliary potential is taken to he a cut-off oscillator 
2 

Vq = - f- (D - rV4b^) in r^a 

= 0 in r >a (2,35) 

with depth D, r.ange a and shape or range parameter b. The 
radial wave function n^(r) is then a solution of the following 
Schro dinger eq_uation 


( 




CL 


dr 


.t( mi) 

2 

r 


+ D + k' 


^2 

■^-)uUr) = 0 (2.36) 
4b^ ^ 


for r4 a and is the free nucleon scattering wave function 
for r_^a. For energy 


-1 ,<2 ^ 2 ,. ^ 2 , 2 

^ = {2n+l- , (2.37) 

mb m m 

the radial part U|-(r) is given by 

Uj^(r) 7= B ®]^^(^) LT ^ a 

= k(cos c^ 0 |^(kr) - sin 6^%(kr) r ^a , (2.38) 

3 ^(kr) and ■r)|^(kr) are the spherical Bessel and leiuaann functions 
defined on page 621 of reference (Koe 53), is the radial part 
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of the hound state v/ave function of an infinite oscillator 
potential.;Uj^(r) has the required asymptotic normalization 
given in equation (2.34). Constants B and 6^ are determined 
hy the matching conditions 

® = li:(cos 6^ 3^(ka) - sin 6^ n^(ka)) , (2.39) 

B^^(a) (cos 6^ j^(ka) - sin 6^ T]^(lca)) 

a) ~ Tcos '6'^ 'f^”(ka) - sin 6^ (ka) ) * (2.40) 

where the prime denotes differentiation with respect to r. 
Combining equations (2.33) and (2.38) v/e have 


( h = - iij tg(6-6g 

mi) 


ii^k 


0 




^ ^ u^(r)j.r^dr. 

(2.41) 


Ihe second term on the right-hand- side .of equation (2.41) is 
referred to as long-range correction. It arises due to the 
fact that i^jj(2^) is not an o.scillator function beyond r >a. 
However, this is very small if a is chosen to be near the 
range of nuclear forces and can be neglected. Ihe matrix 
elements of the total potential ¥ are obtained by adding the 
matrix elements of ¥^ in the range O^r^a in equation (2.41) . 
In the range a:^r<(pc>, V is same as ¥^. Ihus, 


k 
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Y'vfJ 

Q.) H^^(r) r^dr 


o 

mB 


■fcg(&-OQ) 


or more compactly 

\t> = • (2-«) 

‘-‘ ‘■’ mB 

The matrix elements of Y are insensitive to the choice of 
auxiliary potential as long as 6 is close to 6^. Equation 
(2,42) gives matrix elements diagonal in radial quantum 
number. The off-diagonal matrix elements are obtained by 
differentiating the diagonal matrix elements v/ith respect to 
b. The theory described above has been extended to cover the 
mixing of channels so that the matrix elements in the coupled 
channel can also be calculated. The details can be found in 
reference (Elt 63a) . 

Thus, the auxiliary potential method described above 
consists in assuming that the difference between a rou^ly 
chosen, smooth auxiliary potential, and the real potential may 
be treated in Born approximation. In the earlier work of 
Elliott et. al. (Elt 67) and in Kallio method it is assumed 
that the entire potential V is small enou^ to be treated in 
Born approximation. Thus, the auxiliary potential = 0 
and the radial wave function Uj^(r) will be spherical Bessel 



fiinc^ions for all 


Uj^_(r) = (kr) (2.43) 

Lallio uses an approximate relation 'betY/een spherical Bessel 
functions and harmonic oscillator radial functions given^by 

= B'k 3^^(kr) , ‘ (2.44) 

2 2 

E = = ( 2n+t^-^^^xJ . (2.45) 

The above approxim.a.tion is valid at short distances and 
improves with increasing n. Thus^ Kallio matrix elements 
are given by 

<Sni''fl\L>= ^ ^=8* • (2-«) 

Kallio technique breaks down at 6 = %/2 since the matrix 
elements become infinite. In fact Kallio method can be 
regarded as the limiting case of the auxiliary potential 
method described above. 

The Sussex matrix elements are obtained using 
experimental phase shifts of Breit et. al. (Brt 62) and are 
tabulated in reference (Elt 68a) , The range of the oscillator 
length parameter is 1.4 ^■2.6K, At each value of b, the 
complete matrix, i.e. for all (n, n' ) would define the 
potential uniquely. A connection between the n value of a 
diagonal matrix element and the ener^ of the phase shift, 
which governs its value, is defined in equation (2.37). 
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Oonsequently , an upper limit of 300 MeY, up to which experim- 
enttil phase shifts are reliably Imown, imposes an upper limit 
on 11 . Moreover, since the .accuracy of the successive dif- 
ferentiation on a curve determined numerically over a finite 
range goes on decreasing, there is a maximum limit of the 
size of the difference n^-n, for which roliahle matrix elem- 
ents can he deduced. However, the matrix elements which can 
he obtained reliably are the ones which are of importance in 
most of the nuclear structure calculations. 

In obtaining the Sussex matrix elements it was 
assumed that • Y is non-singular. However, the matrix 
elements calculated by the auxiliary potential method give 
a first-order approximation to the reaction matrix elements 
even when the true potential is singular. This has been 
shown explicitly by Mavromatia et.'al, (Has 69). These 
authors start with a reaction matrix (t-matrix) associated 
with a two-body interaction Y and defined as 

t = V + Y § t . (2.47) 

The energy denominator e is , . 

e = Ho - E , , (2.48) 

Sq is the unperturbed Hamiltonian of relative motion 

iIq = T + U (2.49) 

with eigenstate 9 and 1 is energy of the perturbed 



7/111011 cp approximates. Q is the Pauli operator as defined 
earlier. If the normalization is chosen such that 




(2.50) 

then tcp 

= Yvh 

(2.51) 

and <( 9 |t j 9 ^ 


(2.52) 


Ihe matrix elements <^ 9 |-t J 9 ^ remain finite even when Y has 
a hard core. Of interest for nuclear structure calculations 
are those t-matrix elements for which is the harmonic 
oscillator Hamiltonian, 

i m(^r^ (2.53) 

and the basis function cp is the complete set of harmonic 
oscillator wave functions. Such a t-operator is defined as 
t^. In the auxiliary potential method of obtaining reaction 
matrix elements, one divides the .potential V into two parts, 

A 

the low-energy part V for which the exclusion principle is 

• A 

significant, and the hi^-energy part (Y-Y) . for hi^-energy 
part Q is approximated by Q^, the appropriate operator for the 
two particle system, This is justified since intermediate 
states of high-energy are largely -unfilled in the nucleus and 
so exact treatment of Q is not very important. One define 
the reaction matrix associated with the hi ^-energy part as 
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with. e. = !P 4 . Y' - E 


(2.55) 


Keeping terms to second-order in t^ and Y the perturbative 


expansion (Broii 54) of t^ is 


^n - 


h - T . ^ n K i- 1. . t, ^ n t.(2- - % t, 

A 


n 


^A ^A e:.' ^ ^ e,/ ^ 


n 




with e^ = i ^ - E 


( 2 . 56 ) 


(2.57) 


In Elliott method, V = ¥q is taken to be a cut-off harmonic 
oscillator defined in ec[uation (2,55). The basis functions 
9 associated vd.th t^ have the radial part Uj^(r) defined in 
equation (2,38). If the ' radial part of the wave function^ 
is denoted 'as Wj,.(r), the matrix elements 
<Uj^( r) 1 I uj.( r)^ = <u^( r) ) ( V-V^) j x)y 

may be obto.inod from the solution of follov/ing Schrodinger 
equations for U 2 ^(r) and WT^(r): 


[„ -i-- + ^ . 4 . Lkl I + Vq - B] (ru2j.(r)) 

m dr"^ m r*^ 


0 


[- 1^ ^ + V - ElCrwPr)) 

“ dr'^ ¥a 


0 


( 2 . 58 ) 
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One obtains 


= J dr ruj^(r) (Y-Yq) rwj.(r) 


K 

m. 




K 

m 


k tg{6-6^) 


(2.59) 


Thus, using equations (2.38), (2.56), (2,59)and neglecting ■ 
any long-range corrections we have, to first-order in t« and 

ii 


Y. 




*2 k 


m 


^ tg (5-S„) + 


+ higher-order terms. 


(2.60) 


With a, suitable choice of Y^, it is generally possible 
to make 6^=6 for any particular value of k in which case the 
first term in equation (2.60) vanishes J is 

given in first order by even v&ien the true 

potential is singular. This is precisely the prescription 
used in the auxiliary potential method TAiiere, however, it is 
assumed that Y is well behaved and the tv/o first terms on the 
right-hand-side of equation (2,60) together then give the 
potential matrix elements Y fen « Thus, in any case 

the Sussex matrix elements can be regarded as a first-order 
approximation to nuclear reaction matrix elements , 
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There is a direct analogy between the Moszhowski- 

Scott separation method and the auxiliary potential method 

of calculating reaction matrix elements. In the separation 

method Y is divided into a short-range part Y® and the long- 

range part Y^ such that the matrix elements of the reaction 

matrix t associated v/ith Y^ vanish and in first-order t is 
® n 

replaced by Y^. In the auxiliary potential method Y® can be 

_ A ? A 

associated with Y-Y and Y^ with Y, Matrix elements of t, 

A 

associated with (Y-'f) can be made equal to zero in which 

I . 

case t^ IS replaced by Y, However, the separation in the 
auxiliary potential method is not spatial but is energy- 
separation. Unlike the separation method the auxiliary 
potential method can be applied equally well whether the 
potential Y is attractive or repulsive, 

A 

In Zallio method Y = 0 so that 

"^^nthnl^nt^ ^ ^ • (2.61) 

This is also the expression to which right-hand- side of 
equation (2.60) reduces in the limit 6^’^0. In fact Sussex 
matrix elements with 6^^ 0 are substantially in agreement 
with Kallio matrix elements. Koltun phase shift method 
is a calculation of t-matrix for the free two. nucleon system 
in an oscillator basis, but we shall not go into the details 
of his method. ‘The work done by Elliott and co-workers is 
most intensive and complete. The matrix elements have been 
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tested in soae of tlie nuclear structure calculations for 
finite nuclei and found to give reasonable results. In the 
present work these matrix elements have been employed in 
the HI' calculations of light nuclei in the range 4, 4 . A ^40. 

II. 4 S-ummary 

lo summariae, the situation is as follows, 2he 
experimental information about the nucleon-nucleon potential 
comes from the two-body scattering and two-nucleon bound 
state data. Most of the potential models nowadays use the 
results of meson theory, mainly to determine the long-range 
part of the potential. Ihe short-range part is left to 
phenomenology and the uncertainty in it is mainly due to the 
necessity of worhing below an upper energy limit of 300 MeV, 
The adjustable parameters of the model are determined from 
a fit to two-body ■ data. One important characteristic of 
these potentials which fit the two-body data, is that they ^ 
do not fit the nuclear matter data, e.g. in the binding 
energy calculation of nuclear matter Yale and Hamada-Johnston 
potentials give (Brr 62; Eay 63 ) a minimum of -8.3 Me'f per 

particle at a density corresponding to fermi momentum ■ 

—1 "d. 

ISj, = 1.191 and -7.8 MeY per particle at kj, = 1.121 respect- 
ively, compared to the experimental values of B.E. per 
particle = - 15.75 MeY, kj, = 1,51“^. According to Brueckner 
theory, following features of the nucleon-nucleon potential 
are responsible for the smaller values of energy and density; 



(i) large core-radius 

(ii) strong odd-state repulsion. 

(iii) quadratic-spin-orbit tern 

( iv) weak even triplet central force. 

However, tiiege a,re just the same features necessary for a 
good fit to hi^-energy scattering data, labakin potential 
also gives only -*8 MeV per particle, k^, = 1.6P“^ in nuclear 
matter calculations, in first orderi Ihese potentials when 
used in any perturbation theo:^ treatment of finite nuclei 
fail to give adequate binding in first-order (Ken 66; Shn 67b) 
It has been found that second-order corrections for these 
potentials are quite important and reasonable values of 
binding energies of finite nuclei and infinite nuclear 
matter can be obtained only when the contributions of second- 
order terms are taken into account (Ian 64; Pal 66), 

On the other hand there are some effective potent- 
ials, available for use in .binding energy calculations, which 
are computationally easy to handle. She most important 
criteria in choosing these potentials are that they fit the 
nuclear matter data and the second-order corrections to the 
binding energy of nuclear matter are small. Binding energy 
results obtained with these potentials are, in some cases, 
better than those obtained with realistic potentials 
(Das 66: Krr 66) at least in first-order, but then the fit 
to the two-body data is poor. Ihus, the question still 
remains open whether a potential model can be constructed 
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wMcli would fit tlie scattering data as well as the nuclear 
matter and finite nucleus data at one and the same time. 

The fact, that different realistic interactions 
which fit the tvro-hody data q.uantitatively give similar 
results in nuclear structure calculations (Lyh 67) , suggests 
that one can avoid the complicated procedure of constructing 
a potential model and obtaining the parameters from a fit to 
two-body data. Instead, one can obtain the matrix elements 
of nucleon-nucleon potential in some suitable basis directly 
from experimental phase shift data as done by Elliott et. al. 
and others. However, the results of such calculations would 
be expected to be close to those made with realistic potent- 
ials (see Chapter lY) . 



GHAPTEE III 


HAEIESE-EOCK IHEOEY EOR EIIIPS lJUGLEI 

III.l Hartree and Hartree-Eock Equations 

Present microscopic theories of nuclei are based on 
two postulates; 

(i) She degrees of freedom associated with the meson 

fields in nuclei may he replaced hy the potentials 
(including exchange forces) between nucleons. 

(ii) Tv/o-body forces are most important. Ihus, the 

Hamiltonian of the nucleus consisting of A nucleons 
is of the form 


A 


H = > 


Ii +> 

i=l ±'<k 




(3.1) 


■^i th 

where is the kinetic energy of the i nucleon, 


A 


ik 


is the interaction potential between the pair of particles 
i and k, and is the nucleon mass, To understand the 
nuclear structure one has to solve the many-nucleon Schrodinger 
equation 


H^(^, ^2 ... rj, ..-.r^) , (3.2), 

where r^, define the coordinates of the particles, 

Phis is a non-separable differential equation in 3A independ- 
ent variables and needs some simplification for its practical 


treatment. Hartree and Hartree-Fock (HF) approximations 
(sir 60) consist of obtaining approximate solutions of 
equation (3.2) using variational principle. In Hartree 
method, the approximate groimd state solution 'pg of the 
above equation is assumed to be given by the simple product 
of A single particle wave functions (orbitals) 

first A states being regarded as 

occupied so tha.t. 


Vo ...Vp = • 

(3.3) 

The set of functions -® varied to minimize the average 

energy (pp, system. The minimum cond- 
ition E'pg. /* = 0 leads to the following Hartree 

equation for the single particle wave function 




T 


T 


H 


(3-4)..' 


where the summation y ^ns over only the occupied single 
particle states. It can be seen that Hartree product function 
does not satisfy the Pauli principle according to which 
the wave f-unction must be antisymmetric under the exchange 
of any two nucleons. In the HP approximation many -body 
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wave fimction^ is taken to be fully anti symmetrized product 
of an ortbonormal set of single particle wave functions , 


y A! 


.H^^(r^) S>]_(i^) 

■ « « « « 

9 « 

• . • 

♦ mr 



<• 

rfl 

• 

m 

« • 

• 0 

* 

* 

^a(^) 


(5.5) 


are varied to minimize the energy E =<^(p, ^ of the 

system maintaining the subsidiary conditions, namely, the 
normalization of all ^^s and the orthogonality of any two 






(3.6) 


Thus, the HP eq.u.;-tions are gotten by demanding 

5[<*p, = 0 , (5.7) 

The details of the variational calculation are not presented 
here and can be found in any book on many-body theory 
(Len 64; Bron 67). The result is that single particle wave 
functions Uj’s satisfy the following HP equations. 

I (X ■ < 

k 

<a|Tlp;> + ^ <“Tbjii,|PY> = ’ 

■ Y=1 


where the anti symmetrized matrix element of V is 



and 


= <aYj7[p6^ - ^aYjYjdp^ 

<“l^lP> =|s^a^^I) % > 

<aYlTp> = V(fr. 

^ 5 (^ 2 ^} d?3_<ir2dr^dr2'' . (3.8) 

Tile HI' self-consistent single particle potential is defined as 

. A 

<“iuiP2 = (5-3^ 

Y=1 

Tlius, 

<<=^|H|P>_ <“|T|P2 + <“ln|P> = ' 

( 3 . 10 ) 


where h is the EE' single particle Hamiltonian, It can be 
shown that the* miocctipied HF single particle orbitals, which 
correspond to the states ayk, A, also satisfy equations 
similar to (3.10). That is, 

<^a{hjp> = for a>A, p>A , 

but <^ajh[p^ =0, if only one of the two states a and p are 
occupied and the other one is unoccupied. Thus-, the HF 
Hamiltonian h has no matrix elements between occupied and 
unoccupied levels. The ground state energy of the system 
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is then given hy 

V=1 Y=1 

A' 

= >__ [ <y|i;1y> +1 <y|u|y>] • (3.31) 

Y=1 

Equation (3.10) looks like a one-body Schrbdinger equation 
but the potentia.1 U depends upon the very single particle 
wave functions one has set out to find, lo carry the sum- 
mation over Y equation (3.9) we need to know the wave 
functions of all the states which are occupied in^ . Hence f 
the solution of HE single particle equovtion is obtained by 
the method of self-consistency (iterations) over the wave 
functions. These equations can be solved by choosing any 
one of the follo’ving two representations for wave functions 
^^s, described in next section. 

III. 2 Different Ways of Doing HE 

HE equations can be solved in two different repres- 
entations. 

III. 2a Coordinate-Space Representation 

In the coordinate-space representation the single 
particle wave functions are just written as functions of i?: 
equations (3.8) and (3.9) lead to the following Schrbdinger 



equation when Y is local 


P: 


A 


A 








Y=1 -•' 


( 3 - 12 ) 


The first terti is a local one-body potential (the direct, or 
Hartree potential) , whereas the second yields the exchange 
term, a non-local one-body potential (exchange potential). 

If the two-body interaction is already non-local, all one- 
body potentials will also be non-local. 

Oalculations in this representa,tion have been done 
by Kerman and Loclcett (Ken 64), and Vaiitherin and Vene^roni 
(Van 67). However, they have used a purely central two-body 
potential. Inclusion of a spin-orbit term or a non-locality 
in tvro-body potential makes the calculations enormously 
complicated, Ihis representation has also been used in all 
calculations of Brueckner et, al. (Srr 61) and Kohler (Kor 65) 


III.2b Oscillator Representation 

Ihe HP single particle wave functions are expanded 
in a convenient set of orthonormal functions 



Y/liere a and \x denote the angular raonentun and isospin quantum 
nuiahers of single particle states, i,e. a = (aiomx) , 
(p=nE 3 mT): a, ii are radial quantum numbers. Combining equ- 
ations (3.13) and (3.10), we get 


1 ■ 




where ^ ^ <^‘l^ 2 hAHl^^y 2 > 

Pg \^2 ' 


U2^2 = J ■ 


( 3 . 16 ) 


^ is the densit^r of single particle levels and the summation 
runs over the occupied states only. The total energy Eq is 


then given by. 


(3.17) 

The expansion coefficients g“’s can be found by diagonalizing 
a matrix whose elements are <(p^jhtp'^ . Thus, this is also 
known as the matrix method (Bar 63: Net 63) of solving HP , 
equations. The choice of the basis cp in equation ( 3 . 13 ) 



is completely ax-bitrary "but the choice of thxee— dimensional 
haxmonic oscillator wave functions is preferred for two 
reasonss 


(i) The H? single particle functions are expected to be 
close to oscillator ones, so that a few terms of the series 
may be needed to get a good representation. In other words, 
the set of oscillator ftinctions used need not have a large 
dimensionality and instead of solving integro-differential 

e(iuations of the type (5.12) one has to diagonalize small 

0 

real symmetric matrices, Thus, the calculations are much 
ea.o.ier. The coefficients and energies e’s are always 
real since they are the eigenvalues and eigenvectors of a 
real symmetric matrix. 

(ii) In this representation the sepa,ration between the 
relative and centre-of-mass coordinates can be done exactly. 
Thus, the calculation of two-body matrix elements of nucleon- 
nucleon potential is easy since the latter generally depends 
only on the relative coordinates of two nucleons. 


This representation has been used in most of the 
earlier HI calculations described in Chapter I as well as in 
the present work. The function is the normalized solution 
of the three-dimensional isotropic harmonic osoillator: 






f 


( 3 . 18 ) 



m 

is the eigen function of total angular momentum defined 
in (2.14) , A- is an isospin function and the radial function 

_X ' 

is 




2fpJ(n 


) 

^ i 

h 


J^(n+^+^-) ^ 


(Yx) 


5 .+I g-X/2 


a+^) 

' 11-1 


(X) 


(3.19) 


2 

X = Po^ > n = li, 2^ 3 X^(X) is a Laguerre Polynomial, 

defined on page 784 of reference (Moe 53) and 
being the oscillator frequency, ' The energy eigenvalue 
associated with this state is e^^.. ~ ( 2n+i-4-)hui. 


In general, the summation over p in equation (3.13) 
is over all the quant'um numbers, Hov/ever, the HP Hamiltonian 
h can have various symmetries, such as spherical, ellipsoidal, 
axialf’etc. which when considered in the initial step, can 
reduce the vast variational space in expansion (3.13). If 
we consider the IIF field . to be sphericallj’’ symmetric, then 
the single particle HP states have good angular momentum and 
the c“ (and hence A, ,,/) are diagonal in t, j and independent 

ji \ 

of m ( z-pro 3ection of Thus, the expansion (3.13) can be 

rewritten as; 

. N 

\poc(r) G|[( ' (3.20) 


where 1 specifies the maximiim number of terms which are 



JZi 


included in ilie SLuimation in equation (3.13)» and is same as 
tlie maximum number of relative nodes in functions occurring 
in tile summation, liie mixing of neutron and proton states in 
a given orbit iias not been considered so that expansion of 
equation (3.20) is limited to a single value of third comp- 
onent of iso spin x which distinguishes the neutron and proton 
orbits. Also, if there is no coulomb force G^s are indep- 
endent of t , 

I I I . 3 Oalou l atio n of the Matrix Elements 
in 'Osciira)tb~r Basis 

III. 3a One-Body Matrix Elements, "v ^ • 

The matrix elements of the kinetic energy operator 
T in the oscillator basis is given by 

where 

I (Bni+t^^— ■^) , 

f 

i; 

<n^jT(!,3_)|ii'> =^’ X )YmVFey¥T.lni-n'|=l. n=min (n^np 

1 0 ■ ( 5 . 21 ) 

- " 

III. 3b Two -Body Matrix Elements, <^P-iP-2 * 

Since most of the two-body interactions depend 
only on the relative coordinates of the two nucleons, it 
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is convenient to express tlie tsTO-body m-rtrix elements in 
terms of rel^.tive and cent re -of -mass coordinates of the 
two particles.. 'li-e following procedure Ib used, fhe two- 
hody uncoupled st.'tes 112 ^ 2 ^ 2 ^ 2 ^ 2 ^ 

are coupled to a total angular momentuiii J and total iso spin 
^ state with the use of 3-3 symbols; 

I^A^A'^1’ "'^2 '-'2 2 V 



TI16 'two-"P^2?*bicl0 couplsd 81 ^ 3,156 ill coupling sciioniB is 

iransf orined io ilie^v—S coupling sclieiiiB oy means of 9"“3 syffilols 

in the following ways 
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[ t. 




( WE^lEVJLop] ) < i 


^2 7i 1^ 

i S > X 


V 


3l 


3o J 






(3.23) 


Further, Moshinslry transformation is used to write 


n^l/^ngf/g^X^ = ^ <^ntl^L^jn^^^n2i2^^ j n^EL, 7\)> , 

nbl'TL 

(3.24) 

where the <^ nl^lTL Ajn^i^, ^n 2 ^ 2 )\^ Brody and Moshinsky 

transformation hrachets (Bry 60), Ihe quantum numbers (nt) 
and (ETL) specify the relative motion of 'the two particles 
and the motion of the centre-of-mass of tlae pair respectively, 
The total orbital angular momentum is the result of the 
coupling of and 1, A simple recoupling procedure yields 


jnt, BL, i i, S: 

= (-l)J+S+A^ ([>.][J])* 


I A 


jED, niSJ^: 

(3.25) 


The quantity inside the curly-b rackets is the Wigner 6-3 
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symtol. The final expression for the laost general antisym- 
metrized two-body matrix elements of an interaction Y which 
conserves total spin S and isospin 'f" is 


<nl 






= ^ ([3i][32K3il[32l)^ [J'lWCXlEslUlEriy, 

Mt 

M^nn' I, Iffii 

■ t 

X<nthL K|n^t,^n2t;2^> <. ^ n^|^n2L2 ^ . X 





( 3 - 26 ) 



where 


S 

nn'tt- 


<^niSi J^Ms r Ss J M,^'^ 


7572 ^ ^ V?; SJM*'y^ l\'^* 


(3.27) 


In most of the situations, the interaction 7(r) conserves 
total angular momentuitt, parity and isospin. The 3-j» 6-3 and 
9-3 symbols follow the notation of Edmonds (Eds 57). The 
normalization constant for the matrix elements is 


1 




Thus, apart from the geometrical factors occurring in the 
expression of equc^tion (5*26), the two-body matrix elements 
can be obtained in a straightforward manner from the relative 
matrix elements of 7 in an oscillator basis. 


111.3c Matrix-Elements of HE Single Particle 
Potential, TJ | ^ i 

Since the summation in equation (3-20) is restricted 
to the radial qu£..ntum numbers, only the matrix elements of 

the f ^2^2^2“2'^2| ’ “■2^2^ 2“'2'^2/^ 

are req.uired in the calculation of <(^ M’]l| ^ ‘ Tfriting 


explicitly. 







V 


D 


3 2^1 2"^ 2 

( occ) 


312112 




\. 


K-i^lhhVl^ '^ 2 ^ 232 “ 2 ^ 2 > 

(3.28) 


witli 


' “2“2 


/(*'-b3 


2 . 2 ) = H 

a( occ) 




G^/ ( 

^2 


3o) • (3.29) 


2^2 


In tlie case of a doulDly closed shell nucleus all the occupied 
32 “Sii'bshells are completely filled, i.e., each one of them 
contains 2(232+1) particles. Thus, the summation over m 2 and 
X 2 In eq.uation (3.26) can he carried out very easily and the 
orthogonality of the 3-3 symbols loads to the result 


2__ 

1112X2 

= 6(tiii.) 6(3 i3i) 


n2^2^2“^2’^2l'^i:Nl ^i^3]_)3i3_T^, 


(m^^m^) (239+1) I*(B-) 3-^ii-|ii^ > 




^2^2^2^2*’' 2')> 


(3.30) 



u:^ 




i: 

nn'i V NL 


(-)^+^'[j'][A][ X][S][j][r] 


X <X (/ UL /' i 1,2 ^ > 


^'1 ^‘2 Xl I ^a ^2 7 ^ ^ 


J i 


2 


I 3i 3; 


j i 


< i 


1 

2 


a ^2 


S >X 


J J 




J- s 


J 


>n1 


_/ y 


\ 


>L 

is 




I 


J 


(3.31) 


■/ i 


with. 


G’ 


j'r s 

nn^l i-"' 




1 

2 


jX s 


i 6(.3_X) 


T 


-M 


< 


1 

2 


Tl ^2. 


nn'il' 


-M 


T, 


(3.32) 


Thus. 


<njD(53_3i)|n^> , 


(3.33) 
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where 


IT 

n^^=l 


y~ Pn24^^232> ^ 232+1) 

^232(occ) 


P ( ^ ^ 2^2^2^ ' 

(3.34) 


For non-closed-sliell nuclei, which contain partially filled 
j-suhshells, following approximation is used for the density 
p of the occupied levels. It is assumed that the neutrons 
and protons are equivalent and averaging over the m-suh states 
for a given 3 ic done. The average density is given hy 




■C — ^a^^2^2^ 

I 2 ( 232 + 1 ) 

a( occ) 




cS/(i> 32 > 


n 


(3.35) 


where 1 ihoio) = number of nucleons in a given sub- 

sli©H ( ^2^2^ * 

Equations (3.29) to (3-33) remain the same. 


III . 4 Iterative Solutions of HE Equati_ojis 

With equations (3-21) and (3*33) s equation (3.14) 


becomes 



(3.36) 



where N is the number of terms included in expansion (3.20), 
There is one such equation for each (|L']_3]_) * fh-S total ground 
state energy defined in equation (3.17) is 

^0 " '^ZZ -Z - i n-,n:(^^l^l^ 2 ( 23 ^+!) 

^111^=1 ^33_(occ) " 


X [<^i|t(^)|4> + I<^]J u^^i3i)|X>3 


(3.37) 

The factor 2(23^+!) comes from the sum. on m^ and The set 

of equations (3.36) is solved by the following iteration 
procedure: 

(i) An initial set of C^EqOq) is guessed and also which 
orbits are occupied. Usually, ~ ^na' 

(ii) With this set of coefficients the matrix elements 
of HE Hamiltonian occurring in equation (3.36) are 
calculated. 

(iii) Equation (3.36) is solved by diagonalizing a H XN 

Hamiltonian matrix. A new set of coefficients 

G^'Ctn 3 i ) is obtained. Total ground state energy E^ 

1 1 

is computed. 



DO 


These new set of coefficients are recycled throu^ the step 
(ii) for recalculating the matrix elements and another diag— 
onalization of the matrix. This process is repeated until 
successive di agonal izations produce the same set of coef- 
ficients 0^ ^^ 3 ^) and reaches a constant value. The final 
set of » are the correct expansion coefficientsi 

is the total hinding energy of the system and e # . are the 
single particle levels of the nucleus. The calculations are 
repeated for sever^ values of oscillator range parameter 
h = ^(h/m^)), which enters in our expression through the 
expansion (3.20), The minimum in the E^ vs. h curve is 
determined. 

III. 5 Goulomh and Centre-of-Mass Correction s 

There are two important corrections to HE solutions 
arising due to the coulomh force between the protons and 
centre— of-mass motion of the. nucleus. The effect of these 
corrections on the total binding energy can be estimated 
externally if one assumes some simple model for the nucleus. 
Thus, for example, if the nucleus is assmed to have a uniform 
spherical charge distribution, its coulomb energy would be 
given by 

E = 2 ^ Z(Z-l) with R = 1.31^“^^ • (3.58) 

Similarly, if one assumes that the centre-of-mass of the 
nucleus is in the lowest Is state of a harmonic oscillator 



well, the kinetic energy of its motion is 


with-ha>= , (3.39) 

However, the effects of these corrections on various 
nuclear properties calculated in the framework of HH theory 
can he determined only if they are accounted for in a self- 
consistent manner, i.e., these ©orrections are included in 
the Hamiltonian h first and then the HP equations are solved. 
For this, the following procedure is used. 


III. 5a Ooulomh Correction 

The Goulomh force between the nucleons can be 
W£i«ten as 

7=™’- = (l+2i:p ( 1 + 21 : 2 ) . (3.40), 

.yCoul only if loili particles are protons# Tixc HajEiltonian 

12 

,H with the coulomb correction is 


H 


Xi. 




i:*:! i<i5: 






A A. 

■1 ■ r : it 

i=l 


A ft. 


(5.41) 


' ,new _ +T , y'" 

where 


GOUl 


( 3 . 42 ) 
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Thus, the only change which would occur in the EP eq.uations 

is that the two-hody matrix elements '^ik 

will he replaced hy those of These matrix elements 

can still he given hy expression (3.26), hut /, ,9 now will 

nn I L 

consist of two terms, 


1' 


S ^ S 


nh 1 9j 


(nucleus) + (coulomb) (3.43) 


with 


t: 


iJ T S 
nn'Jti' 


|-jn' t' SJM, TM^> 


(coulomb) = (ntSJM, T M, 

1 

6(Tl) 6(M^)6(itt^Js^^(r) 

(3.44) 


Care should be excercised in carrying the summation over Tg 
in equation (3.30) since coulomb matrix elements depend upon . 
the isospin quantum niimhers. However, the sum over ^2 
still he carried over explicitly if we make the time-saving, 
hut non-essential approximation that in calculating , 

is independent of T 2 * Sximmation over m 2 is dona as 
before. Thus, equation (3.31) will now be replaced 

by 


S 

^nn' tl' 


= (nuoleus) + (ooulcD) (3.45) 
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with 


nJ'T S 


(coulomb) = 2 



2 


2 



F 


J'T'S 
nn't i' 



c S 


(3.46) 


Subsequent changes will occur in the expressions that follow 
equation (3.30). 


Ill.Sh Centre-of-Mass Motion 

Considering a nucleus with A particles, we use 3A 
coordinates for A orbitals^ However, the centre-of-mass 
coordinates should be subtracted, so that we really have 3A-3 
independent coordinates. Thus, the Hamiltonian -sdiich describes 
the intrinsic state of the system is obtained by removing the 
centre-of-mass energy in equation (3«1). Thus, the intrinsic 
Hamiltonian is 



(3.47) 


with the total momentum ^ = r (3.48) 

i 

The operator of the oentre-of-mass Mnet^ energy can he 
writtsn as a smn of one— hohy and. two ••ho dy operators y 


' p2 __ y- P j^Pic 

'2m7A ~ 2m. A 


ik 
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The square of the relative 


Bioinenttim is 


(Pik)^ = T l'^. --^l2 


4 ‘ P i - ■^1, 


|((p|)^ + (^)^ - 2^-^) 


+ Pfe) - pj 


1/-2 2 t 2 


ik 


“■Sn^ 


1 


2_2s Pik 


4V - 


V ■ 


(3.49) 


Equation (3.47) no?; becomes 


Hx = 


7 


I “ 2 iT - ^ 2 i 


'■ ov.2 

^Pjk 

k AMj^ ’ 


(3-50) 


which is a pure tv;o-body operator. The secondpterm can be 

c m ^^12 

regarded as an additional potential T^p * = ■ and #an 
be treated as the nucleon-nucleon two-body interaction ’7^2" 

s 


This adds an additional term to E, 
and (3.43) » which is no?; given by 


nn'lt' 


in equations (3-26) 


(o.m.) = 6{tO| <n|l(t) |i/> , (3.51) 

where <n|T(t)tn') is defined in equation (3-21). The solution 
proceeds as iDefore ezcept that now in (3-36) and (3-37) one- 
body term (nj T( £3^) |n^ is missing. 



III. 6 Hartree-'f ocI l Density Distri'butioii 
and R.M«S. Eadii 

III. 6 a Density 


Ihe density of the particles inside the nucleus 


is 













( 3 . 52 ) 


Using expression (3.5) for (p obI integrating oter 
coordinates, one obtains 


(A-1) 


A 




(T^i 

(liF 


(3.53) 


0 , nnd o>s explicitly from equations (3.16) and 

Writing ana y 


( 3 , 18 )» and summi: 


iii,-' over m and x we obtain 


£1^(1)^ ,p -v 

nii^ t 3 

(3.54) 


le normallted so tiat it = A- 


Tliis densx’ty xs 



74 


and neutron charge densities, ^q(^) and ^(r), have the same 
form except that in equation (3.53) the summation is only 


over the proton or neutron orbits respectively, and the factor 

eqiiation (3.54) is replaced by ) • 

I'inally , the three radial densities are related as 

f,,(r) = ppr) + fpi) . (3.55) 

III. 6b R.M.S. Radius 


Ihe mean square radius for the matter distribution 
is given by 



Using equation (3.54) this becomes 


(3.56) 

(3.57) 


= I ^ (2j+l) (^d) 

nn^ 2 



(3.58) 


The integration can be carried over exactly and finally 
^^2^ = ( 23 + 1 ) [fnn(to) (2n+l^|-) 

n-tj 

" ?.^n n+l(^3) fn(n+i^4%). 


(3.59) 
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2 

The x.m.s. radius is '{ y . The mean square radius for 
charge distribution is given by expression similar to equation 
( 3 . 59 ) except that A is replaced by Z and svm is only over the 
proton orbits, 

III. 7 Hg Pote ntial 

The oscillator matrix elements of the HF single 
particle potentia,! U are obtained from equation (3.34) after 
several iterations have determined a self-consistent density 
matrix It is possible to express U in coordinate 

space by the equation, 


3r 


¥(r|rO = 


17 ^ 

^1^1“^ iVi 
( occ) 




(r) 


<P„/ 




(r^). (3.60) 


Using explicit forms for 9 ’s we have 


U(r|?0 


‘'■1 .“^1 


^l^l^l'^l 


^ 1 * ■^1 


(3t6l) 


where 

U 


(b 


,3 (r.rO -4r — — — 

^ ^ ( 3 . 62 ) 



76 


Tliis is a non— loc3,l and siaie— dependent potential* A possible 
local ec^uivalent to tins potential can be obtained in tbe 
static approximation which consists of computing the velocity- 
dependent potential U(r}^ corresponding to non-local TJ at 
0. Thus, 




= \ U ( ? j r / ) dr" 


/ 



^l4 



(r)r^ dr"' , 

(3.63) 


after doing the angular integrals. The integral on r^ can be 
done by using the Laplace transform on page 30 of Magnus and 
Oberhettinger (Mas 49) . The result is 


J E;Q(r) rdr = b5/2(|)lA 


[i(2n-l)!] 


1/2 


2^”^ (n-1)! 


(-) 


n-1 


(3.64) 


The static potential becomes, 


U(?|?= 0) = -^ <n^|U(0i)|n^) b5/2(|)l/4x 




A 1^-3 7: 

2 ^ (n^-1) ! 


(3.65) 
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Choice of I'Tuelei for Spherical HP QaloulatioiLS 

The HP Hamiltonian h does not necessarily have the 
same symmetry property as the nuclear Hamiltonian H. The 
possible symmetries of the HP Hamiltonian are governed by the 
following theorem (Ria 68). 

Theorem: Let E be an operator that commutes with H, 

ii^HR = H . (3.6€) 

Then, if E leaves the set of occupied orbitals invariant, it 
commutes with HP Hamiltonian h, 

E'^hE = h . (3*67) 


The operator E will leave the set of occupied orbitals invari- 
ant if the states obtained by operating E on any occupied 
orbitals a may be expressed as a linear combination of occupied 
orbitals only. That is, 

.A 

0 = 1 


a = 1 


( 3 . 68 ) 


If the nuclear Hamiltonian^ H has spherical ^mmetiy, it com- 
mutes with the three operators 2*'y the total 
angular momentum operator t- »e wish that HP Hamiltonian 
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should also have spiierical symmetry in vrtiioh case it ^ould 
commute with and The se operators will leave 

the set <a = of occupied orbitals invariant if for 

. _ J 

a particular j aJLl m-substates are either filled or empty. 

Hence, spherical s3rniiietry can exist only in nuclei that liave 
right number of particles to form closed shells. Such nuclei 
in the range A 4 4-0 are ^He, ^^0, ^^0 and "^^Ga and, thus, one 
would expect that results of spherical HI’ calculations would 
describe the ground state of these nuclei correctly. If the 
coulomb interaction is neglected and nuclear forces are assumed 
to be charge independent then the HF field is an isospin scalar, 
and neutron and proton orbits are degenerate. 


I I I . 9 Single Osci ll ator Approximatio n 
Limit to HE Solutions 


In the limit v/hen there is only one term in the 

expansion (3.20) , the HP orbitals are represented by pure 

oscillator functions. In this case P,, ,, f will become 

V ^ 2^-2 




2'^2 


z. 


a(occ) 


c“ g“/ 
^2 ^^2 


a( occ) 




= > 


a (occ) 


oci-^2 ^2^2 


(5.69) 


Thus, the density P,, ..f of the levels would be either one qt' 

*^ 2^2 
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zero depending upon whether the state Pg is occupied or unoc- 
cupied in the ground state configuration, fhe single particle 
potential matrix elements are 


<Pl[u[Pl^ 


n j! / ^ ^ 

r ' 2^''2 




= 7 


a(occ) 2 2 2 


21 , 


(3.70) 


^ 2 ( 000 ) 


The total energy is 


E. 


“21 f ^l[^ j^l) "^2 <^^lj^|^l>^ 


^i^i 


1 

2 


H [< ^‘ll ® l^i>, 1 < nl'f K> ^ > 


p^(occ) 


(3.71) 


where the potential energy of the nucleus is, 

= <*^1^2 (^1*^2 > (3*72) 

p^p.2(occ) 

for U = 2 nuclei with IS closed shells equation (3.72) reduces 
to a simple expression. 



P.E. = <n4sj|V|ntSJ> , 

str J 

where the weight-factor , ic given by 


(3.73) 




'npl 


/ 


[^] [J] X 


y~'' ?■'' ■( 

i-'s rv IJ \ ^ > 


^hti^2lz >■ 

(occ) 

The potential energy can also be written as 

P.E. = / Hill, \x) + Will, \l) , 

p(occ) P 

(occ) 


(3.74) 


(3.75) 


where W(p, ) is the contribution to potential energy due to 
interaction of nucleons in the state p, and can be obtained 
with the help of equation (3.72). For doubly closed shell 
nuclei, the single particle potential U, defined in equation 
(3.70), can be written down explicitly as, 

^^1 ^’1 ^ l ^ l "^ 1 ! ^ I ■^ 2'^2 ^ 2 ^ 2 ^ 2 ^ 


n^32^ 
(occ) 


<nA3l> natsSg: 




+(-l) 


|^2^2^2’ ^l^^l* 


(3.76) 



Equation (3.76) defines the single particle potential for 
neutrons. Eor protons, additional coulonh term should he 


added and the result is that Y in equation (3.76) is replaced 
hy V , defined in equation (3.42). I’he kinetic energy 1 of 
a single nucleon in state 1 is given hy 

i = (2n^ + ti - i) . (3.77) 


III. 10 


ration 


Ihe deterninantal v-rave function 


constructed from 


the final self-consistent HE' single particle orhitals cp^s 
cannot he the true nany-hody wave function and hence it is 
relevant to ask the question how good is With the pos- 
sible improvements to one can write an improved as 




(3.78) 


where 





Jph ^p^h 4'hE * 


(3.79) 


a)2P”2h 

THE 


— + + 

^Xp^P2:h^h2 ^Pi^P2%% YHE 


(3.80) 


;h^h2 


p and h refer to the HE unoccupied and occupied states 
respectively. Ihe various amplitudes, i.e,, Aq, 



Eq[uation (3.76) defines the single particle potential for 
neutrons. Eor protons, additional coulonh term should he 
added and the result is that Y in eq.uation (3.76) is replaced 
hy V , defined in equation (3.42). fhe kinetic energy T of 
a single nucleon in state T is given by 


r = (2n^ + - t) . (3.77) 


a 


III. 10 Yalidity of the HP Approximation 


Ihe deterninantal wave fimction 


$ 


EE 


constructed from 


the final self-consistent Hh single particle orbitals 9’s 
cannot be the true many-body wave function and hence it is 

“Y" 

relevant to ask the question how good is With the pos- 
sible improvements to one can write an improved as 


$1 = 


Aq 


'Hi- + 


$ 


Ip-lh 

HE 


+ A, 


, 2p-2h 


2T HE 


(3.78) 


where 


-Ih 



+ 


■^ph ®p^h H'hE 


( 3 . 79 ) 


m2p-2h 

THE 


4- + 




PpPg 

■h^2 


(3.80) 


P and h refer to the HE unoccupied and occupied states 

4^ph» 


respectively. Ihe various amplitudes, i.e., Aq, AtX_ 





'^2^P2_P2” principle calculable. If 1 q is very- 

large compared to tlie otlier amplitudes, is close to ^ j- 

Making use of the fact that the HP Hamiltonian h does not 


have any matrix elements bet-ween the occupied and unoccupied 


states, it can be slioYm that the amplitudes for the Ip-lh 
excitations are identically zero and, thus, the corrections 
io can come only from 2p - 2h and more complicated 

excitations, llie admixture of 2p-2h excitations to is 

measured by the ratio of some typical off-diagonal matrix 


elements of the liamiltonian H, between the wave functions 
corresponding to these excitations and to the smallest 

amount of energy needed for two particle excitations. A small 
ratio results in small admixture. Same is true for higher 
excitations'. This, then means, that the validity of the HP 
approximation is measured by the difference in -the HP single 
particle energies of the unoccupied and occupied states, The 
lowest of these differences is the one betv/een the uppermost 
occupied level aid the lower most unoccupied level. If this 
is significantly large compared to the differences in the 
energy of the -unoccupied states themselves or of the occupied 
states, we can spealc of the existence of a gap of energy. 

The larger this gap is, the nearer is 'to ^ Sometimes, 

a relaxation of the constraint imposed on single particle 
orbitals results in a considerable improvement in the 
approximation and is reflected in increase in the energy gap 
between occupied and unoccupied orbitals. 



Hi' approximation is a variation -.procedure as described 
earlier in Section III.l and provides an upper bound to tbe 
energy of tiie system, The possible corrections to arise 
from 2p-21i and higiier excitations, and can be calculated by means 
of perturbation theory in Hartree-Eocx basis. However, HE sol- 
utions have to be \mderstood in a different way when one is 
dealing with the effective Hamiltonian of Shahin, Waghmare and 
co-wortors. Those authors obtain an effective Ha mi ltonian by 
using the theory of cluster expansion and higher-order terms 
in this expansion are dropped, HE method to such an effective 
Hamiltonian has not been applied as a variation procedure since 
we are not dealing with complete Hamiltonian. We, then, have 
used the HE method as an attempt to find orbitals in this 
harmonic oscillator basis that yield self-consistency so that 
Ip-lh corrections to the energy, calculated with this effective 
interaction are minimized. 



CHAPTER lY 


APPLICATIOHS OP HlRTEEE-POCK TilEDEr MW DISCUSSION 

^"^•1 ' Details of ijhe Calculation 

Spherical Hartree-Pook (HP) - calculations have heen 

carried out in oonsiderable detail for the nuclei ^Be, 

12 IS 

G and 0 wi'lfi effective Yale interaction. Less complete 
calculations have .also 'been done for and ^^0. In view 

of the fact that besides the properties such as binding energy, 
root-mean-square radius, spin— orbit splittings, etc, , we are 
also interested in the distribution of protons inside a nucleus 
and that the motion of the centre-of-mass of the pair of part- 
icles contributes considerably to the binding energy, we add 
the coulomb and centre-of-mass correction terms to' Y while 
calculating the HP potential and wave functions (see Section 
III.5 for details). The effect of these tv;o terms modifies 
the Yv'ave functions through the density matrix. These effects 
have not been included in a self-consistent way by Shakin 
et. al. (Shn 67b) while calculating the properties of ^^0 
and '^^Ca with effective Yale interaction. However, while 
they have included the second-order term in pseudopotential 
(YP) by adding a term YP ^ YP to defined in equation 

(2,27), we have not included it in our calculation, HP 
calculations have also been done for the nuclei %e , ^^0 

and "^^Ca with Sussex interaction. In these calculations also. 



are included in 


coulomb and ceiitre-of-mass mairix elemeii'ts 
a self-consistent rianner in all the cases, except in 
for wbicii tlie effect of these terms is estimated externally. 
Iwo-body matrix elerients are calculated u,sing formula (3.26) , 
Since the nuclear force has a short-range it is necessary 
to use onl;/ a few relative ^ values ill il le calculation of 
two-body matrix elements of the nucleon-iiucleoh -interaction. 

A maximum of relative ‘t=2 states has been used in our calcul- 
ation. Higher H states are not expected to contribute signif- 
icantly for light nuclei studied in the present work (lan 68). 
The configuration space used in our El? calculations of A ^ 16 


nuclei consists of s^/2» P 3/2 /2 states of 

tt = 1, 2 and 3 magor shells for each of these states. In the 
calculation for ‘^^Ca and also for some of the lighter nuclei,- 


a larger configuration space is employed which includes 

and oscillator states of a = 1, 2 and 5 major shells in 

addition to states described -above. Thus, in all the HH 

calculations expansion (3.20) contains only three terms. 

Earlier HE calculations (Das 66: Sve 65: Shn 67 b: Krr 66) 

provide some justification for such a truncation. However, 

this truncation males' our results depend on b( = f ), 

A 

the oscillator range parameter and in the spirit of a proper 
variational calculation the criterion for the best choice of 
,b is the one which minimizes binding energy. After a suf- 
ficient number of iterations have been perfomed, the computed 



q^uan'ti’ties converge to corLstanl; values and tlie number of iter- 
ations needed for sucto. a convergence depend upon the initial 
0-coefficients. In all the calculations v/e start with initial 
and ten iterations have been found to give a 
satisfactory convergence of computed cuantities which finally 
describe the properties of the nucleus. In order to study the 
nature of the HF r/eve functions and also the single particle 
potential, we have concentrated mainly on the results obtained 
for ^^0. Ihe reason for this being, although there is quant- 
itative disagreement between the HP wave functions of various 
nuclei studied here, the general nature of these wave functions 
remains the same for all the nuclei- The discrepancies wher- 
ever encountered have been pointed out. Single oscillator 
configuration approximation to the HP solutions has also been 
discussed. ¥e present the results of otir calculations in the 
following sections. ■ 

Effective Yale Interaction Results 

IV. 2a Binding Energv', Single Particle Energies 
and Root-LIean-Square Radius 

The ground. state total energ;/' E^ of the nucleus is 
calculated as a function of oscillator range parameter b. The 
results of the binding energy as a function of b are shown in 
Pigure IV. 1 for the nuclei ^He, Be, and 0, Ideally, 
the HP results should be independent of b or basis chosen and 
vs. b curve should be a straight line. This would be true 




oo 


only when the III’ v/ave functions are expanded in a complete set. 
The larger the numToer of terms included in expansion ( 3. 20) , 
the weaker would he the dependence on b, Hov/ever, from the 
figure we notice that the curves are reasonably flat near the 
values of b which givo maximum binding, She two-bo{3y matrix 
elements needed for these calculations are evaluated using 
tables of relative matrix elements for the effective Tale 
potential given in reference (Shn 67b). 

In lable lY.l, we list calculated HP properties of 
various nuclei for the values of b which give maximum binding 
(b=1.76l for ^He and b=2.09I' for ®Be, and ^^0). Ihese 

properties are the binding energy per particle (B.E./A), 
single particle energies, matter and charge root-mean-sq.uare 
(r.m.s.) radius, and spin-orbit splittings {l-s splittings). 

Ihe corresponding experimental values wherever available 
are also given for the sake of comparison. Experimentally 
single particle neutron levels have been determined from the 
stripping process, as for example, the (d, p) reaction. 

Analysis of the experimental data on cross-sections for 
excitations of the various states of the final nucleus and 
angular distributions corresponding to these states gives 
information about the shell model single particle neutron 
states. Similarly, neutron hole states are located with 
(p, d) or (d, t) pick-up reactions while proton hole states 
are determined with the help of (p, 2p) and (e, ep) knock-out 



Oaloiilated Hi' properties with effective Yale interaction for various nuclei 
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reactions, Froii tlie TalDle IV. 1 we see that binding energy 
results are in poor agreement with, the experimental values. 

The reason for the discrepancy in the binding energies stemP 
from the fact tr^at we have not tahen into account the second- 

order term in T2. Calculations of Shahin et. al. (Shn 67b) 

show clearly that this term is i*)!5>ortant as far as the bind- 
ing energies of nuclei are concerned. In fact this term 
alone gives a binding of about 3 - 4 HeV per particle. In 
view of this, the theoretical results are in satisfactory 

agreement for ^He and ^^0 but still somevdiat poor for ^Be and 

12 

0. This is, hov/ever, not surprising. From the experimental 
information on transition rates and quadrupole moments these 
nuclei are knovai to be deformed. This conclusion is also 
reached from the calculations of Bassichis et. al. (Bas 67). 

In fact, these authors have made deformed HF calculations 
(mix:ing different orbital angular momentiua states) and obtain 

O ir\ 

binding energies for '^Be and which are lower than spher- 
ical HF estimates. Their results a.re in satisfactory agre- 
ement with the erperiment. Pal et. al. (Pal 67) have also 
recently performed similar calculations for the 8^ A ^40 
nuclei and obtain satisfactory agreement for the deformed 
nuclei in this region. However, these authors have not 
included coulomb .and centre-of-mass corrections in a self- 
consistent manner. On the other hand, if we include more 
orbitals in the exrpansion of basis states our results would 
be somewhat improved. - 
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llie calculated lPjy 2 “ ^^1/2 splitting in ^^0 is 
alout 10 percent snialler than tlie result obtained from quasi 
(ppO scattering (lyii 66), ihe agreement would be still 
better if second-order terns in VP are included (Slm 67b). 

Tile HP single particle levels are in satisfactory agreement 
with the experincntal values. There is a large gap between 
the occupied and the unoccupied levels in the case of ^He and 
^^0 and is a characteristic of HP approximation. The calcul- 
ated charge and matter r.m.s. radii are in fair agreement with 
the experimental data. When the second-order term in VP ia 
taken into account the particles inter ant more strongly and 
it is expected tha,t this would result in an increase in the 
binding of HP single particle levels and decrease in the 
values of charge and matter r.m.s. radii. 

In order to estimate the effects of coulomb (coul) 
and centre-of-mass (c.m.) corrections on various nuclear prop- 
erties we have repeated the HP calculations for the following 
cases. 

Case As Without coulomb and centre-of-mass corrections- 
Case B: With centre-of-mass and without coulomb corrections. 
Case Cs ?/ith coulomb and v/ithout centre-of-mass corrections. 
Case Bs With' both the corrections included- 
Table IV. 2 shows the total ground state energy E^ obtaihed 
in different cases for b=2.09P. The kinetic energy of the 
centre-of-mass motion is determined from the difference 



Effect of coulomb and centre-of-mass corrections on total ground state 
energy of various nuclei obtained with effective Yale inte raction, 

b = 2.0gP. iO-1 energies are in MeV, 
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1 



^etvreen obtained in cases D and G. Sinilarly, the coulomb 

energy of the nucleus is obtained fron the diffe^rence of total- 

energies in cases h and B. from the table, we sec that 

centre-of-mass correction is very important nuclei 

4 8 

He, Be v^rhile coulomb correction is not so important. For 
0, coulomb and centro-of-mass effects tend to cancel each 
other. These corrections to the binding energies have also 
been estimated enternally by assuming that the charge distr- 
ibution of the nucleus is that of a uniformly charged sphere 
and centre-of-mass moves in Is state of a harmonic oscillator 
well (see expressions 3.58 and 3.39) and compa^^ed with the 
values obtained self-consistently , It is se®^ that in all 
the cases the spherical model overestimates the coulomb energy. 


The effect of these corrections on the HE single 

particle energies is sho^m in Table I¥.3 for b=2.09E. The 

energy of the ls^y'2 state for ‘^He almost 

double when the centre-of-mass correction is included. This 

correction is also very important to reproduce the correct 

12 

ordering of Ip^^^o and ^'Di/2 Particle levels in G, 

In Table IV. 4, we give the matter r.n.s. (r^) charge 

r.m.s. (r ) radii obtained in different cases for b=2.09F. 

The coulomb correction increases r^ in all the nuclei, whereas 
the centre-of-mass correction reduces r^ in the case of He 
and ^^0 by approximately 20 percent and s’ percent respect- 
ively. Besides, charge r.m.s. radius is lat‘g®r than matter 
r.m.s'. radius ill all the cases. Effect of coulomb and 



Effect; of coulomb and oentre-of-mass corrections on HE single particle energies 
of various nuclei obtained with effective Yale interaction, b = 2.09E. 

j\ll energies are in MeV. 
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gable I Y.4 


Effect of oouloEb and centre-of-mass ccrrections on 
oharge (tg) and mass (r,^) r.si.s. radius of various 

U 



nuclei 

obtaine 

d v/itii ei 

'f ective 

Yale 




interaction, b = 

= 2.09E. 



Nucleus 

Case A 

u 

Case B 

To (J) 

Case 

ro(E) 

n 

(F) 

r„ 

D 

i-c (S') 


2.2229 

1.9013 

2.2382 

2.2438 

1.9066 

1.9096 

CD 

PQ 

00 

5.0335 

3.3531 

3.0873 

.3.1113 

3.4283 

3.4735 

12. 

^6 

2.7403 

2,7526 

2.79 26 

2.8144 

•2. 8158 

2.8571 

16 o 

^8 

2.5420 

2.515^^ 

2.5734 

2.5876 

2.5486 

2.5699 



centre-of-EiaSvS corrections on HE calculations lias also been 
recently studied by Gunye (Gue 68a) for several nuclei in the 
2s-ld shell. 

HE calculations for Be and have been carried 
out for only one value of b=2.09E. ilie caiculated binding 
energies per particle for ^®Be and are -0.77 MeY and 
-3.33 MeY comp8,red to the experimental values of -6.5 MeY 
and -7.5 MeY respectively. Ihus, ^'^Be is considerably under- 
bound under the assumption of a spherically symmetric HE field 
Our interest in doing the HE calculations for these nuclei is 
not only to study their ground state sta-tic properties but 
also to employ the HE single particle energies and wave 
functions obtained in these calculations in a configuration 
mixing calcula,tion described in a later Chapter. In view pf 
this, we shall onij' present the results of single particle 
energies and expansion coefficients for these nuclei. Ihe 
neutron single pa,rticle HE levels are given in Table IY.5. 

In the case of lPT_y '2 le'^el lies below the lp ^^2 level 

in absence of centre-of-mass correction, an effect similar 

12 

to that encountered in 0. Figure IY.2 shov/s the variation 
of neutron single r-)a,rticle HE levels with the nucleon number 
A. Proton energies v/ould show similar variation. As seen 
from the figure, the levels become more and more bound as 
the nucleon number increases and the effect is largest for 
the lOTj-est /o 


level 
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Table IV. 9 

HF single particle energies of various nuclei for 
effective Yale interaction obtained ?,’itliout 
coulomb and centre-of-raass corrections. 


Nucleus 

States 

(Me?) 

(!ieV) 

“Os 

(Me?) 

lSi/2 

-26.99 

-47.10 

-44.85 

IP3/2 

- 4.91 

-20.99 

-18.95 

lPl/2 

- 2.56 

-16.94 

-13.78 

ld5/2 

- 1.38 

— 0. 66 

0.67 

2V2 

2.50 

0.48 

1.59 

ld^/2 

- 

3.10 

4.70 


3.10 


4.70 


SINGLE PARTICLE ENERGY (MeV) 



fartatioa of HF nmtxxm particle ener^ ol 

with eff eotiTT# lale.' tHe nmleon 

aumber. A* 
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The single particle energies can also be obtained, 
directly by computing the self-consistent total energj’’ S 

©X 

of the A+1 (for unoccupied states) or 1-1 (for occupied states) 
system and talking the difference of E „ and E (the ground 
state energy of me nucleus mth A nucleons) t Such a calcul- 
ation has been carried out for ^^0 v/here a single particle 
was removed from the states ls^^ 2 » ^-' 3/2 added 


1/2 


2 Si /2 state in turn and it is assumed that resulting states 

15 17 

of 0 and 0 still remain spherical. Thiis, for the occupied 

16 *1:11 

levels of 0, the energy of the i level is given by 


= E(^^O) - E^(^^O) , 


(4.1) 


where i indicates a hole in the i level. The unoccupied 
levels are similarly given by 


ej_ = E^(^'^O) - E(^^O) 


(4.2) 


The difference between the single particle energies obtained 
in this way and the corresponding energies obtained in the 
HE calculation io the rearr^angemen-t energy of the system. 

The results are presented in Table IY.6, 

In a similar manner, particle-hole excited states 
of ^^0 are calculated e^ssuming that they a.re still spherical. 
Some such states are given in Table I?. 7. The first excited 
0^ state of ^^0 is experimentally observed at 6,05 MeV and 
is expected to arise from two particle-two hole (2p-2h), 
four particle-four hole (4p-4h) excitations whereas the 
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Table IY.6 


Comp ari so n of slngl e p arti cl e neut ron eM rgi e s 
obtained directly witii eigenvalues of ^*^0 for 
effective Yale interaction, b=2.09P. 

~ ± plus being for 

particle states and mini’s for hole states. 


State 

Total 
energy of 
the state 
E^CMeY) 

Excitation 

energy 

®ex (*’'5 

Eigen- 

value 

e^.(MeT) 

Rearrange- 
ment energy 

l^air^exl 

(MeY) 

G'.S* 

-57.835 

0.000 




-17,061 

-40.773 

-45.596 

4.822 

(lP3/2)"^ 

-38.734 

-19.101 

-21.208 

2.107 

(IPl/g)"^ 

-43.572 

-14.263 

-15.753 

1.490 

( ^2 ) 

-57.921 

- 0.086 

- 0.528 

0.442 
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Table IV. 7 

Particle-hole states in ^^0 for effective Yale inter- 
action^ 1=2.095. 3 is determined from the 

difference of total energies of excited 
state and ground state* 


State 

Total energy 

EQ^eV) 

Excitati on 

energy 

E (MeY) 
ex 

Ground State 

-57.335 

0.000 


(l%/2) 

-50.436 

7.399 

(lPl/ 2 ) 


-43.3S0 

13.455 


(l<i5/2)^ 

-42.100 

15.655 

(iPi/a)"® 


-30.651 

27.184 

(IS1/2) 

( ^2 ^ ^ 

13.669 

71.504 

(lSi/2) 

( 2 Pi/ 2 )^ 

13.894 

71.729 
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lowest calculated 0^ state comes at 15.65 HeY. Thus, it is 

necessairy to introduce departure from spherical symmetry 

or residual interaction to describe this state. In fact, 

it is an erperinental fact that all nuclei which are spherical 

in ground state, such as ^^0, exhibit deformation in their 

excited states and vice versa- Thus, the excited states of 
16 

0 would not be described well under the assiunption of a 

zero order spherical field. Kelson (Ken 65) and Bassichis 

and Eipka (Bas 65) have satisfactorily reproduced the excited 
16 

states of 0 by carrying out deformed HP calculations for 
the excited particle hole states and using projection of 
angular momentum techniq^ue. They have obtained lowest 0"^ 
state very near the experimental value. 

IV. 2b Hartree-Pock Wave Punctions and Density 

After a sufficient number of iterations have been 
performed the single particle wave functions are determ- 

ined from a final set of coefficients by substituting 
them into equation 3.20. These expansion coefficients are 
given in Tables (IY.8 - IY.13) for the miclei ^He, ®Be, 

^^0, ^*^Be and ^“^0 respectively. The results are presented 
for the values of b which give maximum binding (b=1.76P for 

^He, b=2,09P for all other nuclei). We shall discuss here 

16 

the nature of the wave functions for 0 alone and only make 
remarks on the structure of v/ave functions for other nuclei. 
The radial HP wave function (see equation 3*19 for 
detailed expression) for the ^3/2 %/2 
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Table IV.8 


Expansion ooeffioients C^(|,i) for b = 1.76E 

obtained- witii effective Yale interaction. 


ites 


neutron 



proton 


a \ 

1 

2 

3 

1 

2 

3 

1 

0.9684 

0.2224 

0.1129 

-0.9694 

-0.2185 

-0.1123 

r 2 

-0.1573 ^ 

0.8959 

-0.4156 

0.1527 

-0.8940 

0..4213 

3 

0.1936 

-0.5847 

-0.9025 

0.1924 

-0.3912 

-0.8999 

1 

0.8055 

-0.4828 

0.3435 

0.7939 

-0.4966 

0.3510 

5 2 

5 

0.5920 

0..6800 

-0.4325 

0.6078 

0.6658 

-0.4327 

3 

0.0247 

-0..5518 

-0.8336 

0.0188 

-0.5569 

-0.8304 ■ 

1 

0.7137 

-0.5986 

0.3636 

0.7042 

-0.6060 

0.3700 

3 2 

-0.6931 

-0.5286 

0.4901 

-0.7020 

-0.5164 

0.4905 

3 

0.1012 

0.6018 

0.7922 

0.1062 

0.6051 

0.7890 


Tab le IV. 9 

Expansion ooeffioients for ®Be . , b = 2.09P 

n if 

obtained with effeotiYe Yale interaction. 


states 

neutron 


proton 


'i ^ 


1 2 

3 

1 2 

3 


1 

0.9191 0.3255 

0.2220 

0.9204 0.3220 

0.2219 

0 i 

2 

-0.2294 0.9003 • 

-C.37OI 

-0.2216 0.8970 

-0.3825 


3 

0.3204 -0.2893 - 

-c .9021 

0.5222 -0.5028 

-0.8969 


1 

0.8930 -0.2663 

0.3629 

0.6740 -0.3059 

0.3775 

1 i 

2 

0.3953 0.8494 - 

-0.3497 

0.4401 0.8277 

-0.3482 


3 

0.2151 -0.4557 • 

-0.8637 

0.2060 -0.4705 

-0.8580 


1 

0.8025 -0.4665 

0.3721 

0.7816 -0.4898 

0.3862 


1 I 2 0.5921 0.699: 0.3997 0.6203 0.6750 -0.3994 

3 0.0739 -0.5410 -0.8373 0.0651 -0.5518 -0.8314 
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fable I?. 10 

Expansion ooeffioients for ^^Cg, b = 2.09E 

obtained with effective Yale interaction. 


states 


neutron 


proton 

i 3 

a'V 

1 

2 3 

1 

2 3 


1 

0.9357 

0.3140 0.1608 

0.9379 

0.3083 0.1594 

0 1 

2 

-0.2430 

0.9041 -0.3516 

-0.2294 

0.8952 -0.3820 


3 

0.2557 

-0.2899 -0.9 223 

0.2604 

-0.3217 -0.9103 



1 

0.9672 

0.1215 

0.2233 

0.9685 

0.0844 0.2344 

1 i 

2 

-0.0115 

0.8983 

-0 . 439 2 

0.0317 

0.8915 -0.4518 


3 

0.2539 

-0.4222 

-0.8702 

0.2471 

-0.4450 -0.8608 


1 

0.9667 

0.0854 

0.2414 

0.9655 

0.0453 0.2563 

1 1 2 

0.0296 

0.8991 

-0.4367 

0.0768 

0.8912 -0.4471 

3 

0.2543 

-0.4293 

-0 .0666 

0.2487 

-0.4513 -0.8570 
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■Table 17 . 11(a) 

Expansion coefficients G^(tD ) for b = 2.09E 

obtained x^/ith effective Yale interaction. 


states 


neutron 


proton 


1 j 

\ n 

1 

2 3 

1 

2 

3 


1 

0.9472 

0.3033 0.103S 

0.9497 

0.2965 

0.1006 

° 1 

2 

-0.2532 

0.9065 -C.33SO 

-0.2361 

0.8892 - 

-0.3919 


3 

, 0.1966 

-0.2^39 -0.9354 

0.2056 

-0.3484 - 

-0.9145 


1 

0.9282 

0.3214 0.1876 

0.9341 

0.3034 

0.1880 


2 

-0.2118 

0.8706 -0.4441 

-0.1857 

0.8629 

-0.4700 


3 

0.3060 

-0.3725 -0.8761 

0.3048 

-0.4042 

-0.8624 


1 

0.9491 

0.2434 0.1998 

0.9555 

0.2102 

0 . 2068 


2 

-0.1294 

0.8799 -0.4573 

-0.0889 

0.8741 

-0.4775 

3 

0.2872 

-0.4081 -0.8666 

0.2811 

-0.4379 

-0.8540 


1 

0.9486 

0.0570 0.3143 

0.9460 

0.0082 

0.3240 


2 I 2 0.0931 0.9165-0.3890 0.1242 0.9142-0.3858 

0.3024 -0.3983 -0.8660 I 0. 2994 -0. 4052 -0.8638 


3 


Table 17.11 (1)) 

• Expansion ooeffioients ) for ^^Og, b = 2.09E 

obtained in the HE oalcalation with effeotive 




Yale 

interaction w 

centre- of-maes 

1th cut coiO-omb 

coxrecjti ons. 

and 




1 

2 

3 



1 

0.9475 

0 . 30 47 

0.0969 

0 

1 

2 

2 

-0.3008 

G.9522 

-0.0529 



3 

,0.1084 

-0.0210 

-0.9939 



1 

0.9414 

0.30 32 

0.1476 

1 

3 

2 

2 

-0.2547 

0.9261 

-0.2782 



3 

0.2210 

-0 . 2243 

-0.9491 



1 

0.9570 

0.24-98 

0.1476 

1 

1 

2 

2 

-0.1953 

0.9509 

-0.3088 



3 

0.2146 

-0.2S66 

-0.9396 



1 

0.9789 

0.0528 

0.2017 

2 

5 

2 

2 

0.0407 

0.9560 

-0.3496 



3 

0.20 02 

-0.3504 

-0.9150 



3 


0.9390 

0,3302 

0.0960 


-0. 25 55 
0.8569 
-0,4476 


0.2301 

-0.5958 

-0.8891 


ge Me IV. 1 2 

Expansion coefficients for b = 2.09 E 

obtained with effective Yale inter.nction* 



states 

neutron 


proton 



\ 

a 

2 3 

1 

2 3 



1 

0.9339 0.3023 0.1909 

0.9363 

-0.3044 0,1751 

0 

1 

2 

2 

-0.2150 0.9014 -0.3758 

-0.2277 

0.9059 0.3570 



3 

0.2857 -0.3099 -0.9068 

0.2673 

-0.2944 -0.9175 


3 

1 

0.9548 -0.0711 0.2837 

0.9601 

-0.0621 0.2725 

1 

2 

0.1908 0.8911 -0.4118 

0.1760 

0.8917 “0.4169 



3 

0.2280 -0.4483 -0.3643 

0.2171 

-0.4483 -0.8671 



1 

0.9133 -0.2546 0.3178 

0.9472 

-0.1594 0.2780 

1 

1 

2 

2 

0.3759 0.8272 -0.4177 

-0.2708 

0.8621 -0.4283 



3 

0.1566 -0.5010 -0.3512 

0.1714 

-0.4810 -0.8598 



1 

0.7076 -0.5740 0.4121 

-0.7288 

-0.5548 0.4013 

2 

5 

2 

2 

-0.7066 -0.5660 0.4247 

0.6847 

0.5'988 -0.4155 



3 

0.0105 0.5917 0.8061 

0.0097 

-0.5776 -0.8163 
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Table IV.13(a) 

ExpaJision Goeffioients for b = 2. OS'S' 

obtained with effective Yale interaction. 



states 

— — — 

neutron 


proton 

1 

D 

X" 

1 

2 3 

1 

2 3 



1 

0.9420 

0.3185 0.1062 

0.9482 

0.3035 0.0942 

0 

1 

2 

2. ■ 

-0.2715 

O.9O88 -C.3I69 

-0.2525 

0.8996 -0.3564 


3 

0.1975 

-0.2697 -0.9 425 

0.1929 

-0.5124 -0.9296 



1 

0.9133 

0.3552 C.199O 

0.9199 

0.3446 0.1872 

1 

3 

"9 

2 

-0.2368 

0.8610 -0.4502 

-0.2256 

0.855'5 -0.4662 


(L 

3 

0.3313 

-0.3640 -C.8705 

0.3208 

-0.3866 -0.8647 



1 

0.9268 

0.3091 0.2133 

0.9361 

0. 2916 0.1966 

1 

1 

2 

-0.1850 

0.8700 -G.457I 

-0.1693 

0.8637 -0.4748 


c 

3 

0.3269 

-0.3842 -0.8635 

0 . 3083 

-0.4112 -0.8579 


1 0.9392 0.1429 0.3121 0.9411 0.1614 0.2972 

5 2 -0.0052 0.9151 “O. 403 O -0.0277 0.9127 -0.4078 

0.3432 -0.3770 -0.3603 0.3371 -0.3755 -0..8633 


3 



gable IV.l^Cb) 


Expansion coefficients 0^(0 j ) for ^^Og, b = 2.09E 
obtained in the liF calculation with effective 
Yale interaction without coulomb and 
centre-of-mass corrections. 


it 

j 


1 

2 

3 



1 

■0.9499 

0.2991 

0.0903 

O' 

1 

2 

2 

-0.2990 

0.9541 

-0.0152 


3 

0.0907 

0.0126 

-0.9958 



1 

0.9405 

0.3083 

0.1429 

1 

3 

2 

-0. 2651 

0.9287 

-0.2593 


2 

3 

0.2126 

-0.2026 

-0-9552 



1 

0.9520 

0.2701 

0.1442 

1 

1 

2 

-0.2212 

0.9324 

-0.2856 

2 

3 

0.2116 

-0 . 2400 

-0.9474 



1 • 

■ 0.9792 

0.0727 

0.1894 

2 

5 

2 

-0.0028 

0.9382 

-0.3460 


"2 

3 

0.2028 

-0.3383 

-0.9189 



1 

0.9658 

-0.1589 

0.2050 

2 

3 

0 

2 

0.2291 

0.8931 

-0.3872 


c. 

3 

0.1216 

-0.4209 

-0.8989 


states are plotted in figures (IV. 5 - IV. 5). Ilie proton wave 
functions are only slightly different from these as can he 
noticed from the erpansion coefficients for the neutron and 
protons. Pure harnonic oscillator v^ave functions have also 
been calculated for a value of the parameter b which gives 
the same r.m.s. radius as obtained from HP calculations, and "fiiesB 
are plotted in Figures (IV, 3 - IV, 5). from these figures we 
notice that the occupied orbitals for 0 can be well approx- 
imated by the harnonic oscillator functions while the struct- 
ure of the unoccupied orbitals is quite different. Ihe un- 
occupied orbitals are damped inside and have large amplitudes 
at large distances. This is true for all the three states. 
Sinilar effects are observed for the wave functions of other 
nuclei, The reason for a value of b different for Hf and 
harmonic oscillator wave functions for the same value of the ; 
r.m.s. radius is obvious from Tables (IV. 8 - IV.13). In other ' 

I 

words, strong mixing of various orbitals is responsible for 
this difference as well as the nature of the unoccupied 
orbitals, from figures (IV. 3 - IV. 5), we also see that there | 
is an overall shift of the peaks of the Hf wave functions with 
respect to the harmonic oscillator ones to larger values of 

! 

r. In the spherical Hf calculations one mixes different 
orbitals corresponding to the far-lying states (n mixing) 

These states have V7eaker binding and are associated with | 

larger separation from the centre of potential. The range : 

in which nucleons move is considerably large compared to I 

1 

pure states,. Therefore, the average separation of the 
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nucleins from the centre of the potential is increased, result- 
ing in the shift of peaks, I'or ^^0, using the HP values of r.Hi.s. 
radius the calculated value of h for the oscillator fTmctions 
turns out to he 1.699F* However, for this value of h the two 
wave functions do not overlap exactly. Best overlap is obtained 
for the tabiilated h values (Table 17,14) , Bor ^e nucleus, the 
calculated value of b is 2. 4241'. This large b indicates weak 

binding as can be verified from Table 17.1. 

■' ■ ' ■ ' ' ' ' ' ' . ' . ^ ' ■ 

ji 

■ 

HB single particle proton orbitals with and without the ; 
correction teriiis have been plotted in Figures (17.6 - 17.8) for 

1 fi A 

some of the single particle states in 0 and ^He. As seen from ' 
the Figures (17,6 - 17.7), the correction terms have vei^y little 
effect on occupied orbitals in 0 but the unoccupied orbitals | 

have been significantly affected and the centre-of-mass correction 

■ 

is much more important than the coulomb correction. Figua^e 17. 8f 

A 

shows the proton orbital in He with and without the centre-of~. 

. 

mass correction. 

■ ■ ■'■HF.' matte;!!,' and cliarge densities ' li^t :'miole3- ■ ( '.He, ,■ 

®Be, ^^0, ^^0) are shown in Figures (17.9 - 17.10) . Both the 

16 

matter and charge densities of 0 show a dip at the cantre. 

Such a dip is predicted from electron scattering experiments, 

17. 2c Hartree-Fock Potential 

The matrix elements ^nfUjn*^ of the HF single particli 

T'6 

potential U are given in Table 17.15 for -0 neutron states for 
b=2.09F. Since < n| U fn ^ = < n^ , oiily n ^ n' values are 



116 


lalle IV«14 

Yalues of the oscillator parameter (h) giving 
the hest overlap of the H? orbital . for ^^00 
obtained vfith effective Yale interaction and 
harmonic oscillator wave functions. 


^bital 


b(F) 


ls^y '2 1-612 

1 p ^ y '2 1-625 

lPl/2 


1-667 




0*0 


2-0 An 3-0 4.0 5-0 6-0 7*0 

\\ ^ 


n 

n 

A‘A A 

Vw 

w V 




A — *5,y y 

// / 

/X 

■f c 

// I 


V \ // .> 

\vV// 


F1G.GZ:.6- state Ef protoa wave fanotions for obtaiaot ni% 

funn??nnffl OU*T» k OOTTSSpOnds t© WWTO 

opjloab ^ i>«btr«-bf-®a8S mxxmttimM- 

SfS ?L respsetivtly, Cw® » i® obtaSid 

bom m®^ ©©rroetionB iam aoooimt* , ■ 



G.JS'7- p 2^2 e-fe«t;e Kf=,.. proton wav® maotions for ^^0 obtaintS with 
offdOtiT® Tal» iat»raotl®a#-'' Other detail® are earn ae ia 
fie, if.i. 



effectiye Tele 'interaction.* Otixer details are saa® as in 
Fig. I?*6, 


•, i- 
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Ta ble IY.15 


Matrix 
for ^^0 


elements <fn |U({? 3 ){ of tiie in' potential 
g, b = 2. OSS’ obtained with effeoti-ve Yale 
interaction for neutron states. 



n^ ^ 

1 

2 

3 


1 

-23.5188 



0 1/2 

2 

-13-6123 

13.3058 



3 

- 6.4951 

1.6423 

25.7619 


1 

- 0.1265 



1 3/2 

2 

- 7.5401 

15.6957 



3 

- 7.1378 

3.2100 

25.5869 


1 

3-9921 



1 1/2 

2 

- 5.2082 

17.8884 



3 

- 6.5334 

4.4390 

26.8043 


1 

10.0815 



2 5/2 

2 

- 1.8058 

16.9826 



3 

, - 5.4234 

4.4253 

24.0448 



given. The potentials also been evaluated using 

expression (5,62) for different single particle neutron states 
(l,j) of ^^0 for b=2.09I’. These potentials have been plotted in 
Figures (lY.ll - IV .13) for the ^*3/2 ^1/2 

states respectively in ^^0. The fact that HF potential is non- 1 
local and stsite-dependent is obvious from the plots. A local 
potential would give a delta function at r = rC Potentials for 
^3/2 ^1/2 very similar. The HF single particle ■ 

potentials for other nuclei are expected to show similar behaviour 4 


The sta.tic limit of the HF potential has been calculated 
using equation (3.65). As seen from Figure IV. 14, this potent- 
ial does not reserable any of the standard local potential forms 
such as harmonic oscillator, Woods-Saxon, etc . The static limit 
of the potential is most attractive at about r = 1.8F and becomes 
repulsive beyond r = 3.6F. It is not clear what meaning this ^ 
local potential can have. 


We have also attempted to find out the parameters of 
a local potential ¥/hich would give the same single particle 
energies as obtained in the HF calculation. The local potent- 
ial employed is of Woods-Saxon form with spin-orbit and coulomb 
terms added to it. The detailed shnpe of the potential is given 
in Appendix A. Calculations have been carried out for 0. The 
range parameters r^^ and r^ of the potential were kept fixed at 

"I ' 

a value of 1,25F, appropriate for 0 and the strengths V^^ and 
V of the central and spin-orbit parts of the potential are 




FIG . M-i 1 
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varied so as to reproduce the HF single particle levels. For 
the lp^^ 2 ( 21.21 MeY) and lPn^2^15-75 MeY) neutron levels of 
the strengths are found to he Y,r = 58.206 MeY and Y = 8*854 

SO 

MeY, ihe corresponding single particle wave functions are 
plotted in Figures (IY.15 - IY,16) and also compared with the 
HF wave functions. It is seen that the two wave functions (HF 
and equivalent local potential) agree pretty closely in the 
region r43.5ijhut the HF wave function has a much longer tail. 
Hon-local correction to the wave function (see Appendix A) has 
also been applied. She Ip^y^ functions for different 

values of the non-locality range parameter are shown in 
Figure IY.16, As the parameter is increased the tail of the 
wave function approaches the HF wave function, but the disagre- 
ement in the inner region increases. 

It shoxiLd be mentioned that the two-body coulomb and 
centre-of-mass matrix elements, used in the calcxilations reported 
in this section, have been calculated xirith a maximum relative 
1,= 2. This has caused some error in our calculations since 
coulomb force is a long-range force and there would be signif- 
icant contributions coming from larger h values. Later on, 
coulomb and oentre-of-mass matrix elements were calculated with 
a maximum relative i = 10 (maximum allowed I in the smaller 
configuration space) and some of the HF calculations for effect- 
ive Tale interaction were repeated with iliese new coulomb and 
centre-of-mass matrix elements, A few revised HF results have 
been added to Appendix B, It is seen that static properties of 





H 
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nuclei are not very much affected. Ihese new coulomh and centre- 
of-mass matrix elements have been used in the HF calciilations 
with Sussex interaction reported in next section. 

IV. 3 Sussex Interaction Eesults 

HF calculations with Sussex relative matrix elements 
have been done for the nuclei "^He, , ^^0 and Iwo- 

body matrix elements occurring in the HF equations are calculated 
using tables of relative matrix elements given in reference 
(sit 68a),. The tables given in reference (Elt 68a) are not 
complete since only the matrix elements between the states 
whose oscillator energies [i.e. (2n+|— t)^^ differ by not more j 
than 2hUJ are tabulated. Some other matrix elements, such as 
off-diagonal matrix elements of the type <^n'fv|n y when the I 
radial quantum numbers differ by more than one, i.e. (n -n)^2, ! 
are also needed in the present work and have been calculated i 
using a simple recurrence formula (Bit 68a) from the tables I 

given by Elliott et. al. (see Appendix G) . Ihe complete set : 

of matrix elements used in otir calculation i c given in fables 
IV. 16 for some of the relative states, fhe specific form of 
the two-body Sussex interaction is not known but is not needed . 
in any part of the calculation. | 

Results of HF calculations obtained with Sussex inter- I 

I 

action are tabulated in fables (IV. 17 - IV.19) for the nuclei ; 

12q ^^0 and for several values of b. We have not | 



T able IV.l6(a) 


Sussex relative matrix elements for ^S^ channel, 
h = 2. OF. Only n^ ^ n values are given since 

Yjn^S^>= <n^aj_j ¥|rj 


n 

n' 1 

2 

3 ■■ 

4 

5 6 7 

1 

-3.58 





2 

-5.22 

-5.11 




3 

-4.41 

-4.24 - 

3.76 



4 

-3.45 

-3.55 - 

2.93 

-2.53 


5 

-2.47 

-2.46 - 

2.21 

-1.89 

-1.61 

6 

-1.56 

-1.62 - 

1.53 

-1.34 

-1.13 -0.89 

7 

-0.74 

-0.87 - 

0.89 

-0.83 

-0.67 -0.43 -0.07 




Table 

IV.16(b) 


Sussesf relative matrix 

element 

s for channel, 


h = 2. OF. Only 

11^4 n 

values 

are given since 



< n' J V ! 

in\> 

11 

txi 


IX 

n* 1 

2 

5 

4 

5 6 7 

1 

1.20 





2 

1.34 

1.85 




3 

1.40 

1.97 

2.48 



4 

1.41 

2.13 

2.70 

3.14 


5 

1.42 

2.27 

2.87 

3.26 

3.39 

6 

1.47 

2,36 

2.95 

3.20 

3.22 3.21 

7 

1.53 

2.40 

2.88 

3.00 

3.04 0.00 0.00 


gaMa 17 . 16 (g) 


Sussex relative matrix elements for ciiaimel, 
b = 2. CP. Only n'^n values are given since 
<n* ^P^ i V in^P^> = <^n^p^ I Y |n' ^P^>. 


il 

V 1 

2 

J 

4 

5 

6 

7 

1 

0.86 







2 

1.21 

1.78 






5 

1.46 

2.17 

2.71 





4 

1.66 

2.48 

3.09 

3.54 




5 

1.81 

2.72 

3.37 

3.84 

4.14 



6 

1.93 

2.88 

3.56 

4.01 

4.30 

4.50 


7 

2.02 

2.99 

3.65 

4.08 

4.39 

4.63 

0.00 


fable IY.16(d) 


5 

Sussex relative niatrix elements for Dp channel, 
b = 2. OP. Only n4 n values are given sinoe 
<^n' ^D2| Vjn^Dg) = 


n 

n ' 1 

2 

3 

4 

5 

6 7 

1 

- 1.63 






2 

- 1.43 

- 2.01 





3 

- 1.19 

- 2.11 

- 2.46 




4 

- 1.09 

- 1.92 

- 2.33 

- 2.60 



5 

- 1.04 

- 1.59 

- 2.15 

- 2.44 

- 2.74 


6 

- 0.93 

- 1.34 

- 1.87 

- 2.34 

- 2.62 

— 2.86 

7 

- 0.77 

- 1.16 

- 1.62 

- 2.15 

- 2.51 

1 

o 

» 

o 

0 

1 

o 

• 

o 

o 
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T&hle IV.17 


Oaloulated HI properties for obtained with Sussex 

interaction. 


b(P) 

lo 

1.8 

2.0 

2.2 

Experi- 

mental 

B.E./a (MeV) 

- 3.602 

- 3.263 

- 3.134 

- 2.947 

- 7.07^ 

Matter r.m.s, 
radius (P) 

1.695 

1.827 

1.098 

1.999 


Charge r.m.s. 
radius ’ (P) 

1.699 

1.031 

1.901 

2.002 

1.67^ 

Single particle 
energies (Me?) 






neutron 

-24.05 

-21.24 

-20.15 

-18.75 


1 S-| 

proton 

-23.21 

-20.47 

-19.41 

-18.04 

-20.4+0.3° 

neutron 

6.38 

4.62 

3.94 

3.42 

0.96*^ 

Ip, 

proton 

7.36. 

5, 46 

4.72 

4.12 

1.97*^ 

neutron 

8.86 

6.05 

4.99 

4.18 

4.36'^ 

IPl 

proton 

9.76 

6.83 

5.70 

4.81 

5.37'^ 

neutron 

proton 

13.53 

8.65 

6,83 

5.57 


14.47 

9.40 

7.50 

6.18 


neutron 

Ip^—lp^ 

■5 2 

2.49 

1.43 

1.05 

0.76 


splitting 
(Me?) proton 

2.41 

1.36 

0.98 

0.68 



a. Mattauch et, al. (Mah 65) 

b. Hofstadter et, al, (Her 67) 

c. Rlou (Mu 65) 

d. Yashahidze et. al-, (Yae 63) 
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Table IV. 18 


Galoulated 

b(p) 

B.E./a (MeV) 

Matter r.m.s. 
radius (p) 

Okarge r.m.s. 
radius (P) 

Single particle 
energies (MeY) 


HP properties for obtained with Sussex 

inijG raction . 


2.309 - 2.107 - 1.992 - 1.828 
2.227 2.404 2.489 2.609 


Experi- 

mental 

- 7.68^ 


2.238 2.419 2.504 2.625 2.47 


^®1 

neutron 

proton 

-47.73 

-42.03 

-39.67 

-36.51 

2 

-44.67 

-39.25 

-36.99 

-33.94 

Ip^ 

2 

neutron 

-15.34 

-13 . 49 

-12.79 

- 11.90 

proton 

-12.54 

- 10.92 

- 10.30 

- 9.52 


neutron 

- 11.98 

- 11.78 

-11.77 

-11.47 

IPl 




Ju 

2 

proton 

- 8.84 

- 8,85 

- 8.91 

- 9.23 


neutron 

5.04 

2.76 

2.34 

2.18 

2 Si 

2 

proton 

7.71 

5.00 

4.49 

4.17 

IP^-lPl 

neutron 

3.36 

1.72 

1.02 

0.43 

2 2 






splitting 

(MeV) 

proton 

3.70 

2.07 

1.39 

0.29 


-34.2+2^ 

-33.0+3.5 


a. Mattauoh et. al. (Mah 65) 

b. Hofstadter et. al. (Hor 67} 

c. Eiou (RLu 65) 


Calculated HI' properties for 

o 

interaction. 


obtained with Sussex 


b(p) 

1.5 

1.8 

2.0 

2.2 

Experi- 






mental 

B.E./l (lev) 

- 5.111 

- 4.910 

- 4.850 

- 4.772 

- 7.98^ 

Matter r. 

in* s* 

2.217 

2.325 

2.380 

2.478 


radius r^ 
0 

(E) 




Charge r. 

in*s. 

2.227 

2.337 

2.391 

2.479 

2.71^ 

radius r^ 

(p) 




r X B.E./a 

11.231 

11.416 

11.543 

11.825 

11° 

(MeV-P) 







Single particle 






energies 

(lev) 







neutron 

-57.16 

-53.05 

-51.53 

• 

00 

1 

tA -51 

2 

proton 

-53.04 

-49.16 

-47.51 

-44.65 

-34+3.5® 




-44+2 

IP3 

2 

neutron 

-28.82 

-26.57 

-25.77 

-24.44 

— 21.81*^ 

proton 

-24.86 

-22.82 

-22.11 

-21.17 

-18+2.5® 

-19+1 

IPl 

2 

neutron 

-19.90 

-19 .18 

-19 . 20 

-18.90 

-15. 65*^ 

proton 

-16.06 

-15.53 

-15.61 

-16.44 

-13+2° 

-12.4+1 


neutron 

2.12 

0.15 

0.10 

0.31 

- 4.15^ 

2s-|^ 

2 

proton 

5.76 

3.38 

5.15 

3 .12 


H 

1 

neutron 

8.92 

7.59 

6.57 

5.54 

6.16'^ 

? 2 







spli tting 
(MeV) 

proton 

8.80 

7.29 

6.50 

4.73 



a. 

b. 

c. 

d. 

e. 


Mattauch et. al. (Mah 65) 
Hofstadt«ti' et, al, (Hor 67/ 
Angeli et. al, (Ani 69) 

Cohen et. al, (Con 65 a» 

Eiou (Eiu 65) 


63h) 


! 


13T 


presented the results for h 
close to h = 1.5P results. 


- 1.4S’ and h = 1,6F, since they are 
Binding energy per particle vs, h 


curves (Figure I?,l7) are very nearly strai^t lines for all 

these nuclei ( He, “ C and ^^0). We find that the calculated 

binding energies for 0 is about (50 - 40 percent) less than 

the corresponding experimental value and still poorer agreement 

with the experimental data has been obtained in the case of ^He 
12 

and 0., As described in Section IV. 2a, binding energy result 
12 . 

for 0 will improve if the constraint of spherical symmetry 
is relaxed. Moreover, since minimas in the binding energy vs. 
b curves could not be obtained till b = 1.5F, the calculations 
should have been extended to lower values of b in order to know 
how much maximum binding one can obtain in such a calculation, 
Ihe matrix elements in reference (lit 68n) have been calculated 
only up to b = 1.4F since lack of experimental phase shift data 
above energies ^ 300 MeV imposes a lower limit on the value of 
b up to which matrix elements can be determined accurately (see 
Section 11,3). However, as pointed out by Elliott et, al, the 
matrix elements for lower b values can be obtained by a 
polynomial . fit to the matrix element vs. b curve in the known 
region and extrapolation of the fitted polynomial to smaller. b 
values, but the reliability of the matrix elements so obtained 
can be questioned. Moreover, even for the smaller configuration 
space included in the present calculation and for b^l.4i', the 
relative matrix elements of V involving enersies E'^ 300 MeV will 
also enter. Ihe second and higher-order terms (see expansion 
2.56) for' such matrix elements are quite important since the 
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auxiliaF/^ potential can be chosen to approximate the true potent- 
ial only below 300 KeY. In view of this, the calculations have 
not been carried out for smaller values of b. 

As seen from I'able IY..19, the sin,;^le particle energy 
levels in 0 are overbound compared to experimental data- fhe 
coulomb correction shifts the proton levels relative to neutron 
levels but the relative separation between any two proton levels 
is very nearly the same as the corresponding separation between 
the neutron levels. Angeli et. al. (Ani 68), on the basis of 
experimentally measured cross-sections, ha.ve suggested the fol- 
lowing relationship bet?/een the binding energy'" per particle and 
matter r.m.s, radius for 0, . 

B.E./A (Me¥) X r^(P)^ 11. (4.3) 

16 

Our HP calculations for 0 agree with this. 

Assuming tlia,t the Sussex interaction V can be broken 
up into central, spin-orbit and tensor parts, the matrix elements 
of any particular component can be expressed in terms of the 
matrix elements of the full potential o.s described in Appendix 0. 
Such a calculation has been carried out and the results are pres- 
ented in Table IY.20 for some of the relative states of nucleon- 
nucleon pair. In order to test the effects of each of these 
parts on various ground state properties hP calculations have 
jcepeated with the contributions of 
(i) only the central part. 
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Sable lY 

. 2C( a) 



Central, 

spin- orb it end 

tensor components of ■ 

Sussex 

relative 

matrix 

elements^’ n{V|n)>, 

diagonal in 

osci- 

11 at or radial quantmi 

niimber n 

for b = 2. OP. 

Radial Gh.ajmel 

Quantnm 

Number 

Central 

Spin- orbit Sensor 

Sotal 


u 

-0.13 

0.60 

-2.05 

-1.58 

n,=l 

^^1 

-0.13 

0.30 

1.03 

1 . 20 



-0.13 

-0.30 

-0.21 

-0.64 


5p 

^0 

-0.15 

1.34 

-2.67 

-1.48 

CM 

11 

'5p 

^1 

-0.15 

0.67 

1.33 

1.85 



-0.15 

-0.67 

-0.27 

-1.09 



-0.03 

2.11 

-2.92 

-0.84 

n=3 


-0.03 

1.06 

1.46 

2.49 



-0.C3 

-1.06 

-0.29 

-1.38 



-0.39 

0.10 

1.27 

0.98 

n=l 


-0 . 39 

0 . 03 

-1.27 

-1.63 



-0.39 

-0.07 

0 . 36 

-0.10 


Contd. . 
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Tabl e I V. 20(a). (Contd. ) 


Radi al 
Quanirum 
Rymber 

ChanrLel 

Oentral 

Spin- orbit 

Tensor 

Total 



-0.46 

0.28 

1.65 

1.47 

11 


-0 . 46 

0.09 

-1.65 

-2.02 



-0 . 46 

-0.19 

0.47 

-0.18 


\ 

-0.56 

0.43 

2.05 

1.92 

n=3 

\ 

-0 . 56 

0.14 

-2.05 

-2.47 


S 

-0.56 

-0.29 

0.58 

-0.27 
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gable I V. 20(l3) 


Central 

, spin-orbit and 

tensor components of 

Sussex 

relative matrix 

elements n |T|n > 

, off-diagonal in 


oscillator radial 

iu an turn 

number (n?^') 




for b 

= 2. OP. 



Radial 

Channel 

Central 

Spin- orbit Tensor 

Total 

QuaJitnm 

Number 








- 0.15 

0.88 

-2.1 

-1.37 

n = 1 

^^1 

- 0.15 

0.44 

1.05 

1.34 

n'= 2 

^2 

- 0.15 

-0.44 

-0.21 

-0.80 



-0.10 

1.68 

-2.46 

-0.88 

II 


-0.10 

C.84 

1 . 23 

1.97 

n' = 3 







-0.10 

-0.84 

-0.25 

- 1.19 


^0 

-0.13 

1.11 

- 1.95 

- 0.97 

11 

H 





n' = 3 

3p 

^1 

-0,13 

0.55 

0.98 

1.40 



-0.13 

-0.55 

-0.20 

-0.88 



- 0.51 

0.22 

1.19 

1.10 

11 il 

H 

«■ 

-0.31 

0.07 

-1.19 

-1.43 


D 3 

-0.31 

- 0.15 

0.34 

—0,12 


C ont d 
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S able I Y-. 20(13) (Contd. ) 


Radial 

Quantum 

Number 

Ghaimel 

Central 

Spin-orbit 

Sensor 

Total 



-C.48 

C.32 

1.73 

1-57 

n = 2 


-0.48 

0-11 

-1.73 

-2.11 

n = 3 







3 

-0.48 

-0.21 

0.50 

-0.19 



-0.25 

0.25 

1.03 

1.03 

n = 1 


-0.25 

0.08 

-1.03 

-1.20 

li 

'2; 






s 

-0.25 

-0.17 

0.29 

-0.13 


25 


-0.17 


0.29 


-0.13 
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(ii) central 'infi spin-orbit part, 

(iii) central end tensor parts of tie interaction. 

We present the results in Tables (IY.21 -■ 17.22) for the nuclei 

16 12 

0 and C respectively .and for b = 2.01. This is also the 
value of D for -c— .cl spin-orbit splitting for 0 is closest to 
enperimenta.l value (see Table 17.19). As seen from the tables, 
major portion of t'le binding energi' cones from central part of 
the interaction, v.lile the tensor force leads to saturation. In 
the ca,se of ^^0 rliiost all of the Epiii-orbit splitting comes 
from the spin-orjit part of the force. In G, tensor force 
also produces large splitting but of v/rong sign- This is cons- j 
istent with the f'j.ndings of Wong (Wog S3) and Tarbutton (Tan 68)' 
The contribution to the splittings from the spin-saturated ‘ ! 

shoj-ls (in spin~satur.ated shells both tie members of ‘file doubletl 
j = L+ are connletely filled), sucl as in 0, is positive and; 

I' 

comes mainly froii the spin-orbit part of the two-body potential.! 

1 2 , 

In the nuclei witli spin unsaturated s'-ells, such as in C, i 

I 

there is significant contribution to ‘lie splitting from the 
tensor forces, but it always has tie wrong sign. Thus, the 
effect of only partially filling up a j = t ± i doublet is to 
decrease the splitting. The spin-orbit splitting is a sensitive 
function of b and inoreases rapidly v/ith decreasing b. This is! 
eicpected since ‘lie spin-orbit force dominates at higlt-energy . | 

40 ' 

HF calculations have also been carried out for Ga ^ 
and also for ^^0 in the larger coi-ifiguration space which } 





Sable IV. 21 


Iffeot of central, spin-orbit and tensor parts of 
Sussex interaction on calculated Hf properties of 

1 C 




for 

O 

b = 2.0?. 




Only Central 

No Tensor 

N 0 Spin- orbi t 


(MeV) 

-4.828 

-4.924 

-4.917 

Matter 

radius 

(F) 


2.376 

2.365 

2.362 

Charge 

radius 

W 


2.387 

2.367 

2.365 


Single particle 
energies (MeV) 

Neutron 

Proton Neutron Proton 

NeutTOn Proton 

1S3_ 

2 

-51.11 

-47.29 

-51.78 

-47.82 

-51.61 -47.70 

IP3 

-23.60 

-19.96 

-26.13 

-22.63 

-23.86 -20.40 

2 






^Pl 

-23.60 

-19 .96 

-19.15 

-16.25 

- 23.73 - 20.95 

2 






2s^ 

0.15 

3.20 

0.14 

3.18 

0.20 3.21 

2 






IP3-IP1 

2 2 
splitting 
(MeV) 

0.00 

0.00 

6.98 

6.38 

0..13 - 0.53 
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l able IV, 22 

Effect of central, spin-orbit and tensor parts of 
Sussex interaction on calculated HE properties of 

for b = 2. OP. 


Only Central Ho Tensor Ho Spin-orbit 


(MeV) 

-1.506 


-2.627 


- 1.007 

Matter r.m. s, 

radius 

(E) 

2.594 


2.403 


2.712 

Charge r.m.s. 

radius 

(E) 

2.610 


2.415 


2.731 

Single particle 
energies (MeV) 

leutron Proton 

Neutron Proton 

Neutron Proton 

isi 

-34.25 - 

'31.72 

- 41.49 - 

•38.71 

-32.32 - 29.90 

2 






IP 3 

-11.67 - 

■ 9.30 

- 15.51 - 

■12.93 

- 9.34 - 7.07 

2 






IPl 

-10.49 - 

■ 7.72 

- 8.74 - 

• 6.13 

-12.54 - 10.09 

2 







2.63 

4.73 

2. 45 

• 4 . , 'G 2 

2.51 4.58 

2 






Ip^-lPp 

2 2 
splitting 
(MeV) 

1.18 

1.58 

6.77 

6 #80 

-3.20 -3.02 
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includes a^d d^^g oscillator states of n = 1 - 3 major 

shells as v/ell. In these calculations, coulomb and oentre-of- 
mass corrections were not included in a self-consistent manner. 
Results are presented in lables (IY.23 - IY.24) for '^^a and 
°0 respectively and for b = 2. OF. In the case of the 

calculated binding energy is close to experimental value, but 
the single particle levels are somewhat overbound compared to 
experimental data. Expansion coefficients are given in Sables 
(IY.25 - IY.28) for the nuclei ^He, ^*^0 and ^*^Ga respect- 

ively. Sable IY.29 sliov/s variation of expansion coefficients 
for ^^0 v/ith b, 

d'7.4 Discussion 

In our calculations with effective Yale interaction 

4 16 

we find that binding energies for He and 0 are in fair 
agreement with the experimental data when the second-order 

term in YP is fallen into account. Reasonable values for the 

3 12.,' 

binding energies of Be and C can only be obtained if the 
HP field is permitted to be deformed. An important correction 
to the binding energ;>^ arises because of tlie approximate 
treatment of the energy denominator e in the evaluation of 
second-order terms in expression (2,27). She energy denom- 
inator e is defined as the differ^ce in the energy between 
the excited intermediate state and the occupied initial state. 
In reference (Shn G7a) , the intermediate states were taken to 
be plane waves and the energy of the interacting pair was 
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lable IV.23 


/A 

Galoulated HI’ properties for ^ ob tallied vri-th. 

Sussex interaction mthout ooulonb and centre-of — 
mass corrections, b = 2. OF. 



Sussex 

Experimental 

Single particle energies 
(MeV) 

^®l/2 

-93.04 


IP3/2 

-61.70 


^^1/2 

-55.03 


l%/2 

-55.17 

■-21.9^ 

CM 

CM 

-25.88 

-18.2^ 

ld-/2 

-22.20 

-15.6^ 

B.E./A (MeV).. 

- 9-.84 


SooulA (MeV* 

1.85 


Bq.ji/A {MeV.^ 

- 0.20 


B.E./a (MeV) including c.m. 

- 8.19 

- 8.55^ 

and coul. corrections 



lP 3 / 2 “lPi /2 splitting (lleV) 

6,67 


^^5/2”^ *^3/2 splitting (MeV) 

10.97 

7-.0° 

iL- splitting in the same major 

- 2.90 

- 1.1^ 

shell, 2s-ld (MeV) 

J.m.s. radius (F) 

2.789 

3 . 50 


a. Elton et, al. (Eln 67) 

b, Mattauoh et. al. (Mah. 65) 

G. Oohen et. al. (Gon 63a) 

d. Dey et, al. (Dey-P re print) 

e. Hofstadter et. al. (Hor 67) 
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gable IV. 24 

Single particle energies of ^^Og, b = 2. OP obtained 
in the HF cal oi'lation with Sussex interaction 
mthout coulonb and oentre-of-mass oorreotions« 


State Single particle energy 

(MeV) 


Is^ 

-51.08 

2 


IP3 

-23.41 

2 


IPl 

-17.13 

2 


Idc 

5 

2 

- 1.87 

2si 

- 0.68 

2 


Id, 

3.40 


2 
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laol e IV. 25 

Expansion coefficients G^Ctj) for obtained 




VO. th Sussex into 

^raotionj b = 1.8P. 

states 


neutron 

proton 


Xu 

,i\ 

1 

2 

312 3 


1 

0.9700 

0.2152 

C.1123 0.9711 0.2108 0.1121 

0 — 

u 2 

2 

-0.1599 

0 . 915 c - 

• 0.3704 - 0.1547 0.9131 -0.3773 


3 

0.1829 

- 0.3413 - 

• 0.9220 0.1818 -0.3490 - 0.9193 


1 

3 

2 

1 

2 

3 

0.9020 

■ ■ 0.4301 

0.0371 

-0.3627 

0 . 8018 

- 0.4750 

0.2341 

- 0.4150 

-0.8792 

0.8939 

0.4472 

0.0316 

-0.5785 

0.7907 

-0.4812 

0.2401 

-0.4182 

-0.8761 



1 

0.8323 

-0.4949 

0.2495 

0.8239 

- 0.5057 

0 . 2558 

1 

1 

2 

2 

-0.5502 

-0.6G33 

0.4300 

-0.5621 

-0.6718 

0.4823 



3 

-0.0671 

0.5368 

0 

H 

CO 

• 

( 

0.0721 

0.5412 

0.8378 
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ar able IV. 26 

Expansion coefficients C^(ij) for obtained 

I'dth Sussex interaction, b = 2.0E, 



state s 


neutron 


proton 


1 

3 

X" 

1 

2 

3 

1 

2 

3 ■ 



1 

0.9397 

0.3139 

0.1359 

0.9422 

0.3075 

0.1330 

0 

1 

2 

2 

-0.3004 

0.9473 

-0.1109 

-0 . 2888 

0.9467 - 

•0.1427 



3 

0.1635 

- 0 . 0654 

-0.9845 

0.1698 

-0.0960 - 

-0.98(B 



1 

0.9647 

0 . 2140 

0.1532 

0.9688 

0.1954 

0.1525 

1 

5 

2 

2 

-0.1582 

0.9365 

-0.3126 

-0.1360 

0.9335 - 

•0.3318 



3 

0.2104 

-0.2773 

-0.9374 

0.2072 

-0.3007 - 

•0.9309 



1 

0.9675 

0.1999 

0.1547 

0.9714 

0.1801 

0.1546 

1 

1 

2 

2 

-0.1451 

0.94-04 

-0.3074 

-0.1220 

0.9376 - 

-0.3258 



3 

0.2069 

-0.2750 

-0.9389 

0.2036 

-0.2976 - 

■0.9327 
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lab l3 lY. 27(a) 

Expansion ooefficients j ) for obtained 

with Sussex interaction, b = 2.0E. 


states 

neutron 


proton 

h j 

\ 

a 

y n 1 2 3 

\ 

1 

2 3 

0 1 

1 

2 

3 

0.9433 0.3140 0.1079 

0.3163 -0.94S7 -0.0040 

-0.1011 -0.0379 0.9942 

0-9461 

-0.3040 

0.1121 

0.3069 0.1038 

0.0517 -0.0433 

-0.0094 -0.99 3 7 


1 

0.9236 0.3447 0.1677 

0.9286 

0.3328 0.1639 

1 i 

2 

-0.3026 0.9242 -0.2331 

-0.2837 

0.9 218 -0.2643 


3 

0.2353 -0.1645 -0.9579 

0.2391 

-0.1989-0.9504 



1 

0.9445 

0.2882 

0.1577 

0.9500 

0.2712 0.1549 

1 1 

2 

-0.2380 

0.9511 

-0.2765 

-0.2153 

0.9279 -0.3043 


3 

0.2265 

-0.2236 

-0.9480 

0.2262 

-0.2557 -0.9399 
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lal le I Y. 27(b) 


Expansion coefficients for obtained in 

the HE caloiilation with Sussex interaction without 
coulomb and centrc-cf-iiiass correction, h = 2. OP. 


state 

t 2 

a\". 
a \ 

1 

2 

5 

0 1 

1 

2 

3 

0,9398 

0.3275 

-0.0976 

0.3224 

-0.9444 

-0.0640 

0.1132 

-0.0287 

0.9932 


1 

0.9336 

0.3216 

0.1582 

^ i 

2 

-0.2817 

0.9314 

-0.2305 


3 

0.2215 

-0.1706 

-0.9601 


1 

0.9535 

0 . 2609 

0.1483 


2 

-0.2136 

0.9374 

-0.2749 


3 

0. 2107 

-0.2306 

-0.9500 


1 

0.9750 

0.1027 

0.1971 

1 

2 

-0.0325 

0.9428 

-0.3318 


3 

0.2139 

-0.3172 

-0.9225 


1 

0.9592 

-0.1930 

0.2065 

IV) 

lulvij 

2 

0.2G1S ' 

0.3821 

-0.3917 


3 

O.lObu 

-0.4298 

-0.8966 
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T able IV. 28 

Expansion eoefficients C^(y ) for '^^Ca2o obtained in 
the HP oalculation v/ith Sussex interaction without 
ooulomb and oentre-of-mass correction, b = 2.GP. 


states 

\ ^ 





h i 

1 

2 

3 



a \ 






1 

0.9455 

0.3139 

0.0870 


O 

ro|H 

2 

0.3211 

-0.8553 

-0.4108 



3 

-0.0547 

-0.4163 

0.9076 



1 

0.9163 

0.3781 

0.1320 


1 1 

2 

0.3995 

-0.8859 

-0. 2357 



3 

-0.0278 

-0.2687 

0.9628 



1 

0.9167 

0.3740 

0.1407 


^ 1 

2 

0 . 39 27 

-0.9082 

-0.1445 



3 

-0.0738 

-0.1878 

0.9794 


5 

2 2 

1 

0.9058 

0.3886 

0.1687 


2 

-0.3779 

0.9212 

-0.0926 



3 

C.1914 

-0.0201 

-0.9813 



1 

0.9516 

0.-3266 

0.1595 


2 ^ 

ci 2 

2 

-0.2938 

0.9350 

-0.1986 ■ 



3 

0.2140 

-0.1381 

-0.9670 




155 


Table I Y.29 

Expansion coefficients for “"Og obtained with Sussex 
interaction. Only G^(g, j ^ coefficient s 

are given. 


1 

D 

1d(E)— 

coefficient 

1.5 

1.8 

2.0 

2.2 



1 

-0 . .. , ^4 

0.97^ 

0.9433 

0.9162 

0 

1 

C ^ 

^2 

-0 .9653 

0.9708 

-0.9487 

-0.9199 




-0.9674 

-0.9936 

0.9942 

0.9855 




-0.9-J35 

0.9659 

0.9236 

0.8827 

1 

3 

2 

g2 

-0.9243 

0.9373 

0.9242 

0.9063 



^3 

-0.9247 

-0.9544 

-0.9579 

-0.9510 



or 

0.9985 

-0.9809 

0.9445 

0 9043 

1 

1 

2 

Cg 

0.90’ 2 

-0.9357 

0.9311 

0.9170 



c| 

-0.907S 

-0.9395 

-0.9480 

-0.9452 
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replaced loy a ooiistsiii; average value A, . 


i’iius, e was tcilen as 






(4.4) 


T/ith.A = 20 Lev". aoT/ever, A is actually state -dependent and 
generally larger tnan 20 MeV. llie approximate treatment of A 
causes an uncerte,r;.ity of about 20 percent in the matrix 

elements which are the most important ones in the binding energy 
calculation. Actuallj’’, the energy deno;:-:inator should be 


-e. 


m 


2J2^U^ + f 


(4.5) 


where and refer to the HP single particle energies 

of the initial state of the nucleon pair. Shis would mean that 

the second-order terms occurring in the effective interaction 

depend upon the HP solutions, which rn turn cannot be obtained 

?/ithout a pre-Knowledge of the effective interaction. Thus, 

there is a double self-consistency protlem involved. Shakin 

et. al. (Shn S7b) have solved this problem by approximating 

+ E with A. With an intelligent choice of A , e(A ) 

IHP ZaP 

can be made close to and the difference -^)V can be 

tjr ejjp e 

evaluated by perturbation theory in }IP basis. Ihese corrections 
are referred to as spectral corrections and are expected to add 
roughly per particle to the binding energies (Shn 67b). 
Ihus, binding enersr results would further improve when these 
corrections are tahen into account. 
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It is expected that the results v;ould improve if more 
terms are included in tne expansion (5.20) of HP wave functions. 
However, the situation is somewhat ambiguous due to the following 
reason. It has been pointed out by Pal and Stamp (Pal 67) that 
there is a possibility of double coiinting of the contribution 
of high-energy states v/hen HP wave functions are expanded in 
terms of oscillator states with large number of nodes and 
singular interactions are used, fhe strongly repulsive short- 
range part of these interactions gives rise to excitations to 
high-lying intermediate states. Phe contributions of these 
states are included in the calculation' of reaction matrix. How, 
v^hen the HP wave functions are expanded in terms of harmonic 
oscillator states of n = 1, 2, 3 and still higher major shells, 
which are quite high-lying in energy, there is a possibility of 
double counting, Ihus, 'one is left with two alternatives. Pirst 
is, if one wishes to avoid the double counting only low-lying 
oscillator states should be included in the expansion of HP 
orbitals and one lias to be satisfied with .almost pure oscillator 
states. This would not determine the tail of the wave function 
correctly. On the other hand, if more terms are used in the 
expansion (3»20) , li'xX’O HP wave functions are determined more 
acourately but there is always a danger of double counting. 

Since inclusion of oscillator states of n = 1, 2 and 3 major 
shells appears to be a reasonable compromise between these two 
alternatives, we expect that such a double counting would not 
cause any serious error in our calculations with effective Yale 


interaction. 
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In the calculation of Sussex natrix elements (Slt 63a ) , 

as pointed out oy iiavromatis (Mas 69) ^ the higher-order terms in 

the auxiliary potential are left out (see expansion 2.56)» l.he 

most important of these are the second— order Bom term V V and 

e 

the dispersion term t.('^- — §")'t^fi* I^^ ’tne calculation of reaction 

A On e^ A 

matrix elements vrith Hamada-Johnston (HJ) potential Kuo and Brown 
(Kuo 66) have observed that exclusion principle plays an important 
role in reducing the size of the Born term for singlet states. 
However, they have also observed that there is a large contrib- 
ution to this term from the tensor part of the force for triplet 
states. The studies of Dahlblom (Dam 64) and levinger et. al. 

(ler 60) also show that tensor forces give large contributions 
to binding energy when calculated to second-order in perturb- 
ation theory; tiieir contribution being typically ten times larger 
than the corresponding contribution for the pure central forces. 
Thus, there would be considerable improvement in the binding 
energy results when this term is accounted for. Dispersion 
corrections are also expected to contribute significantly. 

In Tables (IY.30 - IV. 33) ^^0 and results of 

present calculation are compared with those of various other 
spherical HB and Bruedmer-Hartree-Bock (BHP) calculations. 
Different types of nucleon-nucleon interactions used in these 
calculations are also specified. It is seen that Yale and HJ 
potentials give very similar results when the second-order term 
in VP is not taken into accoiait in the former. This is hardly 
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Table ly.^p 

Oompariscn of our results for the HI? proton single partlole energies of 

with, other oaloulations# jill nuinbers are in MeV. 




' ■ Table IV, 33 

Oompari&.«in of our results for the HF single particle energies in ^^CiagQ with 
other calculations. Numbers inside rounded brackets are the results obtained 
when seccnd-'order tem in pseudopotential was included in effective Yale inter- 
action. Proton energies are marked with 'f , iCLl numbers are in MeY. 
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su3!prising becaut'e of the similar nature of these t\TO potentials 
and also fox tne lact that they give same quantitative fit to 
the two— body data. Besides, in the calculations referred in 
the table, reaction matrix elements from both these potentials 
were obtained using similar techniques. For example, effective 
interaction froBi Yale potential v/as derived using Moszkowski- 
Scott separation technique and an attractive pseudopotential 
was added for those states where the potential is repulsive. 
Reaction matrix elements from HJ potential were also calculated 
using separation technique and the reference-spectrum method was 
used for the repulsive part of the potential. It is of some 
interest to compare the Sussex interaction results with those 
of potentials obtained from a fit to the phase shift data such 
as HJ, Yale^and Sabakin. It is found that Sussex interaction 
results for the binding energy are better than any other real- 
istic interaction results at least in first-order. Strikingly 
close similarity can be seen in the Sussex and effective Yale 
results when the second-order term in ¥1 is included in latter. 

Very few spherical IIP and BHP calculations have been carried out 
19 A 

for G and He. In the spherical HP calculations of Davies 
et. al. (Das 66) using a very simple velocity-dependent potent- 
ial the calculated B.B./A is -2.40 MeV, which is larger than 
the values obtained with Yale (—1.42 '"'eV) and Sussex (-1.99 MeV) 
interactions (see Sables IV.l and IV, 18). Recently, Krieger 
(Krr-Preprint) has reported HP calculations in a basis of 
cartesian harmonic oscillator wave functions using a velocity- 
dependent interaction especially derived fox use in HP 
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calculations. For %e, lie ottains, B.S./A = -0.825 Me¥ coiapared 
to our calculated values of -3.77 fceY and -3.60 Me? with effect- 
ive Yale and Sussex interactions respectively. Elius, it appears 
that realistic interactions whicli fit the two-hody data give mors 
binding in Sle than the effective interaction of Krieger which 
fits the nuclear matter data- I 


?fe have found that the binding of the lowest Is 


1/2 


state increases very rapidly as the nucleon number changes* A 
local potential with reasonable strength would not give such 
deeply bound states and thus, this phenomenon is attributed 

to the non-locality of the HF potential. There is some experim- ; 

i; 

ental evidence from (e, e^y) experiments (Ami 64) that the - 

binding energy of this level really changes from ij:; 35 to 6'0 Me? i 
as one goes from A = 12 to 28. However, the interpretation of | 
the data is not unambiguous in view of the very wide peak cor- : 


responding to devel. More experimental evidence to settle ; 


this point conclusively seeDBto be necessary. 


HF is the first term in Brueclnier-Goldstone perturb- 
ation expansibh end second and higher-order corrections to it 
will arise frbm 2pi2h (two particle-two hole), 4p-4h, and hi^er 
excitations. The large gap between the occupied and the un- 
occupied levels indicate that these c'irrections would be small 
so that HF solutions would be «lose to the actual solutions. 
However, in order to obtain good wave functions and also 

■ ■ ■ i 

: improved -biading ener^, one must introduce 2p-2h type 
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correlations in tlie fraraeworls: of the theory. The error- 
introduced in choosing the harmonic oscillator wave functions 
as a basis set can be considerably reduced by improving the ' 
wave functions at large distances to give proper tails. The 
results of our analyses also indicate tioat our approximation 
of assiiming an average value of the oscillator parameter b (same 
for all states) nay not be good for malring detailed calculations 
with HP wave functions. However, the oasic characteristics 
indicate that care must be excercised vviiile treating the 
particle-hole types of problems. 

IY,5 Results for Sing--e Oscillator 
Configuration Approximation 

In this approximation, we have calculated the binding 
energies of ^He, ^^C, ^^0 and ^^Ga. In the case when HP wave 

t: 

functions are replaced by pure oseillator wave functions, the 
binding energies of doubly closed LS-shell nuclei (ne, 0, | 

^^Ga) are calculated using expression (3«73) . Wei^t-factors [ 

I 

for ^^0 and ^^Ga a.re given in Tables (IY.54 I¥.35). has [ 

a partially filled p-shell and we have used the following 
approximation to calculate its potential energy. 

P.B.(^^G) = ¥(ls. Is) + I [¥(ls, Ip) + W(lp, lp)l , ; 

(4.6) j 

' ■' ■■ . .i 

where ¥(a, §) gives the contribution to potential energy arising! 

from interaction of nucleons in the states u and p. Binding 

■ , , ■ , ■ ■ ■ ■ ' 

results for these nuclei are presented in Table IY.36 


energy 
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J_a61e I,Y.34 

1 c 

Wei ght-fao tors for Og. We liave used the notation 

2X^1 f 2S+i for relative states. 

J 


State 

Weigh t-f actor 

^^S^(OOOO) 

3 

^^Sq(OOOO) 

3 

^^S^(OOOl) 

9 

^^Sq(OOOI) 

9 

^^S3_(0010)] 

^^S^(0002) 

i 

1 

X 

/ 

^^SgCOOlO)"] 

9 

5^Sq(0002) 

> 

^^1^(0101.) 

6 

^^?q{010L) 

6 

^^P^(OIOL) 

18 

^^PgCoioi) 

50 

^^1)^(0200) 

1.5 

^^BgCOSOO) 

2.5 

^^D^(0200) 

■ 3.5 

^^D2(0200) 

7.5 

^5g^(l000) 

1,5 

5^Sq(1000) 

: 1-5 
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gable 17.35 


Weight-factors for lotation used is same 

as in Table 17.34, Weight-factors for states 

are same as those for states. 


ible I 


State 



^^S3_(0000) 

3 

^^S3_(0001) 

9 

^^S^(0002) 

15 

^^s^(ooio) 

3.0 

^^S^(OOll) 

6.75 

^^8^(0003) 

15.75 

^^S^( 0020 ) 

1.125 

^^S^( 0012 ) 

5.62 

^^8^(0004) 

10.13 

^^ 8 ^( 1000 ) 

3-0 

^^ 8 ^( 1001 ) 

5.25 

^^ 8 ^( 1010 ) 

1.75 

^^ 8 ^( 1002 ) 

3.13 

^^ 8 ^( 2000 ) 

11.12 

^:^p^(oioo) 

3 

^^Pq(OIOO) 

3 

^^P^(OIOO) 

9 

^^p 2 ( 0100 ) 

15 

^^P^(OIOI) 

9 

^^PqCOIOI) 

9 


^^P3_(0101) 

27 

^^PgCOlOl) 

45 

^^P^(0102) 

7.75 

^^Pq(0102) 

7.75 

^^P3_(0102) 

26.25 

^^P2(0102) 

43.75 

^^P^(OllO) 

1.75 

^^Pq(OIIO) 

1.75 

^^P^(OllO) 

5.25 

^^P 2 (0110) 

8.75 

^^P^(Olll) 

2. 25 

^^Pq(OIII) 

2.25 

^^Pt_(0111) 

6.75 

^^12(0111) 

11.25 

^^p^(0103) 

5.25 

^^?q^0103) 

5.25 

5\(0103) 

10.75 

5^P2(0103) 

26 . 25 

^^P^(llOO) 

2.25 

5%q(1100) 

12.25 



gable IY.35 (Contd.) 


State 

3.-eiglit-factor 

Sl'a'.te 

TTeigiit-factor 

^^P^dlGC} 

6.75 

15 , 

03(0210) 

1 m 4 S 


11,25 

^^Pp(C210) 

3.13 

(1101) 

2.25 

13i:_(120C) 

1*12 

53p (1101) 

2.25 

13d2(1200) 

1.38 

53 Pd 1101) 

6.75 

1%)3(1200) 

2.62 

j)3p / 1 1 pi 

Ju. -L. w *mL- j 

11.25 

3132(1200) 

5.62 

■^^11 (0200) 

3 

‘3r-3(C3CO) 

5.25 

^0200) 

5 

^3'o(0300) 

11.25 

13 

D^(0200) 

7 

33 

0300) 

15.75 

31 1 


33 , 


D2(0200) 

15 

10300) 

20.25 

(0201) 

5.25 

^^3(0301) 

5.25 

^3^2 (0201) 

8.75 

33jV:,(05C1) 

11.25 

^3d^(0201: 

12.24 

33:'„(0301) 

15.75 

31 1^(0 201) 

26. 24 

51^(0301) 

20.25 

^^P-(02G2) 

3. OS 

■^^d(0400) 

2.63 

13 '*’ 

K-(020 2) 

5.13 

13 

64(0400) 

3.32 

13d^(020 2) 

7.20 

l%p(C40C) 

4.15 

31d2(02C-) 

15.42 

3d ,(0400) 

10.13 

13dt_(0210) 

0.62 



^^D2^021C) 

1.04 





Table 1Y.36 


Binding energy per particle of Tarious nuclei 
obtained in the single oscillator 
configuration approximation. 


Interaction 


Sussex 


Yale 

nucleus 


12. 



^8 

b(B) 

1.4 

1.8 

1.8 

2.0 

1.76 

B.l./A (MeV) 

0.33 

-3.72 

-4.17 

-7.30 

-3.14. 
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for the values of "b which, give r.m.s, radius close to experin- 
ental values. Results obtained with effective Xale interaction 
are also given .XXX. this table. In the binding energy calculatior 
of Oa with Sussex interaction and for b = 2. OP we find that 
contributions to potential energy of S, F, D and F states are 
93 • 4-^ s 0,’^'jy , 5.3*^ and 0.55^ respectively. The G-^state 
contributes alnoet" negligibly, P-state natrix elements are 
large but they tend to cmcel each others contribution to 
potential energ;^ . On the other hand, although P-state relative 
.matrix elements are relatively smaller, they occur with large 
wei^t-faotors and make a net contribution to potential energy 
comparable to there of P states. Thus, in any precise calcul- 
ations P-state should be included while G— state can he neglecte4 

Single particle neutron and proton energies of 0 

obtained with Sussex and Yale interaction are given in Table 

IY.37 (see Section III. 9 for formulae used). It is seen that 

Sussex levels are more spread out for analler values of b. Singl«f 

40 

particle neutron energies for Ca are given in Table IV. 38 alor^ 
with the experimental values. On comparing these energies with 
self-consistent la? single particles energies (see Table IV. 23), 
it is seen that the effect of self-consistency on the single 
particle levels is to increase the binding of the occupied 
states. Same coiicltision can be drawi on comparing the 0 leved 
g,iven in Table IV. 37 with tliose given earlier (see Tables IV. 1 
and IV. 19) . 



Table IV. 37 

Single particle energy levels of 
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Table IV. 3S 


Neutron single 

particle energy 

for b = 2.0 E. 

levels of 

single 
parti ole 
state 

SUSSEX 

{UeV) 

Experimental^ 

Id^ 

0 

-21.8 

-21.9 

2si 

2 

1 

ro 

« 

O 

-18.2 

Id, 

3 

2 

-lS-85 

-15.6 

“7 

2 

- 5.06 

- 8.4 

CM 

- 6.C3 

- 6.3 

1% 

-6.99 

- 4.3 

2p. 

- 3.07 

- 2.9 




a. Elton et. al. 

(Sin 67) 



173 


On comparing the single oscillator configuration 
approximation results with those obtained in HF calculations, 
we find that the approximation of using pure oscillator wave 
functions is not too had as far as the estimation of binding 
energies and single particle energies are concerned. 
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On comparing the single oscillator configuration 
approximation results with those obtained in HP calculations, 
we find that the approximation of using pure oscillator wave 
functions is not too bad as far as the estimation of binding 
energies and single particle energies are concerned. 



GHAP!EBR V 


aOMPIGURATION MIXING CllOUiAlIOIS WIIE 
liAEIEEE-KJCK WA¥1 FUirailOIS 

V.l 33etails of tlie Galoulations 

In this Chapter we shall describe the results of 
configuration mixing calcilL at ions, with Hartree-Foclc (HE) wave 
fxmctions^ carried out in the spirit of a standard shell model 
spectroscopy calculation for the nuclei and ^^Be. The 

main asstunption of the shell model is that the nucleons move 
independently of each other in an average potential well , com- 
monly taken to he of harmonic oscillator or Woods-Saxbn type. ' 
The more fundamental way of obtaining the average po'tential ! 

field of the shell model, starting from free nucleon-nucleon 
interaction, is the HP method. ¥e have seen in Chapter lY that [ 

i 

the HP orbitals for the occupied levels of low-mass nuclei can I 

be well approximated by pure oscillator wave functions but the | 

i 

two differ appreciably for the unoccupied levels. Since the HP i 
orbitals are obtained in a more fundamental way, one would be 
tempted to believe that they should be closer to the true single I 
nucleon wave functions compared to the pure oscillator wave i 

fxmctions. The single particle wave functions and energies | 

used in our' calculations have been obtained by solving the HP I 

equations seif-consistently as described in Chapters III and lY. 
Each unknown HP sin^e particle orbital is expanded in terms | 
of pure oscillator v/ave functions 

9a * 


(5.1) 
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Here a and !i specify the angular and isospin quantum numbers of 
the Hartree-Foch (HF) and harmonic oscillator (HO) single 
particle states; a = (atfjmx), p = (nijm-r) , where a and n are 
radial quantuia numbers. She HI field is constrained to be 
spherically symmetric (see Chapter III)# and also the forces are 
assumed to be charge independent so tliat the summation in equ- 
ation (5.1) is restricted only to the radial quantum number n. 
Ihree oscillator states of u =1, 2 and 3 ma;}or shells are used 
in each expansion. Ihus, equation (5.1) becomes, 

. 1-T 

4'a = 51 “Pu 

n=l 


or 








(5.2) 


In a spherical field the angular partsof the single particle 

m 

wave functions are the vector spherical harmonics , -.defined 

in equation (2.14) and thus, more explicitly, the wave functions 
'^and f can be written as 

™ m A r 






(5.3) 


where E and^i/ specify the radial (see expression 3-19) and 
isospin parts of the wave functions respectively. Using 
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expressions (5.5) > equation (5.2) can "be rewritten as. 


E 


HE 


S °n(i3) 


(5.4) 


n=l 


18 18 10 

The low-lying excited states of 0, I and ^Be are 
obtained under the assumption that the last two nucleons move 
in a HE field and interact v/ith each other through some 
residual interaction. The remaining (A-2) nucleons take part 
only in producing the average field. The energy levels are 
obtained by cons tmc ting and diagonalizing a Hamiltonian matrix, 
including all the degenerate configurations having the lowest 
possible energi es.The matrix elements are apJ'T , 

where a, p, y? 6 specify the HP single particle states. These 
matrix elements are evaluated analytically using expansion (5.2), 

<apj<r|Y|Y5JT>jjj, = 772 


*^ 1 ^ 2 ^ 1^2 


( 5 . 5 ) 


Writing explicitly. 


Nl^l^l’ 


n-. 


n. 




^^2^4 


^a'^^ti^i) ^a£^^l^l^ 




(5.6) 
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Ihus, each, of the HP matrix element is expressed as a sum of 
several two-body matrix elements in harmonic oscillator basis 
v/ith proper weight-factors (products of different G-coeff icients) , 
jOhere are eighty-one terms to be sumaed in equation (5«6) vdien 
3 states are used in expansion (5.2). She expansion coef- 
ficients are obtained self-consistently and the results have 
already been presented in Chapter I¥, 2he tv/o-body matrix i 

elements of the interaction in harmonic oscillator basis have * 
been calculated using equation (3.26). A maximum of relative 
1 / = 2 states has been used in our calculations. Single part- 
icle energies are added to the diagonal terms. We have carried | 

I 

out the calculations with both the HP sin^e particle energies ; 
calculated in Chapter and the experimental single particle i 
energies. The latter approach has been followed by most of the ; 
authors (Kuo 66 s Git 68) . It has got the advantage that it : 

takes into account the interaction of the last two neutrons 
with the closed core to all orders, but the disadvantage is 
that it may obscure the shortcomings of our calculations which 
might be obvious otherwise. i 

The following two nucleon-nucleon interactions have 
been used in our calculation; 

(i) Effective Yale interaction (see Section II. 2) 

(ii) Sussex interaction (see Section 11,3) 

The prescription used to obtain the effective Yale interaction 
from the Yale potential (Shn 67a), and the matrix elements of | 
the Sussex interaction (Elt 68a) in an oscillator basis 
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directly from experimental phase shift data ha's heen described ^ 
in Chapter II . 

In an ideal case, the shell model calculations should 
be carried out in an infinite dimensional, configuration space 
when the free nucleon-nuolepn interaction is used as residual 
interaction. Ihis is, however, not practicable and the shell 
model calculations can only be carried out in a limited config- 
uration space. She neglected shells renormalize the interaction j 
in the space of the shells considered in one's calculation, Ihe 
problem is then to find out an effective interaction -sdiich would 
give the same results, when used in a limited configuration i 

space, as the true interaction would give in an infinite dimens- ! 

I 

i 

ional configuration space. Several authors have worked in this 

I 

direction of establishing a relationship between free nucleon- I 

nucleon interaction and effective interaction. ¥e mention here | 

the work of Kuo and Brown (Kuo 66, 67) dealing with Hamada- 

Johnston (HJ) potential, and of Cllement and Baranger (Olt 68) 

dealing v;ith labakin potential. In their calculations of P 
18 

and 0, these authors find that there are two important cor- 
rections arising from neglected configurations, to the free 
nucleon-nucleon interaction. First is the second Bom correction 
in the matrix elements of triplet even tensor force. Ihese are 
the corrections arising from the excitation of two valence 
nucleons to intermediate states of very high-energy. Ihese 
cGrrections mainly affect the ^=0 matrix elements. In deriv- 
ing the effective Yale interaction from the Yale potential, 



Shakin et. al. (31m 67a) have ©alculated the second-order Bom 
corrections by conbining Moszkowski-Scott separation technique 
and Kuo and Brorm netliod (Kuo 66) of taking intermediate states 
as plane waves and angle averaged Pauli operator. ITo such 
corrections are included in Sussex interaction. 

2he second important correction arises due to the 

16 

excitation of the particles from 0 core to the intermediate 
states, which lie close in energy to those included explicitly 
in the diagonalisation. Such corrections are known as core- 
polarization effects. Kuo and Brownf and Clement and Baranger 
have calculated the admixture of 3p-lh (three particle-one hole) 
states to the usual 2p state hy using second-order perturbation 
theory in harmonic oscillator basis. They find that these cor- 
rections, although relatively smaller than the Born cor- 
rection, are quite important for'T' = 1 states particularly for 
1 = 0"^ states. 'T" = 0 states are, however, almost unaffected. 

Ho such corrections have been included in our calculations for 
either of the t\70 interactions used* • 


The configuration space iised in our calcuiations for 
1 '= 18 nuclei consists of (ld ^^2 ’^^ 1 / 2 } model states. 

The ground state configuration for both these nuclei is 
(Id^^g)^. The level spectrum that one would observe in 0 is, 

. 2 *, 4 + 


(ld5/2)^ = 

J = 0"^ 


J = O'^ 


3 2"^ 
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ill these states have total isospin ^ = 1. The follov/ing'T"= 0 
levels (besides the ^ = 1 levels mentioned ahov^ would he 
observed in the case of 


{ld5/2)2 

: J = I'*’, 3"^, 5 '^ 


: J = !■*■ 


; J — 2^ , 3^ . 


The ld^y '2 T?hich lies at 5MeY in ^*^0, has not been 

included in our calculations mainly because of computational 

<3-ifficulties . The HP single particle parameters used in our 

1 R 

calculations of 0 should have been obtained from self- 

17 

consistent HP calculations for. 0, The observed groimd state 

17 ' K+ 

spin of '0 is J = •“ . The HP method described in Chapter III 

17 ' ' 

would only describe an intrinsic state of 0 and the ground 

'5+ 

state with spin J = -^ can be obtained by projecting out good 
angular moment-um states. To avoid such complexities, we decided 
to take single particle parameters from solf-consistent HP 

1 ^ 1 Q 

calculations for 0 or 0. The spherical HP calculations 

give the J = O"*” ground state of these nuclei. In taking the 

16 

single particle parameters from 0 liP calculations the contrib- 
utiore of the interaction of outer two neutrons with the closed 

core to the self-consistent average field and wave functions 

IB' 

are missed. In 0 HP calculations, part of the interaction 
between the last two neutrons is already included in the calcul- 
ation of self-consistent field and thus, the residual inter- 
action left over is not just the free nucleon-nucleon inter- 
action but gets modified. However, we think that these 
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corrections would de small and average self-consistent field 
would not be very different for ^^0, ^*^0 and In all our 

calculations the last two neutrons in ^^0, and one neutron and 
one proton in have been treated on eq^ual footing. It is 

assumed that the self-consistent field is the same for neutron 
and proton so t -IC- c same wave fimctions (neutron wave functions) 
describe the notion of both the particles outside the closed 
core. This is, however, not a very serious restriction as v/e 
have seen in Giiapter lY that coulomb force has little effects 
on the single particle HF wave functions of light nuclei. 


10 


The configuration space used for the calculations of 
Be levels consists of lp,/o and lp-> shell model states. 


' 3/2 -"^ 1/2 

O 

The ground state configuration for this nucleus is (Ip^y^) , 

10 rr' 

All low-lying levels of Be would have X = 1 and the following 


spectamm would be observed 


(lP5/2)^ 

4* 

0 

II 

(lp^/2 lPi/2^ 

: 1 = 1"^ 

(IP 1 / 2 )" 

: J = O'*' 


^•2 Results and Discussion 

Results of our calculations obtained with effective 

' ' 18 

Yale interaction are shown in Figures (Y.l - V.2) for 0 and 
respectively. In these figures, YA and YB refer to calcul- 
ations done with single particle wave functiffas obtained from 
HF calculations for 0 without coulomb and centre-of-mass 



FIG-IZ'-I Poeit;tV'3 jparity 0 with wffeatiif 

lale intcxaotip^ aiift HF ^-wawa .f 'Wh.tions » All otataa 
have T *= 




YA YB YC YD YE YF, EXPT- 


FiG.2‘2 Poaiti're pa3?i% atatta-sf Ototadaed wUSa affeotiv® 
Yale iateraotioa aad IfiP wars fuaotioaa. I«®velB. wilSi 
T *s 1 are epeoifiad adilto braefcated tpia. All oldiaa 
larals hme T ■* 0* 
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corrections, and using HP and experimental single particle 
energies respectively (see Tables IT.llCb) and IY.5). Cases YO 
and YD refer to calculations done using uncorrected (witiiout 


coulomb and centre-of-mass corrections) IIP wave functions for 


18 . 


‘0 


and t/itli fiP and experimental single particle energies respectively 
(see Tables I?, 13(a) and I?, 5). Gases YE and YF refer to calcul- 
ations done using experimental single particle energies»and v/ave 
fimctions obtained from HP calculations ?/itb coulomb and centre- 


of-fflass corrections for the nuclei ^^0 and ^^0 respectively 


(see Tables lY. 11(a) and IV, 13(a)). All calculations witb 
effective Yale interaction have been carried out for b = 2»09P 
since for this value of b (=Vh/m^d)) we could obtain minimvim 
total energy for 0 in self-consistent HP calculations (see 
Figure lY.l) . For consistency, calculations for Sussex inter- 
action have been done for b = 2. OF. The experimental single 
particle energy difference between and 2s^y2 state was 

taken to be 0.87 licY, as obtained from experimental levels of 

1 7 1 .R 

'0. This is not justified for F for v/hich an average of the 

neutron and proton single particle energies for 2s^y2 s^^a'bes 
should be taken. However, the discrepancy would be small. The 


no 1 R 

experimental levels for 0 and F, shorn in Figures Y.l and 


Y. 2, have been talieh from references (Hea 64) and (Poi 65) 
respectively. The absolute ground state energies are deterrained 
as follows. Using the binding energy data of Koenig et. al, 
(Kog 62) , the ground state energies of valence nucleons are 
taken to be 
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-B.E.(^®0) - B.E.(^^O) + 2 B.E.(^'^O) = -3.90Me? for 
and 

-B.E.(^®F) - + B.E.(^'^O) + 3.E.(^'^B) = -5-00 MeY 

for . 

As seen iron the figures, except for the 4"^ ^=1^ states 
the calculated levels are in very poor agreement with experim- 
ental results in all the cases. In particular, the observed 
ground state of '^1' is a 1"^ state whereas in all our calcul- 
ations we get 5^ state as lowest state and first 1"^ state lies 
much higher. However, the fact that very similar results are 
obtained in various cases supports our assumption that HI field, 
if constrained to remain spherically symmetric, would not be 
very different for ^^0, ^"^0 and ^^0. Results obtained with 
effective Tale interaction using pure oscillator wave functions 
and experimental single particle energies are shown in Figures 
(Y,3 - Y.4) foxC =1 and'C'=0 states respectively. In these 
figures, HOI and 1102 refer to calculations done with harmonic 
oscillator (HO) wave functions in configuration space consisting 

of ( 145/2 ^®l/2^ ^^^5/2 ^®l/2 ^^ 3 / 2 ^ states respectively. 

The importance of including the Id^/g configuration is obvious 
from the figures. It lowers theC=l levels, thereby improving 
the agreement with experiments; and also the first 1"^ level in 

is pushed below the lowest 5"^ state. In Figures (Y.3 - Y,4) , 
we also compare the calculated C=1 and”0 levels obtained with 
HF wave functions in case YB with the results of Kuo and Brovai 
(referred to as KB), and Clement and Baranger (referred to as 



ENERGY CMeV ) 


HF(YB) HOI H02 


CB2 


FlG'2-3 Comparison of calculated HF and HO positive 

parity ( T » 1) levels of for effective Yal 
interaction witab, the resnlts of Kuo and Brown 
'(EB), and" Clement and Barai^er ( results ^ with 
second-order Bom correction are referred to 
as CB2).* 




FlG.2:-4 aoraparisoa of oai,oiaat^:Hf and 10 posltiT®^ parity 

levels of ^%_for effactive Tale interaotioa with, the 
results of 'Kuo"aad' f and Oleaent ' aiid Bahaages 
(0ii2) . Levels with hiftaohated spin have T *, !• AIL 
ptlsi levels have T 


GB2 when second-ordex Bom term is include d) . Both these calcul- 
ations were done in a configuration space consistirxg of 
^^1/2 ^^3/2^ h.£.rconic oscillator sta,tes and using experimental 
single particle energies. KB and G32 levels shorn here include 
the oontribution of second-order Bom term but not the core- 
polarization effects. On comparing the HF spectrum obtained in 
case YB and the JIOl spectrum (see Figures Y.3 “ Y.4) we find 
that other parameters being held fixed, Iff and HO wave functions 
do not give very different results. 

In Sables ?.l and ¥.2, we present the calculated two- 
body HO and Iff matrix elements in the s-d shell obtained in 
various cases. Results of several other calculations are also 
given for comparison , Federman (Fen 65) and Arima (Ara 68) 

^=1 matrix elements have been obtained from best fit to expe rim- 
ental levels of 0. Arima T=0 matrix elements are obtained 
from best fit to experimental levels of ^®F and ^^le. She 

values obtained in two different cases, first when all but the 
1.7Me¥, 1"^ state in ^^Fwere included from the fit and second 
when this state was also included are referred to as A1 and A2 
in Sable ¥.2. It can be seen that no impressive agreement of 
the calculated Iff matrix elements exists ¥7ith any one of these 
calculations. In particular, the J=0'*’ HF matrix elements are 
considerably smaller than the results of other calculations. 
However, Federman-Salmi (FS) and Arima have obtained positive 
values for J=3’*’ and i"*" ^=l) matrix elements. In our calculations 
j=3'*' ^'^=1] HF matrix element is ailways positive althou^ much 
smaller in magnitude. 
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The results obtained with Sussex interaction are shown 
in Figures (7.5 - 7.6) for and respectively along ?n.th 
the experinenral levels. In these figures, SA and SB refer to 
calculations done using uncorrected ITP wave functions for ^^0 
(see Table lY. 27(b) ),. and HP single particle (see Table IV. 24) 
and experimen-cal single particle energies respectively. HOI and 
H02 are the results obtained with pure oscillator wave functions 
in a smaller (id^y'g 282 ^^ 2 ) and larger (ld^y 2 ^^ 1/2 ^^3/2^ 
configuration spa.ce respectively. Clement and Saranger results 
obtained ¥/ith bare Tab akin interaction without second-order 
corrections (referred to as CBl) are also given for oomparisoia. 
In Table V.3, calculated two-body HP matrix elements for Sussex 
interaction in the s-d shell have been compared with Clement 
and Baranger results for Tabakin interaction (CBl), 


10 

The calculated energy levels of Be for effective Yale 
interaction are shown in Figure V.7 along with the experimental 
data. In this figure, XI and XB refer to calculations done 
using HP wave fuactions for ^^Be (see Table IV,12), and HP single 
particle (see Table IV. 5) and experimental single particle 
energies respectively. XG is the result of calculation using 
HP wave functions of ^e (see Table IV. 9) and experimental single 
particle energies. Again, HOI and H02 are the results of calcul- 
ations using oscillator wave functions with experimental single 
particle energies in configuration space consisting of (Ip^^r^ 
Ip^^g) ^^^3/2 ^^ 1/2 ^*^5/2^ states respectively. The experim- 
ental energy separation between ^^1/2 taken 



FlG*]x*5 OoBSarison of calcmla’ted Hf and fl© poaitiw 

levels of' 6 for' Su®»«r 'irteraotioa with the results 
of Clement aad Baraager.freiRilts without seotaid-order 


of Clement and 3arsD^fr.JrfiWp-t»_^,wlthout »eo<md-order 
Bom correction are to',a®,CBl), and experiment: 
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to be Site Y. All calculations liave been done for b=2.091*» We 
find that the calculated difference between the ground state 
and first 2”^ state is much less compared to experimental 
separation. Hov<ever, the calculated energy separation between 
the first and the second 2 '*' states is 5.89 Me¥ in case XB 
compared to experimental separation of 5.9 MeV. 

In all our calculations we find that calculated HP 
spectra are compressed compared to experimental levels. She 
poor agreement with the experimental results may be due to the 
nature of the nucleon-nucleon interactions used in oux calcul- 
ations. As seen in Chapter IV, both these interactions give 
3-4MeV less binding when used in spherical HP calcvilations. 
However, spin-orbit splittings in various nuclei have been 
satisfactorily reproduced with both these interactions (Shn 67b: 
Elt 67). fhus, it appears that the spin-orbit part of these 
interactions is reasonable but the central part is weah. In 
yet another entirely different approach, Gunye has calculated 
(Gue 68b, 69) the low-lying states of various odd and even mass 
light nuclei by projecting out good angular momentum states 
from an intrinsic deformed HP state. He has carried out the 
calculations for 2s-ld shell nuclei (Gue 68b) with effective 
lale interaction# and for Ip shell nuclei (Gue 69) with both 
effective lale and Sussex interactions. In all the cases the 
pro jected HP spectra for even mass nuclei are found to be compr- 
essed compared to the experimental levels, Por example, in the 
case of ^^0, the calculated separation between the ground state 
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to be 6iue?. All calculations have been done for b=2.09P. We 
find that the calculated difference between the ground state 
and first 2”^ state is much less compared to experimental 
separation. Ko’,;ever, the calculated eziergy separation between 
the first and the second 2'‘' states is 5.89 Me¥ in case XB 
compared to experimental separation of 5.9 KeV. 

In all our calculations we find that calculated HI 
spectra are compressed compared to experimental levels. She 
poor agreement vdth the experimental results may be due to the 
nature of the nucleon-nucleon interactions used in our calcul- 
ations. As seen in Chapter I¥, both these interactions give 
3-4Me¥ less binding when used in spherical HP calculations. 
However, spin-orbit splittings in various nuclei have been 
satisfactorily reproduced with both these interactions (Shn 67b 
Elt 67). Ihus, it appears that the spin-orbit part of these 
interactions is reasonable but the central part is weah. In 
yet another entirely different approach, Gimye has calculated 
(Gue 68b, 69) the low-lying states of various odd and even mass 
light nuclei by projecting out good angular momentum states 
from an intrinsic deformed HP state. He has carried out the 
calculations for 2s-ld shell nuclei (Gue 68b) with effective 
Yale interaction# and for Ip shell nuclei (Gue 69) with both 
effective Yale and Sussex interactions. In all the cases the 
projected HP spectra for even mass nuclei are found to be compr- 
essed compared to the experimental levels. Por example, in the 
case of the calculated separation between the ground state 



^ X 

and first 2 state is 0,51 MeV compared to the experimental 
value of 1.93 MeV. 2he effect of increasing state matrix 
elements arbitrarily by over 30 percent so as to reproduce the 
experimental binding energies of various nuclei is to spread 
out the projected spectra. Since the results obtained by 
projection method are close to shell model results (Reh 58), 
the compressed spectrum obtained in the present shell model 
calculations in rD? basis seems to be consistent with the find- 
ings of Gimye. 

It should be emphasized that the poor agreement of our 
results with the e:iperimental data does not imply the shortcomings 
ings in the HR basis since vie have not made any attempt to 
calculate the important second-order corrections in both the 
interactions arising from core excitations (besides the Born 
term for Sussex interaction). To the best of our hnowledge, no 
attempt till now lias been made to use self-consistent single 
particle energies and wave functions obtained with realistic 
potentials in the shell model oalculations > Kuo and Brown 
(Kuo 66) have calculated single particle energies in 0 with 
their realistic potential and obtained q^ualitative agreement 
with observed energies, but they did not use calculated single 
particle energies in obtaining the spectra of A=18 nuclei. 

Our calculation is a first step towards a more complete shell 
model calculation in the Hy basis. In fact, it has been shown 

TO " 

tliat iQw-lying states of 0 are abundant in deformed components 
(Bron 6& End 65) and the assumption of two nucleons outside 
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a closed '^'^0 core- is not very good, liie transition rate from 
second 0"^ to first 2"^ state is also very large and suggests 
that these states should be mainly composed of deformed comp- 
onents. However, Kuo and Brown (Kuo 66) iiave satisfactorily 
reproduced the energy levels of these states by using simple 
two particle model with p-h corrections. It is expected that 
such core excitation corrections ?/ill spread apart the calcul- 
atedTr=l HB spectrum, by depressing the lov/er states and raising 
the higher levels and thus, the agreement with the experiment 
would improve. Studies of Kuo and Brown, and Clement and 
Baranger show that in the harmonic oscillator basis most imp- 
ortant of these corrections are 3p-lh corrections arising from 
excitation energies of Contributions of higher levels, 

which involve excitation energies are relatively unimp- 

ortant. However, this may not be true in BB basis. Bassichis 
et. al. (Bas 69) have made explicit calculations for the secona- 
order energies and probabilities of 2p-2h components in the 
ground state of ile and 0 using perturbation theory in HB 
basis. They find that the relative phases (signs of expansion 
coefficients) of the HB wave functions of unoccupied states 
play an important role in deciding the importance of a part- 
icular level. It is quite likely that excitations to some of 
the high-lying levels may be more important than to low-lying 
levels contrary to the predictions of simple shell-model 
picture in harmonic oscillator representation. However, no 
conclusive remarks can be made about the magnitude of such 
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corrections in HP 130313 unless explicit calculations are 
made . 

In anotlier calculation carried out for the 
residvial interaction "between the last two neutrons in 
was tal'en to be of the form (a + bOj^*^ 2 ) exp(-r^/rQ), Two- 
body matrix elements in the s-d shell were calculated in a 
basis of HP wave functions and vm tried to determine the 
parameters of the potential which would give a good fit to 
the experimental energy levels of ^^0. But no set of 
parameters could satisfactorily reproduce the experimental 
levels of Thus, it appears that a simple residual 

interaction of the form described above is not suitable 
for a HP basis. 



APPEIIDIX 


Tile local potential; used in our calculations 
single particle energies, is of the following form. 






where 


(A.l) 


f(r) 


= i 1 + exp(-^) Y , 

[ ] 


SOn J 


g(r) 


) Si + exp(- 


lx(r) = i 


for rNIU 


A-(3 


for r E, 


(A. 2) 


Y,-r = Central nuclear potential 


Y = Spin-orhit nuclear potential; 
so > 


dY), 


for 3 = •6+'^- 


(t+1) for 3 = 


(A. 3) 


(xiii) 



^x’ ^x radius and diffuseness parameters of the x"^*" 

part of the potential and R = r 

X X 

^A mass and its charge, 

p, are the particle reduced mass and its charge. 

The radial v;ave functions "the solutions of the 

following Sclirodinger equation, 




The hinding energies of the single particle levels are fed in 
and the strength parameters and Y^^ are obtained by the 
method of successive iterations till the values converge. 

The wave functions are normalized to uni1y. 


A separable form of the local potential is given by, 
Y(r)q;(?) = jY(-^i^) H(| r-^^ t)4^' (?)dr^ (A.5) 


H(j~r-f*^|) is given the Gaussian form, normalized to unity. 


I) 



(a. 6) 


where is the non-locality range parameter. An approximate 
relationship between the wave function calculated in a 

potential of the form, given in equation A.5, and the local 
wave function ^ is given by, . 

u^j(r) = exp(-a^Vj(r)/8)iiJy{r) , where 

V^Cr) =%Tnf(r)+c|^)^^?:3)g(r)-Z^2 ®2h(r)] fA.7) 

J di It so 
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APPEilDIX 3 


It has oeen mentioned before that the t’,vo-body 
coiiloiiib and centre-of-xnass matrix elements used in the HP 
calculations v;itli effective Yale interaction were calculated 
with a maximum relative ?. (see Section IY,1) . This has 
caused some error in our calculations since higher f/ values 
are also quite important for both of these matrix elements. 
•The new coulomb and centre-of-mass matrix elements in the 
configuration space consisting of Sqy'2’ ^3/2 
oscillator states of n = 1, 2 and 3 major shells were 
calculated with maximum relative i =10 . The revised HP 
results obtained by using these new matrix elements are 
shown in Table B.l for the nuclei ^He, ^Be, and ^^0. 

On comparing with the older results (see Table I¥«l) we 
find that the effect of truncation at ^=2 on most of the 
static properties of closed shell nuclei is small. However, 
the wave functions are affected, mainly in the tail region. 




fable B.l 

Iiesults of revised HF calculations with effective Yale interaction 

for various nuclei 
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APPENDIX G 


It is ass-’amed that Sussex interaction V can be broken 
up into centraJ., spin-orbit and texisor parts v/hich are tensors 
of rank zero^ one and tv/o respectiTely in the scalar product 
) . Here refers to the spin variables and 

to the orbita.l motion. Ihus, we can write, 

k=0 


t— 

y. y(k) 


I :=0 


= y(0)+Y(l)+V(2) 


and the matrix elements of V(k) are 


( 0 . 1 ) 


<(n' S fjj jY(l') '> = <^n*S |nStj]> 


(- 1 ) 


J+S+ 


s t J' 




I 


S k! 


<s II II s>X 


X <n^ ll II n I >. ( 0 . 2 ) 


SuiiiM.iiig .oTC.x on "bolili ’tli'8 sides we obtain,. 


^ j/ fs t j 

<n' 5fc'j|7|nSf,J > ='^(-1 ) i 

k ^ 


S k 


<s II s> X 
X <n't' |(T^‘^|hl>- 
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Multiplying both the sides by (-D^+S+t yijy i < 

1 ^ S k I 

aud sumsiiiiy over J we get 



Thus, in equation (C.4) the uncoupled matrix elements of any 
particular component are expressed as a sum of matrix element 
of the total potential ¥ in different channels. The coupled 
matrix elements of ¥(lc:) are obtained by substituting equation 
(G.4) in equa,tion (G,2) . 

Tile off-diagonal matrix elements of the Sussex 
interaction ¥, with radial quantum numbers differing by 
more than one, have been calculated using the folloY/ing 
expression, given by Elliott et. al. {Elt 68a), 


(xviii) 



j ¥ |n J 2(r/ -n-1) n'' -1 f¥j n y +^/ (n+l) (n+^i+l^ <^n-l |7| 


f i/n(n+ 1 +|-) <;^n^ -1 j¥| n-1 > -/(n' -l) ( rJ +6 ) y 


X < n-2 |Vj n >+<^ n-1 j [r^,v]/2t^j n %jj/nin+ i+|^ 


(G.5) 


Ihe relation (0.5) connects the matrix elements with 
n^ -n = 2 tc those with X -n4 1 and is true for the matrix 
elements of any general interaction Y in an oscillator basis. 
It is assumed that interaction V is raomentum independent so 

p 

that comnuts/uor [r , Y] vanishes. 
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